利用台阵面波资料研究华北克拉通北缘燕山褶皱带及其邻区的地壳与上地幔结构

The continent of eastern China, especially the North China Craton (NCC), has endured intensive tectonic renovation during Mesozoic and Cenozoic, with the presence of widespread magmatism, high heat flow and development of large sedimentary basins and mountain ranges. The cratonic lithosphere of the region has been destroyed remarkably, which is characterized by not only a significant reduction in thickness but also complex modifications in physical and chemical properties of the lithosphere. As for the tectonic regime controlling the evolution of the NCC, various models have been put forward, including the impingement of mantle plumes (“mushroom cloud” model), the collision of south China block and north China block, the subduction of the Pacific plate, etc. Lithosphere delamination and thermal erosion were proposed as the two end-member mechanisms of the lithospheric thinning. However, given the paucity of the data, deep structural evidence is currently still scarce for distinguishing and testifying these models. To better understand the deep structure of the NCC, from 2000 to the present, temporary seismic array observations have been conducted in the NCC by the Seismological Laboratory of the Institute of the Geology and Geophysics, Chinese Academy of Sciences under the North China Interior Structure Project (NCISP). Many arrays extend from the North China Craton and the off-craton regions, and traverse a lot of main tectonic boundaries. A total of more than 300 broadband seismic stations have been deployed along several profiles that traversed the major tectonic units within the craton’s interior, at the boundary areas and in the neighboring off-craton regions. These stations recorded abundant high-quality data, which provides an unprecedented opportunity for us to unravel the deep structural features of the NCC using seismological methods. Among all the seismological methods, the surface wave method appears to be an efficient and widely adopted technique in studying the crustal and upper mantle structures. In particular, it can provide the absolute values of S-wave velocity that are difficult to obtain with other methods. Benefiting from the deployment of dense seismic arrays, progresses have been made in improving the spatial resolution of surface wave imaging, which makes it possible to resolve the fine-scale velocity structures of the crust and upper mantle based on surface wave analysis. Meanwhile, the differences in the S-wave velocities derived from Rayleigh and Love wave data can provide information on the radial anisotropy beneath the seismic arrays. In this thesis, using the NCISP-III broadband data and based on phase velocity dispersion analysis and inversion of fundamental mode Rayleigh and Love waves, I investigated the lateral variations in the S-wave velocity structure of the crust and uppermost mantle beneath the Yanshan Belt and adjacent regions at the northeastern boundary of the NCC. Based on the constructed structural images, I discussed possible deep processes of the craton destruction in the study region.

利用接收函数方位变化研究地壳各向异性

Because of its sensitivity to the velocity discontinuity of the earth, receiver function technique has become a routine procedure used to probe interior structure of the earth. Receiver functions contain anisotropic information of the earth’s interior, however, traditional receiver function techniques such as migration imaging and waveform inversion method, which are based on isotropic media assumption, can not effectively extract the anisotropy information contained in the azimuth variation pattern. Only by using the anisotropic media, e.g. a model with symmetric axis of arbitrary orientation, computing the response, can we obtain the detailed anisotropy information hidden in the radial and transversal receiver function. Focusing on the receiver function variation pattern changing wtih different back azimuths, we introduced different kinds of symmetric systems of seismic anisotropy used often, and summarized some possible causes of anisotropy formation. We show details about how to calculate the response of a stratified anisotropy model with symmetric axis of arbitrary orientation. We also simulated receiver functions among different models and analyzed how the changing of anisotropic parameters influence the azimuth variation pattern of receiver functions. The anisotropy study by receiver function analysis was applied to Taihang Mountain Range (TMR) in North China in this thesis. The maximum entropy spectrum deconvolution technique was used to extract radial and transversal receiver functions from the waveforms of 20 portable seismic stations deployed in TMR. Considering the signal-to-noise ratio and the azimuth coverage, we got the variation pattern of receiver functions for 11 stations. After carefully analyzing the pattern of the receiver functions that we got, we obtained the reliable evidence on the existence of anisotropy in the shallow crust in TMR. Our results show that, although the thickness of the upper crustal layer is only about 1 km, the layer shows a strong anisotropy with magnitude of 8~15%; in the deeper of crust, the magnitudes of anisotropy is about 3%~5%, showing a pattern with fast-symmetric-axis. The crust anisotropy beneath TMR in North China obtained in this study also shows a significant difference in both the lateral and vertical scale, which might imply a regional anisotropy characteristic in the studied region.

