地震采集观测系统的共聚焦分析

Focal beam analysis is a method for assessment of acquisition geometries that is directly linked to pre-stack migration. About dealing with the complex subsurface structures, the conventional survey design methods which do not take into account the subsurface are no longer valid. Based on the Fourier finite-difference (FFD) large-step wave field extrapolation and Born-Kirchhoff (BK) small-step wavefield interpolation, the thesis presents a rapid resolution analysis of 3D seismic survey design by focal beams in complicated media. Subsequently, The SEG/EAEG salt model is used to illustrate the method. Based on the focal beam resolution definition, each kind of influence factor is discussed. The focal beam analysis usually is carried out in a single frequency, but the actual seismic waves always contain a frequency bandwidth. In this thesis, theoretical relationship between focal beam analysis and frequency is derived. Since the effects of focal beam analysis are linear with frequency simply, the multi-frequency focal beam analysis using interpolation is developed. At the same time, the resolution of different frequency bandwidth is interconvertible in accordance with Signal uncertainty principle. The resolution of all frequency bands can be calculated by using only a few focal beam analysis for a seismic survey. In the last section of this thesis, I propose a new approach to predicting acquisition footprint, based on the assumption of Common-Middle-Point stack without constructing a special velocity model. The approach is a simplistic analytical method in which the acquisition footprint pattern is a weighted, linear summation of limited-offset fold-of-stack plots. Because the value of acquisition can be got by quantificational and rapidly calculating, we can exactly do a comparative analysis among different plans of seismic survey by this method.

反射地震转换波成像方法研究

Multi-waves and multi-component get more and more attentions from oil industry. On the basis of existent research results, My research focuses on some key steps of OBC 4C datum processing. OBC datum must be preprocessed quite well for getting a good image. We show a flow chart of preprocess including attenuation of noise on multi-component datum、elimination ghost by summing P and Z and rotation of horizontal components. This is a good foundation for the coming steps about OBC processing. How to get exact converted point location and to analyze velocity are key points in processing reflection seismic converted wave data. This paper includes computing converted point location, analyzing velocity and nonhyperbolic moveout about converted waves. Anisotropic affects deeply the location of converted wave and the nonhyperbolic moveout. Supposed VTI, we research anisotropic effect on converted wave location and the moveout. Since Vp/Vs is important, we research the compute method of Vp/Vs from post-stack data and pre-stack data. It is a part of the paper that inversing anisotropic parameter by traveltime. Pre-stack time migration of converted wave is an focus, using common-offset Kirchhoff migration, we research the velocity model updating in anisotropic media. I have achieved the following results: 1) using continued Fractions, we proposed a new converted point approximate equation, when the offset is long enough ,the thomsen’s 2 order equation can’t approximate to the exact location of converted point, our equation is a good approximate for the exact location. 2) our new methods about scanning nonhyperbolic velocity and Vp/Vs can get a high quality energy spectrum. And the new moveout can fit the middle and long offset events. Processing the field data get a good result. 3) a new moveout equation, which have the same form as Alkhalifah’s long offset P wave moveout equation, have the same degree preciseness as thomsen’s moveout equation by testing model data. 4) using c as a function of the ratio offset to depth, we can uniform the Li’s and thomsen’s moveout equation in a same equation, the model test tell us choice the reasonable function C can improve the exact degree of Li’s and thomsen’s equation. 5) using traveltime inversion ,we can get anisotropic parameter, which can help to flat the large offset event and propose a model of anisotropic parameter which will useful for converted wave pre-stack time migration in anisotropic media. 6）using our pre-stack time migration method and flow, we can update the velocity model and anisotropic parameter model then get good image. Key words: OBC, Common converted Point (CCP), Nonhyperbolic moveout equation, Normal moveout correction, Velocity analysis, Anisotropic parameters inversion, Kirchhoff anisotropic pre-stack time migration, migration velocity model updating

叠前地震多角度三参数同步反演方法

