Data processing of the Kuiyang-ST2000 deep-towed high-resolution multichannel seismic system and application to South China Sea data

The Kuiyang-ST2000 deep-towed high-resolution multichannel seismic system was designed by the First Institute of Oceanography,Ministry of Natural Resources(FIO,MNR).The system is mainly composed of a plasma spark source(source level:216 dB,main frequency:750 Hz,frequency bandwidth:150-1200 Hz)and a towed hydrophone streamer with 48 channels.Because the source and the towed hydrophone streamer are constantly moving according to the towing configuration,the accurate positioning of the towing hydrophone array and the moveout correction of deep-towed multichannel seismic data processing before imaging are challenging.Initially,according to the characteristics of the system and the towing streamer shape in deep water,travel-time positioning method was used to construct the hydrophone streamer shape,and the results were corrected by using the polynomial curve fitting method.Then,a new data-processing workflow for Kuiyang-ST2000 system data was introduced,mainly including float datum setting,residual static correction,phase-based moveout correction,which allows the imaging algorithms of conventional marine seismic data processing to extend to deep-towed seismic data.We successfully applied the Kuiyang-ST2000 system and methodology of data processing to a gas hydrate survey of the Qiongdongnan and Shenhu areas in the South China Sea,and the results show that the profile has very high vertical and lateral resolutions(0.5 m and 8 m,respectively),which can provide full and accurate details of gas hydrate-related and geohazard sedimentary and structural features in the South China Sea.

Accurate Source-Receiver Positioning Method for a High-Resolution Deep-Towed Multichannel Seismic Exploration System

The near-seabed multichannel seismic exploration systems have yielded remarkable successes in marine geological disaster assessment,marine gas hydrate investigation,and deep-sea mineral exploration owing to their high vertical and horizontal resolution.However,the quality of deep-towed seismic imaging hinges on accurate source-receiver positioning information.In light of existing technical problems,we propose a novel array geometry inversion method tailored for high-resolution deep-towed multichannel seismic exploration systems.This method is independent of the attitude and depth sensors along a deep-towed seismic streamer,accounting for variations in seawater velocity and seabed slope angle.Our approach decomposes the towed line array into multiline segments and characterizes its geometric shape using the line segment distance and pitch angle.Introducing optimization parameters for seawater velocity and seabed slope angle,we establish an objective function based on the model,yielding results that align with objective reality.Employing the particle swarm optimization algorithm enables synchronous acquisition of optimized inversion results for array geometry and seawater velocity.Experimental validation using theoretical models and practical data verifies that our approach effectively enhances source and receiver positioning inversion accuracy.The algorithm exhibits robust stability and reliability,addressing uncertainties in seismic traveltime picking and complex seabed topography conditions.

High-resolution crustal structure of the Yinchuan basin revealed by deep seismic reflection profiling:implications for deep processes of basin

The Yinchuan basin, located on the western margin of the Ordos block, has the characteristics of an active continental rift. A NW-striking deep seismic reflection profile across the center of Yinchuan basin precisely revealed the fine structure of the crust. The images showed that the crust in the Yinchuan basin was characterized by vertical stratifications along a detachment located at a two-way travel time（TWT） of 8.0 s.The most outstanding feature of this seismic profile was the almost flat Mohorovicˇic′ discontinuity（Moho） and a high-reflection zone in the lower crust. This sub-horizontal Moho conflicts with the general assumption of an uplifted Moho under sedimentary basins and continental rifts, and may indicate the action of different processes at depth during the evolution of sedimentary basins or rifts.We present a possible interpretation of these deep processes and the sub-horizontal Moho. The high-reflection zone, which consists of sheets of high-density, mantlederived materials, may have compensated for crustal thinning in the Yinchuan basin, leading to the formation of a sub-horizontal Moho. These high-density materials may have been emplaced by underplating with mantlesourced magma.

High-resolution crustal velocity imaging using ambient noise recordings from a high-density seismic array:An example from the Shangrao section of the Xinjiang basin,China

A profile of shallow crustal velocity structure(1–2 km) may greatly enhance interpretation of the sedimentary environment and shallow tectonic deformation.Recent advances in surface wave tomography, using ambient noise data recorded with high-density seismic arrays, have improved the understanding of regional crustal structure. As the interest in detailed shallow crustal structure imaging has increased, dense seismic array methods have become increasingly efficient. This study used a high-density seismic array deployed in the Xinjiang basin in southeastern China, to record seismic data, which was then processed with the ambient noise tomography method. The high-density seismic array contained 203 short-period seismometers, spaced at short intervals(～ 400 m). The array collected continuous records of ambient noise for 32 days. Data preprocessing,cross correlation calculation, and Rayleigh surface wave phase-velocity dispersion curve extraction, yielded more than 16,000 Rayleigh surface wave phase-velocity dispersion curves, which were then analyzed using the direct-inversion method. Checkerboard tests indicate that the shear wave velocity is recovered in the study area, at depths of 0–1.4 km,with a lateral image resolution of ～ 400 m. Model test results show that the seismic array effectively images a 50 m thick slab at a depth of 0–300 m, a 150 m thick anomalous body at a depth of 300–600 m, and a 400 m thick anomalous body at a depth of 0.6–1.4 km. The shear wave velocity profile reveals features very similar to those detected by a deep seismic reflection profile across the study area. This demonstrates that analysis of shallow crustal velocity structure provides high-resolution imaging of crustal features.Thus, ambient noise tomography with a high-density seismic array may play an important role in imaging shallow crustal structure.