消除采集脚印的处理方法研究

Seismic technique is in the leading position for discovering oil and gas trap and searching for reserves throughout the course of oil and gas exploration. It needs high quality of seismic processed data, not only required exact spatial position, but also the true information of amplitude and AVO attribute and velocity. Acquisition footprint has an impact on highly precision and best quality of imaging and analysis of AVO attribute and velocity. Acquisition footprint is a new conception of describing seismic noise in 3-D exploration. It is not easy to understand the acquisition footprint. This paper begins with forward modeling seismic data from the simple sound wave model, then processes it and discusses the cause for producing the acquisition footprint. It agreed that the recording geometry is the main cause which leads to the distribution asymmetry of coverage and offset and azimuth in different grid cells. It summarizes the characters and description methods and analysis acquisition footprint’s influence on data geology interpretation and the analysis of seismic attribute and velocity. The data reconstruct based on Fourier transform is the main method at present for non uniform data interpolation and extrapolate, but this method always is an inverse problem with bad condition. Tikhonov regularization strategy which includes a priori information on class of solution in search can reduce the computation difficulty duo to discrete kernel condition disadvantage and scarcity of the number of observations. The method is quiet statistical, which does not require the selection of regularization parameter; and hence it has appropriate inversion coefficient. The result of programming and tentat-ive calculation verifies the acquisition footprint can be removed through prestack data reconstruct. This paper applies migration to the processing method of removing the acquisition footprint. The fundamental principle and algorithms are surveyed, seismic traces are weighted according to the area which occupied by seismic trace in different source-receiver distances. Adopting grid method in stead of accounting the area of Voroni map can reduce difficulty of calculation the weight. The result of processing the model data and actual seismic demonstrate, incorporating a weighting scheme based on the relative area that is associated with each input trace with respect to its neighbors acts to minimize the artifacts caused by irregular acquisition geometry.

KEL-TI介质的地震波场响应与成像方法研究

To deal with the problems in multi-component converted seismic wave exploration in coal fields, the wave propagating features and imaging methods of multi-component converted waves in coal measure strata are researched in this thesis firstly. The relations between viscoelasticity and anisotropy in coal measure strata are analyzed to build KEL-TI model, and which seismic wave propagating and attenuating features are researched. The disadvantages of converted wave imaging methods based on common converted point gather are analyzed and constant velocity no NMO converted wave imaging method based on common scattering point gather is put forward, according to Huygens-Fresnel principle, which applicabilities in the elastic isotropic, elastic TI and KEL-TI situations are discussed. To different model simulation data, the common scattering point gathers’ and stacked profiles’ features are analyzed. The results show that the method can image compressional waves and converted waves with high precision. Secondly, the resolution enhancing theories and methods of converted wave are researched by Rayleigh wave suppressing, converted wave static correction and poststack inverse-Q filtering. 1) The polarization filter is designed by the instantaneous polarization information of seismic waves, and the Rayleigh wave suppressing method is researched. From the spectrum analysis before and after filtering, it can be derived that the amplitudes are kept relatively. 2) To constant velocity no NMO converted wave imaging method, the static correction method based on common equivalent offset point gather is put forward and tested to the actual converted waves. 3) The relation between equivalent quality factor of converted wave, compressional wave quality factor and the ratio of compressional to shear wave velocity is derived. The compressional wave quality inversion method by first arrivals of none-offset VSP is researched, and which is then transformed to the equivalent quality factor to perform inverse-Q filtering of actual converted waves. The result has shown that the method can recover the high frequency energy of converted waves. At last, the theories and methods researched in this thesis are practiced to the 3D3C seismic exploration in Guqiao coal mine in Huainan and achieve good results.

