Pleistocene-Holocene sediments interpreted by seismic refraction and wash-bore sampling, Plum Island - Castle Neck, Massachusetts /

"July 1973."

地震波谱元法数值模拟和偏移研究

With the development of both seismic theory and computer technology, numerical modeling technology of seismic wave has achieved great advancement during the past half century. The current methods under development include finite differentiation method (FDM), finite element method (FEM), pseudospectral method (PSM), integral equation method (IEM) and spectral element method (SEM). They exert their very important roles in every corner of seismology and seismic prospecting. Large quantity of researches towards spectral element method in the end of last century bring this method to a new era, which results in perfect solution of many difficult problems. However, parts of posterior works such as seismic migration and inversion which base on spectral element method have never been studied widely at least up to the present whereas are of importance to seismic imaging and seismic wave propagation. Based on previous work, this paper uses spectral element method to investigate the characteristics and laws of the seismic wave propagation in isotropic and anisotropic media. By thoroughly studying this high-accuracy method, we implement a kind of reverse-time pre- and post-stack migration based on SEM. In order to verify the validity of the SEM method, we have simulated the propagation of seismic wave in several different models. The simulation results show that: (1) spectral element method can be used to model any complex models and the computational results are comparable with the expected results and the analytic results; (2) the optimum accuracy can be achieved when the rank is between 4 and 9. When it is below 4, the dispersion may occur; and when it is above 9, the time step-length will be changed accordingly with the reducing space step-length in order to keep the computation stability. This will exponentially increase the computation time and at the same time the memory even if simulating the same media. This paper also applies explosive reflection surface imaging technology, time constancy principle of wave-filed extrapolation and least travetime raytracing technology of surface source to SEM pre- and post-stack migration of isotropic and anisotropic media. All imaging results derived by the above methods agree well with the real geological models and the position of interface and inflexions can also return to their right location well. This indicates that the method proposed in this paper is a kind of technology with high accuracy and robust stability. It can serve as an alternative method in real seismic data processing. All these work can boost the development of high-accuracy seismic imaging, and therefore have significant inference value.

基于广义Lippmann-Schwinger方程的地震波场模拟算法与应用

In the last several decades, due to the fast development of computer, numerical simulation has been an indispensable tool in scientific research. Numerical simulation methods which based on partial difference operators such as Finite Difference Method (FDM) and Finite Element Method (FEM) have been widely used. However, in the realm of seismology and seismic prospecting, one usually meets with geological models which have piece-wise heterogeneous structures as well as volume heterogeneities between layers, the continuity of displacement and stress across the irregular layers and seismic wave scattering induced by the perturbation of the volume usually bring in error when using conventional methods based on difference operators. The method discussed in this paper is based on elastic theory and integral theory. Seismic wave equation in the frequency domain is transformed into a generalized Lippmann-Schwinger equation, in which the seismic wavefield contributed by the background is expressed by the boundary integral equation and the scattering by the volume heterogeneities is considered. Boundary element-volume integral method based on this equation has advantages of Boundary Element Method (BEM), such as reducing one dimension of the model, explicit use the displacement and stress continuity across irregular interfaces, high precision, satisfying the boundary at infinite, etc. Also, this method could accurately simulate the seismic scattering by the volume heterogeneities. In this paper, the concrete Lippmann-Schwinger equation is specifically given according to the real geological models. Also, the complete coefficients of the non-smooth point for the integral equation are introduced. Because Boundary Element-Volume integral equation method uses fundamental solutions which are singular when the source point and the field are very close，both in the two dimensional and the three dimensional case, the treatment of the singular kernel affects the precision of this method. The method based on integral transform and integration by parts could treat the points on the boundary and inside the domain. It could transform the singular integral into an analytical one both in two dimensional and in three dimensional cases and thus it could eliminate the singularity. In order to analyze the elastic seismic wave scattering due to regional irregular topographies, the analytical solution for problems of this type is discussed and the analytical solution of P waves by multiple canyons is given. For the boundary reflection, the method used here is infinite boundary element absorbing boundary developed by a pervious researcher. The comparison between the analytical solutions and concrete numerical examples validate the efficiency of this method. We thoroughly discussed the sampling frequency in elastic wave simulation and find that, for a general case, three elements per wavelength is sufficient, however, when the problem is too complex, more elements per wavelength are necessary. Also, the seismic response in the frequency domain of the canyons with different types of random heterogeneities is illustrated. We analyzed the model of the random media, the horizontal and vertical correlation length, the standard deviation, and the dimensionless frequency how to affect the seismic wave amplification on the ground, and thus provide a basis for the choice of the parameter of random media during numerical simulation.