Purposeless repeated acquisition time-lapse seismic data processing

In China, most oil fields are continental sedimentation with strong heterogeneity, which on one side makes the reservoir prospecting and development more difficult, but on the other side provides more space for searching residual oil in matured fields. Time-lapse seismic reservoir monitoring technique is one of most important techniques to define residual oil distribution. According to the demand for and development of time-lapse seismic reservoir monitoring in China, purposeless repeated acquisition time-lapse seismic data processing was studied. The four key steps in purposeless repeated acquisition time-lapse seismic data processing, including amplitude-preserved processing with relative consistency, rebinning, match filtering and difference calculation, were analyzed by combining theory and real seismic data processing. Meanwhile, quality control during real time-lapse seismic processing was emphasized.

A novel technique for automatic seismic data processing using both integral and local feature of seismograms

A novel technique for automatic seismic data processing using both integral and local feature of seismograms was presented in this paper. Here, the term integral feature of seismograms refers to feature which may depict the shape of the whole seismograms. However, unlike some previous efforts which completely abandon the DIAL approach, i.e., signal detection, phase identifi- cation, association, and event localization, and seek to use envelope cross-correlation to detect seismic events directly, our technique keeps following the DIAL approach, but in addition to detect signals corresponding to individual seismic phases, it also detects continuous wave-trains and explores their feature for phase-type identification and signal association. More concrete ideas about how to define wave-trains and combine them with various detections, as well as how to measure and utilize their feature in the seismic data processing were expatiated in the paper. This approach has been applied to the routine data processing by us for years, and test results for a 16 days＇ period using data from the Xinjiang seismic station network were presented. The automatic processing results have fairly low false and missed event rate simultaneously, showing that the new technique has good application prospects for improvement of the automatic seismic data processing.

Optimization of a precise integration method for seismic modeling based on graphic processing unit

General purpose graphic processing unit （GPU） calculation technology is gradually widely used in various fields. Its mode of single instruction, multiple threads is capable of seismic numerical simulation which has a huge quantity of data and calculation steps. In this study, we introduce a GPU-based parallel calculation method of a precise integration method （PIM） for seismic forward modeling. Compared with CPU single-core calculation, GPU parallel calculating perfectly keeps the features of PIM, which has small bandwidth, high accuracy and capability of modeling complex substructures, and GPU calculation brings high computational efficiency, which means that high-performing GPU parallel calculation can make seismic forward modeling closer to real seismic records.

A systematic printability study of direct ink writing towards high-resolution rapid manufacturing

Direct ink writing(DIW)holds enormous potential in fabricating multiscale and multi-functional architectures by virtue of its wide range of printable materials,simple operation,and ease of rapid prototyping.Although it is well known that ink rheology and processing parameters have a direct impact on the resolution and shape of the printed objects,the underlying mechanisms of these key factors on the printability and quality of DIW technique remain poorly understood.To tackle this issue,we systematically analyzed the printability and quality through extrusion mechanism modeling and experimental validating.Hybrid non-Newtonian fluid inks were first prepared,and their rheological properties were measured.Then,finite element analysis of the whole DIW process was conducted to reveal the flow dynamics of these inks.The obtained optimal process parameters(ink rheology,applied pressure,printing speed,etc)were also validated by experiments where high-resolution(<100μm)patterns were fabricated rapidly(>70 mm s^(-1)).Finally,as a process research demonstration,we printed a series of microstructures and circuit systems with hybrid inks and silver inks,showing the suitability of the printable process parameters.This study provides a strong quantitative illustration of the use of DIW for the high-speed preparation of high-resolution,high-precision samples.

Online Data Processing Methods in Xlet Applications for Seismic Early Warning and Emergency Services on Interactive Digital Television