煤系地层三维三分量地震数据解释方法研究

Conventional 3D seismic exploration cannot meet the demand of high yield and high efficiency safe production in coal mine any more. Now it is urgent to improve the discovery degree of coal mine geological structures for coal production in China. Based on 3D3C seismic exploration data, multi-component seismic information is fully excavated. First systematic research on 3D3C seismic data interpretation of coal measure strata is carried out. Firstly, by analyzing the coal measure strata, the seismic-geologic model of coal measure strata is built. Shear wave logging is built by using regression analysis. Horizon calibration methods of PP-wave and PS-wave are studied and the multi-wave data are used together to interpret small faults. Using main amplitude analysis technology, small faults which cannot be found from PP-wave sections can be interpreted from the low frequency PS-wave sections. Thus, the purpose to applying PS-wave data to fine structure assistant interpretation is achieved. Secondly, PP- and PS-wave post-stack well constrained inversion methods of coal measure strata are studied. Joint PP- and PS-wave post-stack inversion flow is established. More attribute parameters, which are applied in fine lithology interpretation of coal measure strata, are obtained from combinations of the inversion results. Exploring the relation between rock with negative Poisson’s ratio and anisotropy, fracture development in coal seam are predicted. Petrophysical features of coal measure strata are studied, and the relations between elastic parameters and lithology, fluid and physical properties are established. Inversions of the physical parameters such as porosity, permeability and water saturation, which reflect lithology and fluid property, are obtained. Finally, the approaches of shear wave splitting and Thomsen parameters inversion, which provide new ideas for seismic anisotropy interpretation of coal measure strata, are studied to predict fracture development. The results of practical application indicate that the methods in this paper have good feasibility and applicability. They have positive significance for high yield and high efficiency safe production in coal mine.

南襄盆地泌阳凹陷断裂斜坡带构造特征与油气成藏

The faulted slope zone of Biyang depression, a multiple hydrocarbon accumulation zone lying in a rich oil depression of Nanxiang basin, is a structural-sedimentary compounded slope, which is developed in Yanshanian period and has an area of 500 km2. From the ‘bottom up’, the developed strata may be divided into Yuhuangding formation in Neogene, Dacang Fang, Hetao-yuan plus Liaozhuang formations and Fenghuangzheng plus Pingyuan formation in Neogene, while Hetao-yuan formation is the main hydrocarbon-bearing target. Because of transtensional stress fields formed by persistent action of large-scale faulting in the south of the depression, sedimentary differential compaction in different stages, and tectonic inversion in later developing stage of the depression, a series of nose structure zones cut by different strike faults are developed. Therefore, the reservoir migration and accumulation are controlled by the complex faulted-nose structural zone, reservoir types are dominated by faulted-noses, faulted-blocks and fault-lithology, while lithology and stratigraphic unconformable reservoirs are locally developed. In combination with demands of practical production, applying with a new research approach of systematology and a combination with dynamic and static modes, guided by modern petroleum geologic theory, and based on previous data and studies, new techniques, methods of geophysical exploration, various computer simulation and forecasting techniques are applied in the new research of this paper. Starting from the structural features and formation mechanism, the forming mechanism of faulted structure, conditions and controlling factors of hydrocarbon accumulation, as well as various space-time allocation relationships in the process of accumulation are analyzed in the research. Besides that, the hydrocarbon migration, accumulation mechanism and dynamic evolution process are also discussed in the paper. Through the research, the accumulation rule of the faulted slope zone in faulted lake basin, the distribution and enrichment regularity of different reservoir controlling factors are systematically summarized. The summarizations indicate that the faulted slope is a favorable orientational zone, hydrocarbon is accumulated in nose structures and enriched in the main body of nose structures, faulted transformation zone and the ascent direction of laddering faulted blocks, the faults are the controlling factors, hydrocarbon accumulation zones controlled by fault-lithology are distributed along the faulting direcion. In the end, hydrocarbon migration and accumulation models of complex faulted-nose blocks are established. 1) Down cut model—‘flank-sheet’: the hydrocarbon is migrated like ‘sheet’ along a series of faults with parallel distribution and accumulated in the flank of nose structures; 2）Cross cut --‘axis-string’ model: the hydrocarbon cutting across the faults is migrated like ‘string’ and accumulated in the axis of nose structures. In view of different distribution models, reservoir forming combination patterns are divided and hydrocarbon reservoir evaluation exploration is carried out, which achieves good results in application. Key words: faulted slope zone; migration and accumulation model; reservoir controlling mechanism; reservoir-forming combination