Pre-stack seismic inversion has become the emphasis and hotspot owing to the exploration & exploitation of oil field and the development of seismic technology. Pre-stack seismic inversion has the strongpoint of making the most of amplitude versus offset compared with the post-stack method. In this dissertation, the three parameters were discussed from multi-angle reflectance of P-wave data based on Zoeppritz’s and Aki & Richard’s equation, include P-wave velocity, S-wave velocity, and density. The three parameters are inversed synchronously from the pre-stack multi-angle P-wave data, based on rockphysics model and aimed at the least remnant difference between model simulation and practical data. In order to improve the stability of inversion and resolution to thin bed, several techniques were employed, such as the wavelet transform with multi-scale function, adding the Bayesian soft constraint and hard constraints (the horizon, structure and so on) to the inversion process. Being the result, the uncertainty of the resolution is reduced, the reliability and precision are improved, the significance of parameters becomes clearer. Meeting to the fundamental requirement of pre-stack inversion, some research in rockphysics are carried out which covered the simulation and inversion of S-wave velocity, the influence of pore fluids to geophysical parameters, and the slecting and analyzing of sensitive parameters. The difference between elastic wave equation modeling and Zoeppritz equation method is also compared. A series of key techniques of pre-stack seismic inversion and description were developed, such as attributes optimization, fluid factors, etc. All the techniques mentioned above are assembled to form a technique sets and process of synchronous pre-stack seismic inversion method of the three parameters based on rock physics and model simulation. The new method and technology were applied in many areas with various reservoirs, obtained both geological and economic significance, which proved to be valid and rational. This study will promote the pre-stack inversion technology and it’s application in hidden reservoirs exploration, face good prospects for development and application.

叠前时间偏移：各向异性、多次波衰减及应用

At present, in order to image complex structures more accurately, the seismic migration methods has been developed from isotropic media to the anisotropic media. This dissertation develops a prestack time migration algorithm and application aspects for complex structures systematically. In transversely isotropic media with a vertical symmetry axis (VTI media), the dissertation starts from the theory that the prestack time migration is an approximation of the prestack depth migration, based on the one way wave equation and VTI time migration dispersion relation, by combining the stationary-phase theory gives a wave equation based VTI prestack time migration algorithm. Based on this algorithm, we can analytically obtain the travel time and amplitude expression in VTI media, as while conclude how the anisotropic parameter influence the time migration, and by analyzing the normal moveout of the far offset seismic data and lateral inhomogeneity of velocity, we can update the velocity model and estimate the anisotropic parameter model through the time migration. When anisotropic parameter is zero, this algorithm degenerates to the isotropic time migration algorithm naturally, so we can propose an isotopic processing procedure for imaging. This procedure may keep the main character of time migration such as high computational efficiency and velocity estimation through the migration, and, additionally, partially compensate the geometric divergence by adopting the deconvolution imaging condition of wave equation migration. Application of this algorithm to the complicated synthetic dataset and field data demonstrates the effectiveness of the approach. In the dissertation we also present an approach for estimating the velocity model and anisotropic parameter model. After analyzing the velocity and anisotropic parameter impaction on the time migration, and based on the normal moveout of the far offset seismic data and lateral inhomogeneity of velocity, through migration we can update the velocity model and estimate the anisotropic parameter model by combining the advantages of velocity analysis in isotropic media and anisotropic parameter estimation in VTI media. Testing on the synthetic and field data, demonstrates the method is effective and very steady. Massive synthetic dataset、2D sea dataset and 3D field datasets are used for VTI prestack time migration and compared to the stacked section after NMO and prestack isotropic time migration stacked section to demonstrate that VTI prestack time migration method in this paper can obtain better focusing and less positioning errors of complicated dip reflectors. When subsurface is more complex, primaries and multiples could not be separated in the Radon domain because they can no longer be described with simple functions (parabolic). We propose an attenuating multiple method in the image domain to resolve this problem. For a given velocity model，since time migration takes the complex structures wavefield propagation in to account, primaries and multiples have different offset-domain moveout discrepancies, then can be separated using techniques similar to the prior migration with Radon transform. Since every individual offset-domain common-reflection point gather incorporates complex 3D propagation effects, our method has the advantage of working with 3D data and complicated geology. Testing on synthetic and real data, we demonstrate the power of the method in discriminating between primaries and multiples after prestack time migration, and multiples can be attenuated in the image space considerably.

人工源宽角地震资料中的广角效应及其对地壳速度成像的影响