In this paper, I described the methods that I used for the creation of Xlets, which are Java applets that are developed for the IDTV environment;and the methods for online data retrieval and processing that I utilized in these Xlets. The themes that I chose for the Xlets of the IDTV applications are Earthquake and Tsunami Early Warning;Recent Seismic Activity Report;and Emergency Services. The online data regarding the Recent Seismic Activity Report application are provided by the Kandilli Observatory and Earthquake Research Institute (KOERI) of Bogazici University in Istanbul;while the online data for the Earthquake and Tsunami Early Warning and the Emergency Services applications are provided by the Godoro website which I used for storing (and retrieving by the Xlets) the earthquake and tsunami early warning simulation data, and the DVB network subscriber data (such as name and address information) for utilizing in the Emergency Services (Police, Ambulance and Fire Department) application. I have focused on the methodologies to use digital television as an efficient medium to convey timely and useful information regarding seismic warning data to the public, which forms the main research topic of this paper.

On electro-acoustic characteristics of a marine broadband sparker for seismic exploration

The construction of major marine infrastructure projects and the exploration and development of deep-sea mineral resources require fine imaging of seabed strata and structures.The highresolution marine seismic exploration based on a high broadband sparker source is an important approach to reveal seabed stratum and reservoir structure,and identify geohazard.To optimize the performance of sparker seismic source,we investigated the electro-acoustic characteristics of spark discharge under conditions of different charging voltages and electrode numbers.Results show that the sound source level increased with the increase of the charging voltage,whereas the main frequency decreased when the charging voltage increases.In addition,it was found that the charging capacitance had more obvious influence on the main frequency than the sound source level did.Although the load energy decreased with increasing electrode number,the sound source level still increased but the main frequency decreased.Meanwhile,the primary to bubble(P/B)ratio increased with the increase of the electrode number.To gain a deeper insight into the electro-acoustic characteristics,we investigate the relationship between sound source level and power peak,from which a good correlation was observed.A more practical statistical analysis on the rise rate of current was processed,and a perfect logarithmic function was derived.Furthermore,we found that the main frequency was most possibly subjected to the electrical energy,especially the charging energy per electrode.The results indicate that the charging energy per electrode less than 10 J could increase the main frequency to above 300 Hz.At last,the main frequency could be reduced to 20 Hz when the charging energy of a single-electrode discharge was enhanced to over 4 kJ.This study shall be helpful in developing a sparker seismic source and improving the performance for marine engineering exploration and geohazard assessment.

Research on Seismic Data Processing and Interpretation System Model Based on CSCW

Through analyzing the needs of seismic data processing and interpretation,a system model based on CSCW is designed.Using the technology of CSCW to build the environment of cooperative work allows the field data acquisition to possess the functions of remote real-time guidance by experts and remote real-time processing of the data.The model overcomes the influences and barriers existing in the areas

Pre-processing and post-processing method for geostress simulation using seismic interpretation results

In order to carry out numerical simulation using geologic structural data obtained from Landmark(seismic interpretation system), underground geological structures are abstracted into mechanical models which can reflect actual situations and facilitate their computation and analyses.Given the importance of model building, further processing methods about traditional seismic interpretation results from Landmark should be studied and the processed result can then be directly used in numerical simulation computations.Through this data conversion procedure, Landmark and FLAC(the international general stress software) are seamlessly connected.Thus, the format conversion between the two systems and the pre-and post-processing in simulation computation is realized.A practical application indicates that this method has many advantages such as simple operation, high accuracy of the element subdivision and high speed, which may definitely satisfy the actual needs of floor grid cutting.

Seismic Data Processing Approaches for the Study of Gas Hydrates in the East China Sea

A comprehensive study of the data profiles, including the 2D seismic data, single channel seismic data, shallow sections, etc., reveals that gas hydrates occur in the East China Sea. A series of special techniques are used in the processing of seismic data, which include enhancing the accuracy of velocity analysis and resolution, estimating the wavelet, suppressing the multiple, preserving the relative amplitude, using the DMO and AVO techniques and some special techniques in dealing with the wave impedance. The existence of gas hydrates is reflected in the subbottom profiles in the appearance of BSRs, amplitude anomalies, velocity anomalies and AVO anomalies, etc. Hence the gas hydrates can be identified and predicted. It is pointed out that the East China Sea is a favorable area of the gas hydrates resource, and the Okinawa Trough is a target area of gas hydrates reservoir.

A real time processing system of seismic waves using personal computers-Function and characteristics

Seismological Bureau of Sichuan Province, Chengdu 610041, China2) Center for Analysis and Prediction, State Seismological Bureau, Beijing 100036, China3) Observation Center for Prediction of Earthquakes and Volcanic Eruptions, Faculty of Sciences, Tohoku University, Sendai 98077, Japan

Study on Anti-Disturbance and High-Resolution Shallow Seismic Exploration of Active Faults in Urban Regions