松辽盆地南部腰英台区块低渗透储层整体压裂改造技术研究

The Yaoyingtai Block is located within the northeastern Changling Depression of southern Songliao Basin, where the reservoir sandstones are petrophysically characterized by very low permeability, which results in the low success probability of artificial fracturing, and the low oil yield by water injection in the course of oil production. In order to improve the situations as stated above, this research aims to work out an integral fracturing technology and strategy applicable to the low permeable reservoirs in Yaoyingtai Block. Under the guidance of geological theory, reservoir engineering and technology, the subsurface occurrences of natural and hydraulic fractures in the reservoirs are expected to be delineated, and appropriate fracturing fluids and proppants are to be optimized, based on the data of drilling, well logging, laboratory and field experiments, and geological data. These approaches lay the basis of the integral fracturing technology suitable for the low permeable reservoir in the study area. Based on core sample test, in-situ stress analysis of well logging, and forward and inversion stress field modeling, as well as fluid dynamic analysis, the maximum in-situ stress field is unraveled to be extended nearly along the E-W direction (clustering along N85-135°E) as is demonstrated by the E-W trending tensional fractures. Hydraulic fractures are distributed approximately along the E-W direction as well. Faulting activities could have exerted obvious influences on the distribution of fractures, which were preferentially developed along fault zones. Based on reservoir sensitivity analysis, integrated with studies on rock mechanics, in-situ stress, natural fracture distribution and production in injection-production pilot area, the influences of primary fractures on fracturing operation are analyzed, and a diagnostic technology for primary fractures during depressurization is accordingly developed. An appropriate fracturing fluid (hydroxypropyl guar gum) and a proppant (Yixing ceramsite, with a moderate-density, 0.45-0.9mm in size) applicable to Qingshankou Formation reservoir are worked out through extensive optimization analysis. The fracturing fluid can decrease the damage to the oil reservoir, and the friction in fracturing operation, improving the effect of fracturing operation. Some problems, such as sand-out at early stage and low success rate of fracturing operations, have been effectively solved, through pre-fracturing formation evaluation, “suspension plug” fracturing, real-time monitoring and limited-flow fracturing. Through analysis of fracture-bearing tight reservoir with variable densities and dynamic analysis of influences of well patterns on fracturing by using numerical simulation, a fracturing operation scheme for the Qingshankou Formation reservoir is proposed here as being better to compress the short factures, rather than to compress the long fractures during hydraulic fracturing. It is suggested to adopt the 450m×150m inverted 9-spot well pattern in a diamond shape with wells placed parallel to fractures and a half fracture length of 60-75m.

青藏高原东部长周期大地电磁测深研究