The seismic wide-angle reflection/refraction method is the one of the most effective method for probing the crustal and upper mantle structure. It mainly uses the wide-angle reflection information from the boundary in the crust and the top boundary of the upper mantle to rebuild the crust and upper mantle structure. Through analyzing the reflection and transmission coefficients of various incident waves on the interface, we think relative to the pre-critical angle reflection information the post critical angle reflection information that received by wide-angle seismic data exists a time-shift effect with the offset variation, and then it must cause the error for velocity analysis and structure image. The feature of the wide-angle seismic wave field of the fourteen representative crust columns tell us that the wide-angle effects in the different representative tectonic units for the interface depth and the interval velocity in crust. We studied the features of the wide-angle seismic wave field through building the crust model and inverse its travel time by GA method to know the wide-angle influence on crustal velocity image. At last we finished the data processing of the Tunxi-Wenzhou wide-angle seismic profile. The results are as following: (1) Through building crust model, we labeled the travel time for all the phases by ray tracing method and remove wide-angle effects method, it revealed the wide-angle effect exists in the seismic data. (2) The travel time inversion by GA method can tell us that the depth by traditional ray tracing method is shallower than the result by remove wide-angle effects method, the latter can recover the crust structure model in effect. (3) We applied the two method mentioned before to the fourteen representative crust columns in China. It indicates that the removed wide-angle effect method in travel time inversion is reasonable and effective. (4) The real data processing from Tunxi-Wenzhou wide-angle seismic profile give us the basic structure through the two ways. The main influence exhibits in the difference of the interval velocity of the curst, and the wide-angle effects in shallow interface are stronger than the deep interface.

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

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.

三维非均匀各向异性介质中的纵波走时反演研究

Many observations show that seismic anisotropy is very common in the crust and upper mantle of the Earth. Seismic anisotropy can provide some clue about the changing and transporting process inside the earth. in recent years, abundant earthquake travel time data are accumulated, computers become more powerful, and these make the inversion of earthquake travel time data practical. In this thesis we studied the theory of elastic wave in anisotropic media, some formule for travel time inversion were derived. We present an iterative procedure to determine 21 elastic parameters from qP wave travel times. No a priori assumptions about heterogeneity and anisotropy of the model are made. The procedure is suitable for the case when we know nothing about the symmetry of anisotropy of the media, as well as for the case of earthquake travel time inversion which may contain various symmetry of anisotropy. The procedure is tested with a synthetic multiple-source offset VSP experiment. The results proved that the formulae are correct, and the procedure is practical. The results and the related theory indicate that the anisotropic inversion needs more rays than isotropic case. For a 2-D weak anisotropic (WA) medium, we need at least 5 rays in different directions to retrieve the elastic parameters on one grid point, and for a 3-D WA medium we need at least 15 rays in different directions to retrieve the elastic parameters on one grid point. The results also indicate that the starting background velocity has no influence on the final results, at least for the model we specified. Our results also show that insufficient illumination coverage will slow down the convergence rate, and make the results more sensitive to noise. We apply the procedure to a set of field travel time data. The data is from an artificial seismic observation. This observation is for locating micro-seismic events around a tunnel, its purpose is to find out if the digging process and the stress condition around the tunnel can generate micro-cracks. The size of this area is around 100m. The anisotropy derived from qP travel times is the same as the anisotropy showed by apparent velocities, and is also consistent with the anisotropy derived from S-wave splitting phenomena.

大别－苏鲁造山带与华北陆块东部岩石层三维速度结构