The significance of detection of urban active faults and the general situation concerning detection of urban active faults in the world are briefly introduced. In a brief description of the basic principles of anti-disturbance and high-resolution shallow seismic exploration, the stress is put on the excitation of seismic sources, the performance of digital seismographs, receiving mode and conditions, geometry as well as data acquisition, processing and interpretation in the anti-disturbance and high-resolution shallow seismic exploration of urban active faults. The study indicates that a controlled seismic source with a linear or nonlinear frequency-conversion scanning function and the relevant seismographs must be used in data acquisition, as well as working methods for small group interval, small offset, multi-channel receiving, short-array and high-frequency detectors for receiving are used. Attention should be paid to the application of techniques for static correction of refraction, noise suppressing, high-precision analysis of velocity, wavelet compressing, zero-phasing of wavelet and pre-stacking migration to data processing and interpretation. Finally, some cases of anti-disturbance and high-resolution shallow seismic exploration of urban active faults are present in the paper.

Wavelet analysis and its application to signal processing

The construction of basic wavelet was discussed and many basic analyzing wavelets was compared. Acomplex analyzing wavelet which is continuous, smoothing, orthogonal and exponential decreasing was presented, andit was used to decompose two blasting seismic signals with the continuous wavelet transforms (CWT). The resultshows that wavelet analysis is the better method to help us determine the essential factors which create damage effectsthan Fourier analysis.

Seismic risk evaluation for a planning mountain tunnel using improved analytical hierarchy process based on extension theory

Seismic risk evaluation(SRE) in early stages(e.g., project planning and preliminary design)for a mountain tunnel located in seismic areas has the same importance as that in final stages(e.g.,performance-based design, structural analysis, and optimization). SRE for planning mountain tunnels bridges the gap between the planning on the macro level and the design/analysis on the micro level regarding the risk management of infrastructural systems. A transition from subjective or qualitative description to objective or quantitative quantification of seismic risk is aimed to improve the seismic behavior of the mountain tunnel and thus reduce the associated seismic risk. A new method of systematic SRE for the planning mountain tunnel was presented herein. The method employs extension theory(ET)and an ET-based improved analytical hierarchy process. Additionally, a new risk-classification criterion is proposed to classify and quantify the seismic risk for a planning mountain tunnel. This SRE method is applied to a mountain tunnel in southwest China, using the extension model based on matter element theory and dependent function operation.The reasonability and flexibility of the SRE method for application to the mountain tunnel are illustrated.According to different seismic risk levels and classification criteria, methods and measures for improving the seismic design are proposed, which can reduce the seismic risk and provide a frame of reference for elaborate seismic design.

Automatic Extraction of Urban Road Centerlines from High-Resolution Satellite Imagery Using Automatic Thresholding and Morphological Operation Method

The commercial high-resolution imaging satellite with 1 m spatial resolution IKONOS is an important data source of information for urban planning and geographical information system (GIS) applications. In this paper, a morphological method is proposed. The proposed method combines the automatic thresholding and morphological operation techniques to extract the road centerline of the urban environment. This method intends to solve urban road centerline problems, vehicle, vegetation, building etc. Based on this morphological method, an object extractor is designed to extract road networks from highly remote sensing images. Some filters are applied in this experiment such as line reconstruction and region filling techniques to connect the disconnected road segments and remove the small redundant. Finally, the thinning algorithm is used to extract the road centerline. Experiments have been conducted on a high-resolution IKONOS and QuickBird images showing the efficiency of the proposed method.

Application of digital image processing to the determination of spatial distribution of earthquakes

This paper has proposed a practical method for determining the spatial distribution pattern of earthquakes by means of digital image processing and has given some calculation results. This method can overcome the artificial arbitrariness which is usually inevitable in determining the spatial distribution of earthquakes. Meanwhile, the form of seismic gaps and the method for identifying seismic gaps have also been discussed. It should be pointed out that the method proposed in this paper is a new idea in this respect. However, this method is still unable to determine the seismic gap uniquely. In order to identify the real seismic gap, comprehensive analysis of the variation of other relevant parameters (e. g., the b-value, etc. ) should also be made. The results are as follows; Thephenomenon that an area of high seismic activity is surrounded by low seismicity areas, i. e., the seismic clustering pattern is the general form of spatial distribution of earthquakes; whereas a low seismicity area surrounded by high seismicity areas, i. e., the seismic gap pattern, is only an important form of spatial distribution of earthquakes. Earthquakes often occur in a certain seismic clustering area or on its margin.

DEVELOPMENT OF FAST TWO-DIMENSIONAL HIGH-RESOLUTION DIRECTION-FINDING TECHNIQUES

A fast separable approach based on a cross array is presented, which has coarsegrained parallelism. Its computational load is far less than that of the two-dimensional (2-D) direct processing method and other existing separable approaches. In order to compensate for the performance degradation due to separable processing, two postprocessing schemes are also proposed. Some computer simulation results are provided for illustration in the end.