As a complement to conventional MT, Long-period Magnetotellurics (LMT) has been developed at low frequency for soundings of deep electric structures. Eastern Himalayan Syntaxis (EHS) and surrounding area is a key place for the study of dynamics of the uplift of Tibetan plateau. Experiments in the pioneer studies for EHS3D project showed that the study area shares an unusual low resistive crust and upper mantle. Conventional MT could not provide sufficient information about the deep structures of the study area that requested long period MT measurement to be complemented. This thesis presents the LMT studies in eastern Tibet along the EHS3D-3 Profile from Xiachayu to Yushu including data acquisition, processing, inversion and interpretation. The effective period of the measured LMT signals extend from 10s up to 30000s for the duration more than one week measurement. The resulting model shows that the LMT sounding coincides with the MT data in overlapped periods. Especially the induction arrows and tippers derived from LMT data provide more information about the base of the conductors beneath the plateau with higher resolution. Anomalous induction coefficients and 2-D model suggest extensive conductive bodies beneath Lhasa block and Qiangtang terrain which would be a possible evidence for partial melt and fluids at depth.

青藏高原东部电性结构特征及其构造意义

Eastern Himalayan Syntaxis (EHS) and its surroundings (eastern margin of Tibet) is one of the most complicated tectonic areas in the world. As the exhaust opening of the balanced materials of the Tibetan Plateau during the collision of Indan and Eurasian plates, the deep structure beneath EHS surrounding region is referred to as the key to the study of the dynamics of the plateau. EHS3D project, sponsored by NSFC, has been proposed to explore the deep electric features of the area. During the first stage of EHS3D(2006-2008), MT+LMT measurements have been conducted along two lines from Chayu to Qingshuihe (EHS3D-3) and Chayu to Ruoergai (EHS3D-2). This paper will discuss the MT models of EHS3D-3 line. By the data procrssing, including distortion analysis, Robust estimation and strike decomposition, rotated apparent resitivities and phases have been obtained for each station. Then conventional 2-D inversion algorithms (NLCG and RRI) were employed to produce 2-D models. The final preferred 2-D model suggests that the upper crust consists of resistive blocks while in mid-lower crust there are two extensive conductive bodies beneath Lhasa block and Qiangtang terrain respectively. Jinshajiang suture is a gradient belt and Bangong-Nujiang suture appear a conductive belt dipping to the north. . We concluded that the formation of the two conductive bodies attributed to the partial melt and fluids in the lower crust. The regional electric strike derived from decomposition analysis indicates that the crust and upper mantle move in different manners. The upper crust moves like slips of rigid blocks along major slip faults while the lower crust creeps as a flow in the conductive channels.

华北及周边地区层析成像研究

Our study deals with the high resolution body wave tomography in North china and adjacent areas(30°N-43°N，100°E-130°E), where earthquakes occurred many times in history and has a very complicated geological structure. 6870 events recorded at 273 digital seismic stations from CDSN during 1996－2002 and stations settled by Seislab of IGCAS in Bohai Bay area, including 1382 local earthquakes and 5488 teleseismic earthquakes are used in this study. In the data we used, the average number of received stations is greater than 5, the error of picking up direct arrival time is 0.1－0.5s. Before the inversion, we use Checkerboard method to confirm the reliability of result of Local events; use Restoring Resolution Test to confirm the reliability of result of teleseismic events. We also analyzed the effect of different parameters in the inversion. Based the analysis above, the model used in this paper is divided into small blocks with a dimension of 0.33°in the latitude and longitude directions and 5km、15km、30km in depth, and initial velocity model. Using pseudobending method to calculate the ray traveling path, LSQR algorithm to inverse, finally, we got the body velocity images below 25km and above 480km in this area using Joint- inversion with local events and teleseismic events. We made the conclusion at last: (1)at top zone of the south of Sichuan Basin , there exits low velocity anomalies, below 40km is the high velocity zone extend to 300km; (2) Above the 40km of Ordos block exits low velocity zone, while below 40km until 240km, the high velocity anomalies are interlaced by low velocity anomalies. Below 300km, the anomalies are unclear any more; (3) On the whole, the velocity structure below 400km on the mantle transition zone of Eastern China area shows its changes from low velocity to high velocity.