The 3-D velocity images of the crest and upper mantle beneath the region of 112° -124°E, 28°-39°N including the Dabie-Sulu orogenic belt are reconstructed by using 36405 P-wave arrivals of 3437 regional and 670 distant earthquakes during the period from 1981 to 1996, and gridding the area of 0.5° * 0.5°. The results of tomography demonstrate that: 1. The results of tomographic imaging show a broad heterogeneity in P wave velocity structure for the lithosphere beneath the Dabie-Sulu orogenic belt. 2. In the Dabie orogenic belt, the velocity patterns in the crust are different among various tectonic units. The Dabie and Qinling orogenic belts are remarkable in the tomographic images, and in mm the Hongan and Dabie blocks in the Dabie orogenic belt are also imaged very distinguishably. 3. A velocity (about 5.9～6.0 km/s) layer exists in the Dabie block at depth between 15～25 km, which is coincident with the low-resistance layer at the depth of 12-23 km, being inferred to be the tectonic detachment zone and suggesting that the extension detachment structure was formed in the middle crust. Beneath the southern and northerm Dabie tectonic units, the north-dipping high-velocity (at level of 6.5 ～ 6.6 km/s) block was developed in the crust, which might be correlated with the UHP rockswith low content of the meta-ultramafic rocks. This result is in agreement with the geological observation on the surface. 4. The velocity image at 40 km depth reveals the features at the top of mantle and the configuration of the Moho discontinuity. The depth of the Moho changes slightly along the trend of the orogenic belt. It in Hongan block is less than 40 km, but it is different in the western and eastern parts of the Dabie block, the former is more than 40 km, and the latter less than or equal to 40 km. The remnant of the mountain root exists between the Shangcheng-Macheng fault and the line of Huoshan-Yuexi-Yingshan in the Dabie orogenic belt, and beneath the southern and northern Dabie tectonic units. However, the thickness of the Moho is about 40 km and there is no obvious changes, which suggest that the Dabie orogenic belt has been experienced quite in the gravity equilibration. The Moho's depth in the Sulu is less than 40 km. 5. There is a dipping slab-like high-velocity body in the uppermost mantle. It is sandwiched by slow velocities and exists beneath the Dabie-Sulu orogenic belt in the range of depths between the Moho discontinuity and 110 km at least. This high-velocity body outlines a picture of the slab interpreted as the remnant of the Triassic subducted YZ. 6. The Sulu orogenic belt displays "crocodilian" velocity structure, the upper crust of the Yangtze thrusted over the Huabei crest, and the Huabei crust indented into the Yangtze crust, where the ancient subduction zone of the Yangtze lithosphere located. Based on the previous geological data, this structure is not related with the collision between the Yangtze and Sino-Korean Blocks, but caused by the sinistral offset of the Tan-Lu Fault. Studied on the velocity structure of the eastern Huabei lithosphere indicates: 1. The 'present-day' lithosphere of the eastern Huabei is between 40-100 km thick with greatly thinned lithosphere around the Bohai Sea. Generally, thickness of the lithosphere in this region decreased eastwards. 2. The attenuation of the lithosphere is attributed to the strongly uplift of the asthenosphere. In the area between the Taihang Mountains and the Tan-Lu Fault, there is a 'lever' with red low velocity belt, it is clearly defined, transverse continuity, depth between 100-150 km, local variations visible, and an upwards trend towards the Bohai Sea. Generally, the velocity structure in the mantle beneath the lithosphere displays irregular column-shape consisting of alternating high and low velocities, and when cold high velocity ancient lithosphere connects with the hot low velocity mantle materials forming precipitous compact structure. More heat pathways from the mantle occur towards the Tan-Lu Fault. 3. The strongly irregular characteristics of the contact between the asthenosphere and the lithosphere is induced by the long-term hot, chemical erosion and alteration on the contact. 4. There are still preserved high velocity lithosphedc root beneath Huabei with 'block-shape' distribution and surrounded by hot materials. Results of our studies indicate that the evolution models of the eastern China mantle are characterized by the direct contact between the uplifted lithosphere and the Huabei Craton accompanying the upwelling of the deep mantle materials. At the contact betwen the lithosphere and the asthenosphere, the upwelled mantle materials replaced and altered the lower lithosphere forming the metasome through the hot and chemical modifications impacted on the Craton lithosphere, and changed it into the lithosphere gradually, resulting in the lithospheric thinning. Thus, the lithospheric thinning is the result of the upwelling of the asthenosphere.

末次盛冰期以来中国陆地生态系统碳储量变化研究