隐伏金矿床预测中的地震勘探险方法

Based on brief introduction of seismic exploration and it's general developing situation, the seismic exploration method in field work implementation and some problems frequently encountered in field, which should be pay attention to, are analyzed in detail. The most economic field work techniques are emphasized. Then the seismic data processing flow and it's interpretation technique about the processing results are presented. At last four examples of seismic prospecting in gold deposits are showed. The main conclusions of our research are: 1. Seismic prospecting technique is a very efficient method in the prediction of concealed gold deposits. With appropriate application, it can absolutely reflect the detail underground geological structure in the condition of rugged area and complicated geological environment. 2. The field geometry should be designed and changed according to different kinds of objective exploration depth and ground situation. The best field implementing parameters which include offset, the distance between two adjacent traces, the quantity of dynamite and the depth of hole for explosion, should be determined with examination. Only this way, the high quality original seismic data can be gotten. 3. In seismic data processing, the edition of invalid trace and source gather, signal enhancement, velocity analysis and migration are the key steps. It has some different points with conventional processing and needs a new processing flow and methods which is suitable to the data acquired in rugged area and complicated geological environment. 4. The new common reflection area stacking method in crooked line data processing is an efficient method to improve the signal to noise ratio of seismic data The innovations of our research work are: 1. In the areas which were considered to be forbidden zone, we implement the seismic exploration in several gold deposits in China through our application. All acquire distinguished effects. This show the seismic exploration method is a new effective method in the prediction of concealed gold deposits. 2. We developed a set of seismic field work techniques and data processing which is suitable to complex environment, especially find a effective method in stacking and noise elimination in crooked line data processing. 3. In the field of seismic profile interpretation, through our research work, we are convinced of that: in different kinds of geological condition, the seismic reflection character are not same. For example the lava, the intrusion rock and sediment layers are different in the character of reflection structure and strength. So we accumulate some experience about seismic data interpretation in the area of gold deposits.

济阳地区深部结构与含油气性讨论