微动勘查方法在煤矿采区勘探中的应用研究

The conventional microtremor survey is based on the single point of exploration, which includes collecting field data，estimating the phase velocity, investing the dispersion curve and obtaining the S–wave velocity structure. In the case of large-scale exploration, and when making the two-dimensional velocity section, the inversion is quite time-consuming, laborious and its precision depends on the subjective interpretation, which makes the results differently from person to person. In fact，we do not need the S-wave velocity values but only need the relative variation of velocity. For these reasons, this paper is desired to calculate the apparent S-wave velocity (Vx) to replace the S-wave velocity inversion and to obtain the relative variation of the S-wave velocity. Using this method, we can decrease the analysts’ effect, shorten the data processing time and improve work efficiency. The apparent S-wave velocity is a variable of the surface wave property, which can clearly reflect the downcast columns, mined-out areas and other unusual geological bodies. In this paper, Matlab is used to establish the three-dimensional data volume of the apparent S-wave velocity, from which we can get any apparent S-wave velocity section we need. Through the application case, the designed method is proved to be reliable and effective. The downcast columns, mined-out areas and other unusual geological bodies can be clearly showed in the apparent S-wave velocity section. And from the contour of the apparent S-wave velocity, the interface shape of the major target layers can be controlled basically.

井间地震速度层析成像方法与软件系统

Crosshole Seismic tomography has been broadly studied and applied in the fields of resource exploration and engineering exploration because of its special observing manner and better resolution than normal seismic exploration. This thesis will state the theory and method of Crosshole Seismic tomography. Basing on the previous studies，the thesis studied the initial velocity model，ray-tracing method, and developed the three-dimension tomography software. All the cells that a ray passes through are of the same velocities if the paths from transmitters to receivers are straight. The cells that the each ray passes through are recorded, and rays that pass through each cell are calculated. The ray average velocity which passes through a cell is set as the cell velocity. Analogously we can make a initial node velocity model because the velocity sum is calculated on the all cells which own to a certain node, and the cell number is summed about each nodes，the ratio of the velocity sum to the all cells number is set as the node velocity. The inversion result from the initial node velocity model is better than that of the average velocity model. Ray-bending and Shortest Path for Rays (SPR) have shortcomings and limitations respectively. Using crooked rays obtained from SPR rather than straight lines as the starting point can not only avoid ray bending converging to the local minimum travel time path, but also settle the no smooth ray problem obtained by SPR. The hybrid method costs much computation time, which is roughly equal to the time that SPR expends. The Delphi development tool based on the Object Pascal language standard has an advantage of object-oriented. TDTOM (Three Dimensions Tomography) was developed by using Delphi from the DOS version. Improvement on the part of inversion was made, which bring faster convergence velocity. TDTOM can be used to do velocity tomography from the first arrival travel time of the seismic wave, and it has the good qualities of friendly user interface and convenient operation. TDTOM is used to reconstruct the velocity image for a set of crosshole data from Karamay Oil Field. The geological explanation is then given by comparing the inversion effects of different ray-tracing methods. High velocity zones mean the cover of oil reservoir, and low velocity zones correspond to the reservoir or the steam flooding layer.

旋转轴对称弱各向异性介质模型的射线追踪方法