Carbon is an essential element for life, food and energy. It is also a key element in the greenhouse gases and therefore plays a vital role in climatic changes. The rapid increase in atmospheric concentration of CO_2 over the past 150 years, reaching current concentrations of about 370 ppmv, corresponds with combustion of fossii fuels since the beginning of the industrial age. Conversion of forested land to agricultural use has also redistributed carbon from plants and soils to the atmosphere. These human activities have significantly altered the global carbon cycle. Understanding the consequences of these activities in the coming decades is critical for formulating economic, energy, technology, trade, and security policies that will affect civilization for generations. Under the auspices of the International Geosphere-Biosphere Programme (IGBP), several large international scientific efforts are focused on elucidating the various aspects of the global carbon cycle of the past decade. It is only possible to balance the global carbon cycle for the 1990s if there is net carbon uptake by terrestrial ecosystems of around 2 Pg C/a. There are now some independent, direct evidences for the existence of such a sink. Policymarkers involved in the UN Framework Convention on Climate Change (UN-FCCC) are striving to reach consensuses on a 'safe path' for future emissions, the credible predictions on where and how long the terrestrial sink will either persist at its current level, or grow/decline in the future, are important to advice the policy process. The changes of terrestrial carbon storage depend not only on human activities, but also on biogeochemical and climatological processes and their interaction with the carbon cycles. In this thesis, the climate-induced changes and human-induced changes of carbon storage in China since the past 20,000 years are examined. Based on the data of the soil profiles investigated during China's Second National Soil Survey (1979-1989), the forest biomass measured during China's Fourth National Forest Resource Inventory (1989-1993), the grass biomass investigated during the First National Grassland Resource Survey (1980-1991), and the data collected from a collection of published literatures, the current terrestrial carbon storage in China is estimated to -144.1 Pg C, including -136.8 Pg C in soil and -7.3 Pg C in vegetation. The soil organic (SOC) and inorganic carbon (SIC) storage are -78.2 Pg C and -58.6 Pg C, respectively. In the vegetation reservoir, the forest carbon storage is -5.3 Pg C, and the other of-1.4 Pg C is in the grassland. Under the natural conditions, the SOC, SIC, forest and grassland carbon storage are -85.3 Pg C, -62.6 Pg C, -24.5 Pg C and -5.3 Pg C, respectively. Thus, -29.6 Pg C organic carbon has been lost due to land use with a decrease of -20.6%. At the same time, the SIC storage also has been decreased by -4.0 Pg C (-6.4%). These suggest that human activity has caused significant carbon loss in terrestrial carbon storage of China, especially in the forest ecosystem (-76% loss). Using the Paleocarbon Model (PCM) developed by Wu et al. in this paper, total terrestrial organic carbon storage in China in the Last Glacial Maximum (LGM) was -114.8 Pg C, including -23.1 Pg C in vegetation and -86.7 Pg C in soil. At the Middle Holocene (MH), the vegetation, soil and total carbon were -37.3 Pg C, -93.9 Pg C and -136.0 Pg C, respectively. This implies a gain of-21.2 Pg C in the terrestrial carbon storage from LGM to HM mainly due to the temperature increase. However, a loss of-14.4 Pg C of terrestrial organic carbon occurred in China under the current condition (before 1850) compared with the MH time, mainly due to the precipitation decrease associated with the weakening of the Asian summer monsoon. These results also suggest that the terrestrial ecosystem in China has a substantial potential in the restoration of carbon storage. This might be expected to provide an efficient way to mitigate the greenhouse warming through land management practices. Assuming that half of the carbon loss in the degraded terrestrial ecosystem in current forest and grass areas are restored during the next 50 years or so, the terrestrial ecosystem in China may sequestrate -12.0 Pg of organic carbon from the atmosphere, which represents a considerable offset to the industry's CO2 emission. If the ' Anthropocene' Era will be another climate optimum like MH due to the greenhouse effect, the sequestration would be increased again by -4.3 - 9.0 Pg C in China.

地震波场延拓+菲涅尔体叠加方法及成像效果