Oil and gas exploration of marine strata in China's Pre-Cenozoic residual basins is regarded as a worldwide puzzle because of existent problems and cruxes. Objectively speaking, the subsurface geologic structure is complicated, and the surface conditions of some areas are tough. On the other hand, there are still many problems to be solved in oil and gas exploration technologies of Pre-Cenozoic marine fades, and theoretic cognition about petroleum geology is not profound yet. Therefore, it is principal to explore integrated geophysical research ways of Pre-Cenozoic residual basins. Seismic prospecting and geophysical integrated interpretation technologies aimed at middle Paleozoic marine facies with deeper burial and complicated geologic conditions have not formed due to bad quality of deep strata data. Pre-Cenozoic strata, and especially extension, thickness and internal structure of Paleozoic strata can not be recognized from seismic profiles, so it is hard to systematically cognize structural features and oil-gas resources prospect of Pre-Cenozoic basins. To further investigation of fabric and structural features, basin prototype, formation and evolution pattern of Pre-Cenozoic basins, and also their control over formation, migration and aggregation of oil and gas, will play a guiding and promotive role in developing new surveying areas, selecting advantageous zones and predicting oil-gas resources.This paper follows the modem macrocontrol theory of "Region constrains local, deep strata controls shallow ones", and uses the integrated geophysical method of "One guide, two hinges, three combinations and multi feedbacks'*. Based on several years of geological and geophysical results of the Shengli Oilfield, and 14 newly-joint regional seismic profiles, deep structure and oil-gas bearing capacity of the Jiyang area are discussed and new cognitions are drawn as below.Seismic identification marks Tr, Tg, Tgl and Tg2 are established for importantPre-Cenozoic geological interfaces, and promoted to the whole Jiyang area.Through area-wide tracking and clogging of important seismic reflection marker,the isochronic framework of pre-Tertiary basin is set up in the Jiyang area for the firsttime, which is vital for basin research.Integrated with geological and geophysical research results, the Jiyang area isdivided into four first-order tectonic sequences- basement, lower tectonic layer,upper tectonic layer, and top tectonic layer. The basement and lower tectonic sequence which are related to Pre-Cenozoic are studied with emphasis.Through the research of regional seismic profiles, the point of view is given thatthe Kongdian Formation of Jiyang is structural transition period. The top-bottomunconformable interface of the Kongdian Formation is found out for the first time,and the basin model is determined primarily, which lay a basis for prototype basinresearch of the Jiyang Kongdian Formation.The distribution status of Middle-Paleozoic is delineated in the Jiyang area.The maximum thickness of Paleozoic lies in the top of the south declivity of half-graben. The thickness gets thinner towards the center of Mesozoic and Cenozoic half-graben basin, and even disappears. Structural action in the west-north affects the distribution of Paleozoic residual strata.6. The features of second-order tectonic sequence of the Jiyang depression isstudied and its evolution history of is rebuilt.Combined with the 5-stage evolution history of the China continent and structure evolution features of the Jiyang area, the structure sedimentary process since Paleozoic is divided into 5 periods - basement forming , Indosinian orogenic, Yanshan negative reversal, Himalayan extension and Neogene subsidence period.Combined with the research results of gravity, magnetic surveying and regionalprofiles, this paper brings forward the idea for the first time that the western boundaryof the Jiyang depression is the Ningjin-Yangpan fracture zone, and forms aside-column assemblage with the Wudi fracture zone.The opinion that under Middle-Cenozoic basins in the middle Jiyang area theremight superimpose an old residual basin is given for the first time. And if it is provedto be true, a new exploration space will be pioneered for Jiyang and even north China.There exists many types of tectonic-stratigraphic traps formed under piezotropy,extension and compound action in Pre-Cenozoic Jiyang. Therein all kinds of burialhills are the most important oil-gas trap type of Pre-Cenozoic, which should besurveyed layeredly according to the layout of oil sources.As such a new challenging project and field, the paper systematically analyses different geophysical responses of the Jiyang area, frames the deep structure of the area, and preliminarily recognizes the Pre-Cenozoic residual basins. It breaks through to a certain extent in both theory and practice, and is expected to provide new geophysical and geotectonic clues for deep exploration in Shengli.

金属矿集区深层结构高分辨率地震探测——以安徽铜陵金属矿为例

The reflection seismic prospecting technique is an important and a widely used method in the petroleum and coal surveying, and has been developed to a perfectly mature technique from the aspects of data acquisition, data processing to data interpretation. However, the metallic mine seismic prospecting, especially the high resolution seismic prospecting technique are being still in the course of studying and probing up to now. In this paper, the basic theory and the present situation of study on metallic mine seismic reflection are expatiated, the basic theory, the improving measure, the converging velocity and the ability on the integrating global optimization method are also illuminated in detail at first. Then the basic theory, the realization process and the practicing effects of the vector suppressing noise algorithm are also introduced. On the basis of studying of applying the integrating global optimization method to static correction and the vector suppressing noise algorithm, we elaborate processed the seismic data of Tongling metallic mine. We introduced the processing flow, the key steps and the processing effects. Basing on the processing results, we analyzed the major reflection characteristics, the geological interpretation results and the earth's crust top reflection structure and the space distribution status of Wutong set, the space shape of part lithological body and the contacting relations of horizonsunveiled.