Ray tracing is a rapid and effective method for wave field calculation. Not only in the field of seismic-wave theory, but also in the field of seismic inversion and migration imaging,the seismic ray tracing method has become one of the most important methods. In anisotropic media, group velocity and phase velocity have different propagation directions. The seismic wave propagates along the direction of group velocity , it does not depend on the direction of phase velocity. Ray angle is a complex function with respect to phase angle, it is difficult to measure and calculate. But most rocks are weak anisotropic, so the expression of phase velocity can be simplified greatly. Based on the approximate expression of phase velocity this thesis for rotating axisymmetric weak anisotropic media deduces an expression of the partial derivative of phase velocity and an expression of group velocity with the method of linear approximation. This paper uses the fourth order Runge-Kutta method together with the two-dimensional interpolation and linear interpolation to obtain the parameters of the physical locations. At last the paths of seismic wave in rotating axisymmetric weak anisotropic media are computed. According to the analysis of the computational results, it indicates that the method developed in this paper has strong adaptability, high computational efficiency and high accuracy for rotating axisymmetric weak anisotropic media.

贝尔凹陷苏德尔特构造带古潜山构造与储层研究

With the develop ment of oil and gas exploration, the conventional struc ture exploration era has gradually been substituted by the concealed reser voir exploration technology. Hill poll becomes one of the most important areas in the future exploration. This paper is based on the three-dimensional seismic interpretation of Sudeerte structure. In terms of the overall character istics of Sudeerte structure, we use the coherent cube and the time slice to interpret the fault in plane. For the interpretation，we especially used the well to adjust the interpretation results. The results of seismic attribution analysis, spectrum decomposition and post- stack seismic inversion forecast that hill pools reservoir are dist ributed in several bands along the north-northeast to northeast-east. Xing'anling Group shows that the potential reser voirs are mainly distributed along Bei 14 –Bei 40 and De (99-212) - Bei (16-1) and Budate Group distributed along Bei 14 –Bei 40 in northwestern direction and De (99-212) - Bei (16-1) in north eastern direction. At the same time, by analyzing the structure and the reservoir, and combining with other data, three models are built. The characteristics of reservoirs dist ribution are concl uded，and potent ional favorable exploration dire ctions are predi cted.

数据驱动型多次波衰减方法的研究

In the practical seismic profile multiple reflections tend to impede the task of even the experienced interpreter in deducing information from the reflection data. Surface multiples are usually much stronger, more broadband, and more of a problem than internal multiples because the reflection coefficient at the water surface is much larger than the reflection coefficients found in the subsurface. For this reason most attempts to remove multiples from marine data focus on surface multiples, as will I. A surface-related multiple attenuation method can be formulated as an iterative procedure. In this essay a fully data-driven approach which is called MPI —multiple prediction through inversion (Wang, 2003) is applied to a real marine seismic data example. This is a pretty promising scheme for predicting a relative accurate multiple model by updating the multiple model iteratively, as we usually do in a linearized inverse problem. The prominent characteristic of MPI method lie in that it eliminate the need for an explicit surface operator which means it can model the multiple wavefield without any knowledge of surface and subsurface structures even a source signature. Another key feature of this scheme is that it can predict multiples not only in time but also in phase and in amplitude domain. According to the real data experiments it is shown that this scheme for multiple prediction can be made very efficient if a good initial estimate of the multiple-free data set can be provided in the first iteration. In the other core step which is multiple subtraction we use an expanded multi-channel matching filter to fulfil this aim. Compared to a normal multichannel matching filter where an original seismic trace is matched by a group of multiple-model traces, in EMCM filter a seismic trace is matched by not only a group of the ordinary multiple-model traces but also their adjoints generated mathematically. The adjoints of a multiple-model trace include its first derivative, its Hilbert transform and the derivative of the Hilbert transform. The third chapter of the thesis is the application for the real data using the previous methods we put forward from which we can obviously find the effectivity and prospect of the value in use. For this specific case I have done three group experiments to test the effectiveness of MPI method, compare different subtraction results with fixed filter length but different window length, invest the influence of the initial subtraction result for MPI method. In terms of the real data application, we do fine that the initial demultiple estimate take on a great deal of influence for the MPI method. Then two approaches are introduced to refine the intial demultiple estimate which are first arrival and masking filter respectively. In the last part some conclusions are drawn in terms of the previous results I have got.