In the increasingly enlarged exploration target, deep target layer(especially for the reservoir of lava) is a potential exploration area. As well known, the reflective energy becomes weak because the seismic signals of reflection in deep layer are absorbed and attenuate by upper layer. Caustics and multi-values traveltime in wavefield are aroused by the complexity of stratum. The ratio of signal to noise is not high and the fold numbers are finite(no more than 30). All the factors above affect the validity of conventional processing methods. So the high S/N section of stack can't always be got with the conventional stack methods even if the prestack depth migration is used. So it is inevitable to develop another kind of stack method instead. In the last a few years, the differential solution of wave equation was hold up by the condition of computation. Kirchhoff integral method rose in the initial stages of the ninetieth decade of last century. But there exist severe problems in it, which is are too difficult to resolve, so new method of stack is required for the oil and gas exploration. It is natural to think about upgrading the traditionally physic base of seismic exploration methods and improving those widely used techniques of stack. On the other hand, great progress is depended on the improvement in the wave differential equation prestack depth migration. The algorithm of wavefield continuation in it is utilized. In combination with the wavefield extrapolation and the Fresnel zone stack, new stack method is carried out It is well known that the seismic wavefield observed on surface comes from Fresnel zone physically, and doesn't comes from the same reflection points only. As to the more complex reflection in deep layer, it is difficult to describe the relationship between the reflective interface and the travel time. Extrapolation is used to eliminate caustic and simplify the expression of travel time. So the image quality is enhanced by Fresnel zone stack in target. Based on wave equation, high-frequency ray solution and its character are given to clarify theoretical foundation of the method. The hyperbolic and parabolic travel time of the reflection in layer media are presented in expression of matrix with paraxial ray theory. Because the reflective wave field mainly comes from the Fresnel Zone, thereby the conception of Fresnel Zone is explained. The matrix expression of Fresnel zone and projected Fresnel zone are given in sequence. With geometrical optics, the relationship between object point in model and image point in image space is built for the complex subsurface. The travel time formula of reflective point in the nonuniform media is deduced. Also the formula of reflective segment of zero-offset and nonzero offset section is provided. For convenient application, the interface model of subsurface and curve surface derived from conventional stacks DMO stack and prestack depth migration are analyzed, and the problem of these methods was pointed out in aspects of using data. Arc was put forward to describe the subsurface, thereby the amount of data to stack enlarged in Fresnel Zone. Based on the formula of hyperbolic travel time, the steps of implementation and the flow of Fresnel Zone stack were provided. The computation of three model data shows that the method of Fresnel Zone stack can enhance the signal energy and the ratio of signal to noise effectively. Practical data in Xui Jia Wei Zhi, a area in Daqing oilfield, was processed with this method. The processing results showed that the ability in increasing S/N ratio and enhancing the continuity of weak events as well as confirming the deep configuration of volcanic reservoir is better than others. In deeper target layer, there exists caustic caused by the complex media overburden and the great variation of velocity. Travel time of reflection can't be exactly described by the formula of travel time. Extrapolation is bring forward to resolve the questions above. With the combination of the phase operator and differential operator, extrapolating operator adaptable to the variation of lateral velocity is provided. With this method, seismic records were extrapolated from surface to any different deptlis below. Wave aberration and caustic caused by the inhomogenous layer overburden were eliminated and multi-value curve was transformed into the curve.of single value. The computation of Marmousi shows that it is feasible. Wave field continuation extends the Fresnel Zone stack's application.

共反射面元叠加

The seismic survey is the most effective prospecting geophysical method during exploration and development of oil/gas. The structure and the lithology of the geological body become increasingly complex now. So it must assure that the seismic section own upper resolution if we need accurately describe the targets. High signal/noise ratio is the precondition of high-resolution. As one important seismic data processing method, Stacking is an effective means to suppress the records noise. Broadening area of surface stacked is more important to enhance genuine reflection signals and suppressing unwanted energy in the form of coherent and random ambient noise. Common reflection surface stack is a macro-model independent seismic imaging method. Based on the similarity of CRP trace gathers in one coherent zone, CRS stack effectively improves S/N ratio by using more CMP trace gathers to stack. It is regarded as one important method of seismic data processing. Performing CRS stack depends on three attributes. However, the equation of CRS is invalid under condition of great offset. In this thesis, one method based on velocity model in depth domain is put forward. Ray tracing is used to determine the traveltime of CRP in one common reflection surface by the least squares method to regress the equation of CRS. Then we stack in the coherent seismic data set according to the traveltime, and get the zero offset section. In the end of flowchart of implementing CRS stack, one method using the dip angle to enhance the ratio of S/N is used. Application of the method on synthetic examples and field seismic records, the results of this method show an excellent performance of the algorithm both in accuracy and efficiency.

利用宽频带台阵资料研究华北地区地壳及上地幔间断面结构