宽频带地震数据采集系统的研制

The seismic data acquisition system is the most important equipment for seismic prospecting. The geophysicists have been paying high attention to the specification of the equipment used in seismic prospecting. Its specification and performance are of great concerned to acquire precisely and accurately seismic data, which show us stratum frame. But, by this time, limited by the technology, most of the Broad-band Seismic Recorder (BSR) for lithosphere research of our country were bought from fremdness which were very costliness and maintained discommodiously. So it is very important to study the seismic data acquisition system.The subject of the thesis is the research of the BSR, several items were included, such as: seismic data digitizer and its condition monitor design.In the first chapter, the author explained the significance of the implement of BSR, expatiated the requirement to the device and introduced the actuality of the BSR in our country.In the second chapter, the collectivity architecture of the BSR system was illustrated. Whereafter, the collectivity target and guideline of the performance of the system design were introduced. The difficulty of the system design and some key technology were analyzed, such as the Electro Magnetic Compatibility (EMC), system reliability technology and so on.In the third chapter, some design details of BSR were introduced. In the recorder, the former analog to digital converter (ADC) was separated from the later data transition module. According to the characteristic of seismic data acquisition system, a set high-resolution 24-bit ADC chip was chosen to the recorder design scheme. As the following part, the noise performance of the seismic data channel was analyzed.In the fourth chapter, the embedded software design of each board and the software design of the workstation were introduced. At the same time the communication protocol of the each module was recommendedAt the last part of this thesis, the advantages and the practicability of the BSR system design were summarized, and the next development items were suggested.

Amplitudes e padrões de polarização de pulsos em meios anisotrópicos

Extrair informações litológicas da subsuperfície através de dados sísmicos constitui-se num grande desafio à prospecção sísmica, pois a hipótese de estratificações formadas por camadas isotrópicas se mostra insuficiente para representar o comportamento do campo elástico em levantamentos com grandes afastamentos entre fonte e receptor, geofones multicomponentes, medidas de VSP tridimensional, entre outros. Sob este panorama, a prospecção sísmica passa a considerar modelos anisotrópicos de subsuperfície para, por exemplo, caracterizar reservatórios. O objetivo deste texto é apresentar um formalismo para modelar o espalhamento de pulsos a partir de ondas planas incidentes em interfaces planas horizontais que separam meios anisotrópicos. Este espalhamento é obtido primeiramente, através da formulação explícita dos campos de deformação e tração como função das matrizes propagadoras, de polarização e de impedância do meio. Em seguidaeste formalismo é usado para a obtenção das matrizes dos coeficientes de reflexão e transmissão através de uma interface plana horizontal para posteriormente, ser generalizado para o espalhamento através de múltiplas camadas. Finalmente, inserem-se ao campo da onda incidente as amplitudes de um pulso analítico para calcular o espalhamento do pulso através de estratificações.

A comparison of cross-hole electrical and seismic data in fractured rock

Cross-hole anisotropic electrical and seismic tomograms of fractured metamorphic rock have been obtained at a test site where extensive hydrological data were available. A strong correlation between electrical resistivity anisotropy and seismic compressional-wave velocity anisotropy has been observed. Analysis of core samples from the site reveal that the shale-rich rocks have fabric-related average velocity anisotropy of between 10% and 30%. The cross-hole seismic data are consistent with these values, indicating that observed anisotropy might be principally due to the inherent rock fabric rather than to the aligned sets of open fractures. One region with velocity anisotropy greater than 30% has been modelled as aligned open fractures within an anisotropic rock matrix and this model is consistent with available fracture density and hydraulic transmissivity data from the boreholes and the cross-hole resistivity tomography data. However, in general the study highlights the uncertainties that can arise, due to the relative influence of rock fabric and fluid-filled fractures, when using geophysical techniques for hydrological investigations.