Intense tectonic renovation has occurred in the eastern continent of china since Mesozoic, as evidenced by the high heat flow, widespread magma extrusion and volcanic activities, and development of large sedimentary basins. To explain the cause and mechanism for the tectonic process in this period, some researchers have put forward various models, such as mantle plume, subduction of the Pacific slab, Yangtze Block-North China Block collision, etc. Their seismological evidence, however, is still scarce..During the period from 2000 to 2003, large temporary seismic arrays were established in North China by the Institute of the Geology and Geophysics, Chinese Academy of Sciences. Total 129 portable seismic stations were linearly emplaced across the western and eastern boundaries of the Bohai Bay Basin, and accumulated a large amount of high-quality data. Moreover, abundant data were also collected at the capital digital seismic network established in the ninth five-year period of national economic and social development. These provide an unprecedented opportunity for us to study the deep structure and associated geodynamic mechanism of lithospheric processes in North China using seismological techniques.Seismology is a kind of observation-based science. The development of seismic observations greatly promotes the improvement of seismologic theory and methodology. At the beginning of this thesis, I review the history of seismic observation progress, and present some routine processing techniques used in the array seismology. I also introduce two popular seismic imaging methods (receiver function method and seismic tomography).Receiver function method has been widely used to study the crustal and upper mantle structures, and many relevant research results have been published. In this thesis I elaborate the theory of this method, including the basic concept of receiver functions and the methodology for data pre-processing, stacking and migration. I also address some problems often encountered in practical applications of receiver function imaging.By using the teleseismic data collected at the temporary seismic arrays in North China, in particular, the traveltime information of P-to-S conversion and multiple reverberations of the Moho discontinuity, I obtain the distributions of the crustal thickness and the poisson ratio at the northwest boundary area of the Bohai Bay Basin and discuss the geological implications of the results.Through detailed intestigations on the crustal structural feature around the middle part of the Tanlu fault, considerable disparity in poisson ratios is found in the western and eastern sides of the Tanlu fault. Moreover, an obvious Moho offset is coincidently observed at the same surface location. A reasonable density model for the Tanlu fault area is also derived by simulating the observed gravity variations. Both receiver function study and gravity anomaly modeling suggest that the crustal difference between the western and eastern sides of the Tanlu fault is mainly resulted from their different compositions.With common conversion point imaging of receiver functions, I estimate the depths of the upper and lower boundaries of the mantle transition zone, i.e., the 410 and 660 km discontinuities, beneath most part of the North China continent The thickness of the transition zone (TTZ) in the study area is calculated by subtracting the depth of .410 km discontinuity from that of the 660km discontinuity. The resultant TTZ is 10-15 km larger in the east than in the west of the study area. Phase transitions at the 410 km and the 660 km discontinuities are known to have different Clapeyron slopes. Therefore, the TTZ is sensitive to the temperature changes in the transition zone. Previous studies have shown that the TTZ would be smaller in the mantle plume areas and become larger when the remnants of subducted slabs are present The hypothesis of mantle plume cannot give a reasonable interpretation to the observed TTZ beneath North China, Instead, the receiver function imaging results favor a dynamic model that correlates the thermal structure of the mantle transition zone and associated upper mantle dynamics of North China to the Pacific plate subduction process.

注意在视觉刺激位置编码中的作用

There is extensive agreement that attention may play a role in spatial stimmlus coding (Lu & Proctor, 1995). Some authors investigated the effects of spatial attention on the spatial coding by using spatial cueing procedure and spatial Stroop task. The finding was that the stroop effects were modulated by spatial cueing. Three hypotheses including attentional shift account, referential coding account, and event integration account were used to explain the modulation of spatial cueing over the spatial Stroop effects. In these previous studies, on validly cued trials, cue and target not only appeared at the same location, but also in the same object, which resulted in both location and object cued. Consequently, the modulation of spatial attentional cueing over spatial Stroop effects was confounded with the role of object-based attention. In the third chapter of this dissertation, using a modification of double rectangles cueing procedure developed by Egly, Driver and Rafal (1994) and spatial Stroop task employed by Lupiáñez and Funes (2005), separate effects of spatial attention and object-based attention on the location code of visual stimuli were investigated. Across four experiments, the combined results showed that spatial Stroop effects were modulate by object-based attention, but not by location-based attention. This pattern of results could be well explained by event integration account, but not by attentional shift account and referential coding account. In the fourth chapter, on the basis of the prior chapter, whether the modulation of attentional cueing on location code occurred at the stage of perceptual identification or response choice was investigated. The findings were that object-based attention modulated spatial Stroop effects and did not modulate the Simon effects, whereas spatial attention did not modulate Stroop and Simon effects. This pattern of results partially replicated the outcome of the previous chapter. The previous studies generally argued that the conflicts of spatial Stroop task and Simon task respectively occurred at at the stage of perceptual identification and response choice. Therefore, it is likely to conclude that the modulation of attention over spatial Stroop effect was mediated by object-based attention, and this modulation occurred at the stage perceptual identification. Considering that the previous studies mostly investigated the effects of attention captured by abrupt onset on the spatial Stroop effects, few studies investigated the effects of attention captured by offset cue on the spatial Stroop effects. The aim of the fifth chapter was to investigate the role of attention induced by offset and abrupt onset cue in the spatial Stroop task. These results showed that attention elicited by offset cue or abrupt onset cue modulated the spatial Stroop effects, which reconciled with event integration account.