CN1158539C - Optimizing design method of 3D seismic observation system based on geologic geophysical model - Google Patents

Optimizing design method of 3D seismic observation system based on geologic geophysical model Download PDF

Info

Publication number
CN1158539C
CN1158539C CNB011204451A CN01120445A CN1158539C CN 1158539 C CN1158539 C CN 1158539C CN B011204451 A CNB011204451 A CN B011204451A CN 01120445 A CN01120445 A CN 01120445A CN 1158539 C CN1158539 C CN 1158539C
Authority
CN
China
Prior art keywords
earthquake
seismic
imaging results
geophysical model
geologic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CNB011204451A
Other languages
Chinese (zh)
Other versions
CN1337583A (en
Inventor
牟永光
狄帮让
魏建新
顾培成
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CNB011204451A priority Critical patent/CN1158539C/en
Publication of CN1337583A publication Critical patent/CN1337583A/en
Application granted granted Critical
Publication of CN1158539C publication Critical patent/CN1158539C/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Abstract

The present invention relates to a method of optimization design of a three-dimensional earthquake observing method based on a geologic geophysical model. The method comprises the steps that (A) a priori geologic geophysical model is established; (B) two or more types of observation systems of three-dimensional earthquake data acquisition are designed; (C) data acquisition is completed by the combined forward modeling of three-dimensional wave equation numerical simulation of a virtual spectral method and three-dimensional earthquake geophysical model simulation of a method; (D) depth migration before stacking is carried out on acquired data on the basis of a common reflective point, and some imaging results of earthquake treatment are obtained; (E) the imaging qualities of the imaging results of the earthquake treatment are comprehensively compared and evaluated, and one optimal earthquake imaging result is determined through comparison. The method has the advantages of wide range of application, and good applied effect, and the method has operation significance of guiding earthquake prospecting actual production of oil and natural gas, and enhancing oil-gas prospecting successful rate.

Description

Optimizing design method of 3 D seismic observation system based on geologic geophysical model
Technical field
The present invention relates to a kind of technical design optimization method of 3-d seismic exploration in petroleum gas and other mineral resources and the Geological Engineering, particularly a kind of optimizing design method of 3 D seismic observation system based on geologic geophysical model.
Background technology
Seismic prospecting three big links are to gather, handle and explain.Collection is the basis; The receiving system that activating system that epicenter excitation point position is formed on the ground and wave detector acceptance point position are formed must be by correct geometric distributions, and promptly the seismic observation system of You Huaing just can collect effectively complete subsurface geological structure information.Therefore, recording geometry plays an important role in seismic prospecting, and its design quality is directly connected to the success or failure of exploration task.Because the 3-D seismics workload is big, investment is high, three-dimensional work area often is arranged in the complex structure area again, so the stereo observing system optimal design is important doubly just.
So far, both at home and abroad mostly the stereo observing system design software that widely adopts is based on horizontal layer geologic model and common midpoint (CMP) superposition principle, and they are not suitable for the 3-D seismics technical design on complex structure, high spud angle stratum.Under the complex structure situation, use this class design proposal, not only can not obtain high-quality seismic imaging, even the seismic imaging result that may make mistake.
Summary of the invention
Therefore, purpose of the present invention makes provides a kind of optimizing design method of 3 D seismic observation system based on geologic geophysical model.
This method is based on the three-dimensional geological and the geophysical model of exploratory area priori, adopt numerical simulation and physical simulation to unite positive artistic skills art, according to common reflection point (CRP) principle, be imaged as evaluation criterion with the earthquake of geologic objective high-quality, the image quality of three-dimensional seismic acquisition techniques design proposal and geological data processing is contacted directly, realized optimizing stereo observing system with this and design this acquisition technique target.
Optimizing design method of 3 D seismic observation system based on geologic geophysical model provided by the invention may further comprise the steps:
A) according to the existing exploration achievement in exploratory area, promptly complex geological structure and complicated geological volume morphing and geophysical parameters thereof are set up the geologic geophysical model of priori; B) be exploration targets with above-mentioned geologic geophysical model, design the recording geometry that two or more 3D seismic datas are gathered; C) with this model as the 3-D seismics object of observation, design desired technical parameter and construction figure in the open-air earthquake data acquisition process of laboratory simulation according to recording geometry, obtain the open-air original seismic data of simulation, the concrete practice is to adopt the simulation of the three-dimensional fluctuation of empty spectrometry equation number value just drilling with two kinds of method associatings of 3-D seismics physical model simulation to finish this data collection task; D) data of gathering are carried out 3D seismic data and handle,,, obtain earthquake and handle several imaging results the data preceding depth shift processing that change wherein based on common reflection point; E) with the geologic geophysical model of priori as estimating and contrast standard, with signal to noise ratio (S/N ratio), resolution and structural feature feature as the technical Analysis main contents, the comprehensive comparative evaluation that several imaging results are carried out image quality is handled in earthquake, through contrast, determine one of them best earthquake imaging results.
Step e wherein) also comprise:
E1) if imaging results does not also reach the target call of seismic prospecting, then revise original stereo observing system parameter, again carry out the work of 3D seismic data acquisition process subsequently, described stereo observing system parameter comprises 3-D seismics observation width degree, bin size, covering or illumination number of times, illuminance, seismic trace sum, maximum offset, wherein observes width degree N be defined as:
N = θ 2 π ( C 1 L c L i + C 2 n c n i )
Wherein: L c/ L iThe horizontal vertical ratio of-arrangement sheet,
n c/ n iThe horizontal vertical ratio of-degree of covering,
C 1-with arrange the relevant coefficient of the horizontal vertical ratio of sheet,
C 2-with the relevant coefficient of the horizontal vertical ratio of degree of covering, and C 2=1-C 1, generally get C 1=C 2=0.5,
θ-center, subarea bin azimuthal distribution scope (radian);
E2), then export 3 D seismic observation system optimal design process flow diagram if imaging results has reached the target call of seismic prospecting.
The present invention is described in detail below in conjunction with accompanying drawing, so that the advantage of further understanding above-mentioned purpose of the present invention, characteristics and bringing thus.
Description of drawings
Fig. 1 is the process flow diagram of the optimizing design method of 3 D seismic observation system based on geologic geophysical model of the present invention;
Fig. 2 (a) is the Qianmiqiao buried hill structural physical model according to the priori data make;
Fig. 2 (b) is the two dimension 337 survey line log sheets that North and South direction is passed the top, buried hill;
Fig. 3 adopts the three-dimensional Wave equation forward modeling of empty spectrometry respectively and adopts the resulting seismologic record of three-dimensional seismic physical simulation on 337 surveys line;
Fig. 4 is the horizontal superposition imaging results of geological data (337 line) that three kinds of recording geometrys are gathered, and wherein Fig. 4 (a) is vertical 1 recording geometry of 24 horizontal strokes of 8 lines, 2 big guns; Fig. 4 (b) is vertical 2 recording geometrys of 24 horizontal strokes of 8 lines, 4 big guns; Fig. 4 (c) is vertical 2 recording geometrys of 24 horizontal strokes of 8 lines, 8 big guns;
Fig. 5 is the earthquake data offset imaging results (the dropping cut slice figure of cima top imaging) of different width degree three-dimensional acquisitions, and wherein, Fig. 5 (a) is the earthquake data offset of narrow three-dimensional acquisition, and Fig. 5 (b) is the earthquake data offset of wide three-dimensional acquisition.
Embodiment
It is target since 1996 at exploratory area concrete complex geological structure that physical prospecting key lab of China National Petroleum Corporation (CNPC) explored a kind of, based on the CRP seismic imaging, new three-dimensional acquisition technical design method, this method workflow as shown in Figure 1.
At first set up the geologic geophysical model of priori, set up the seismic velocity model of underground medium usually according to the existing exploration achievement in exploratory area (complex geological structure and complicated geological volume morphing and geophysical parameters thereof); Second step was an exploration targets with above-mentioned geologic geophysical model, designed the recording geometry that two or more 3D seismic datas are gathered; The 3rd goes on foot with this model as the 3-D seismics object of observation, design desired technical parameter and construction figure in the open-air earthquake data acquisition process of laboratory simulation according to recording geometry, this method characteristics are to adopt mathematical computations and two kinds of methods of Physical Experiment to simulate open-air earthquake construction, obtaining the open-air original seismic data (referring to Fig. 3) of simulation, specifically is to adopt the simulation of the three-dimensional fluctuation of empty spectrometry equation number value just drilling with two kinds of method associatings of 3-D seismics physical model simulation to finish this data collection task; The 4th step was the data of gathering to be carried out 3D seismic data handle, and handling principle is based on common reflection point (CRP), and core technology is repeatedly preceding depth shift processing, i.e. a kind of seismic imaging technology of current popular, thus obtaining earthquake handles several imaging results; The 5th step with the geologic geophysical model of priori as estimating and contrast standard, with signal to noise ratio (S/N ratio), resolution and structural feature feature as the technical Analysis main contents, the comprehensive comparative evaluation that several imaging results are carried out image quality is handled in earthquake, through contrast, determine one of them best earthquake imaging results; If this imaging results does not also reach the target call of seismic prospecting, then to revise original stereo observing system parameter (design proposal), carry out the work of 3D seismic data acquisition process subsequently again; If this imaging results has reached the target call of seismic prospecting, then its pairing stereo observing system is exactly the end product of this area's 3-d seismic exploration optimal design, exports 3 D seismic observation system optimal design workflow diagram subsequently.At this moment, whole optimal design work also declares to be finished simultaneously.
The 3 D seismic observation system main design contents is that three-dimensional observation template (several lines, a few process of preparing Chinese medicine) is selected and the stereo observing system parameter designing.The stereo observing system parameter mainly contains 3-D seismics observation width degree, bin size, covering or illumination number of times, illuminance, seismic trace sum, maximum offset etc.When definite 3-D seismics is observed width degree N, the following formula that utilizes the applicant to propose:
N = θ 2 π ( C 1 L c L i + C 2 n c n i )
Wherein: L c/ L IThe horizontal vertical ratio of-arrangement sheet,
n c/ n IThe horizontal vertical ratio of-degree of covering,
C 1-with arrange the relevant coefficient of the horizontal vertical ratio of sheet,
C 2-with the relevant coefficient of the horizontal vertical ratio of degree of covering, and C 2=1-C 1, generally get C 1=C 2=0.5,
θ-center, subarea bin azimuthal distribution scope (radian).
When: N<0.5 is narrow orientation recording geometry;
N 〉=0.5 o'clock is the wide-azimuth recording geometry;
N 〉=0.85 o'clock is the omnibearing observation system;
Oil-gas exploration task in conjunction with China the western and eastern complex area, physical prospecting key lab of China National Petroleum Corporation (CNPC) (CNPC) repeatedly is applied to this method in the actual 3-D seismics technical design, huge port Qianmiqiao buried hill model for example, the contrary zone of fracture model of covering in car row area, Xinjiang, Xinjiang Yi Qikelike fold thrust belts model, the complex fault block models of Jurong, Jiangsu and Hongze etc. have all been obtained good effect.Be example with huge port Qianmiqiao buried hill structure 3 D seismic observation system optimal design below, this method is done simple the introduction.
Fig. 2 (a) is the Qianmiqiao buried hill structural physical model according to the priori data make, and Fig. 2 (b) passes the two dimension 337 survey line log sheets at top, buried hill for North and South direction.Fig. 3 adopts the three-dimensional Wave equation forward modeling of empty spectrometry respectively and adopts the resulting seismologic record of three-dimensional seismic physical simulation on 337 surveys line, relatively these two records can think that both are similar.Adopted several different recording geometrys to carry out the 3D seismic data collecting work to this area,, obtained different seismic imaging result (seeing Fig. 4 and Fig. 5) by handling.Subsequently with geologic objective as standard, various imaging results are estimated, optimize the recording geometry design.For example the signal noise ratio (snr) of image of Fig. 4 (c) should think then that than the height of Fig. 4 (b) and Fig. 4 (a) the corresponding recording geometry of Fig. 4 (c) is wherein best; For example the top imaging of Fig. 5 (b) cima is relatively good again, and focusing quality is than Fig. 5 (b) height, and explanation should select the wide-azimuth stereo observing system to carry out the earthquake construction.Some above-mentioned practical applications show: this method scope of application is extensive, and effect is good, has realistic meaning to instructing petroleum gas seismic prospecting actual production, raising oil-gas exploration success ratio.

Claims (1)

1, a kind of optimizing design method of 3 D seismic observation system based on geologic geophysical model may further comprise the steps:
A) according to the existing exploration achievement in exploratory area, i.e. complex geological structure and complicated geological volume morphing and geophysical parameters thereof are set up the geologic geophysical model of priori;
B) be exploration targets with above-mentioned geologic geophysical model, design the recording geometry that two or more 3D seismic datas are gathered;
C) with this model as the 3-D seismics object of observation, design desired technical parameter and construction figure in the open-air earthquake data acquisition process of laboratory simulation according to recording geometry, obtain the open-air original seismic data of simulation, the concrete practice is to adopt the simulation of the three-dimensional fluctuation of empty spectrometry equation number value just drilling with two kinds of method associatings of 3-D seismics physical model simulation to finish this data collection task;
D) data of above-mentioned collection are carried out 3D seismic data and handle,,, obtain earthquake and handle several imaging results the data of the gathering preceding depth shift processing that change wherein based on common reflection point;
E) with the geologic geophysical model of priori as estimating and contrast standard, with signal to noise ratio (S/N ratio), resolution and structural feature feature as the technical Analysis main contents, the comprehensive comparative evaluation that several imaging results are carried out image quality is handled in described earthquake, through contrast, determine one of them best earthquake imaging results; Step e wherein) also comprise:
E1) if imaging results does not also reach the target call of seismic prospecting, then revise original stereo observing system parameter, again carry out the work of 3D seismic data acquisition process subsequently, described stereo observing system parameter comprises 3-D seismics observation width degree, bin size, covering or illumination number of times, illuminance, seismic trace sum, maximum offset, wherein observes width degree N be defined as:
N = θ 2 π ( C 1 L c L i + C 2 n c n i )
Wherein: Lc/Li-arranges the horizontal vertical ratio of sheet,
The horizontal vertical ratio of nc/ni-degree of covering,
C1-and the relevant coefficient of the arrangement horizontal vertical ratio of sheet,
C2-and the relevant coefficient of the horizontal vertical ratio of degree of covering, and C2=1-C1 generally get C1=C2=0.5,
θ-center, subarea bin azimuthal distribution scope (radian);
E2), then export 3 D seismic observation system optimal design process flow diagram if imaging results has reached the target call of seismic prospecting.
CNB011204451A 2001-07-13 2001-07-13 Optimizing design method of 3D seismic observation system based on geologic geophysical model Expired - Fee Related CN1158539C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB011204451A CN1158539C (en) 2001-07-13 2001-07-13 Optimizing design method of 3D seismic observation system based on geologic geophysical model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB011204451A CN1158539C (en) 2001-07-13 2001-07-13 Optimizing design method of 3D seismic observation system based on geologic geophysical model

Publications (2)

Publication Number Publication Date
CN1337583A CN1337583A (en) 2002-02-27
CN1158539C true CN1158539C (en) 2004-07-21

Family

ID=4664139

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB011204451A Expired - Fee Related CN1158539C (en) 2001-07-13 2001-07-13 Optimizing design method of 3D seismic observation system based on geologic geophysical model

Country Status (1)

Country Link
CN (1) CN1158539C (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102147480A (en) * 2011-01-12 2011-08-10 中国石油集团川庆钻探工程有限公司 Test method of three-dimensional observation system

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7286690B2 (en) * 2003-10-23 2007-10-23 Pgs Americas, Inc. Method for imaging of pre-stack seismic data
CN1299127C (en) * 2004-03-26 2007-02-07 中国石油天然气集团公司 Laminated dielectric double-focusing method for earthquake observation system optimized design
CN100429527C (en) * 2005-03-25 2008-10-29 中国石油天然气集团公司 Three-dimensional seismic data processing quality monitoring technology
CN101625417B (en) * 2008-07-08 2011-09-07 中国石油集团东方地球物理勘探有限责任公司 Method for optimizing design of vertical seismic profile observation system
CN102023310B (en) * 2009-09-19 2013-05-15 中国石油化工集团公司 Uniformity analysis method of three-dimensional observation system
CN102053271B (en) * 2010-11-17 2012-09-26 中国石油天然气股份有限公司 Earthquake collection quality analysis method and device
CN102169187B (en) * 2010-12-01 2013-01-23 中国科学院地质与地球物理研究所 Quantitative analysis method for overall balance of three dimensional seismic observation system
CN102495943B (en) * 2011-10-31 2015-04-22 国家卫星海洋应用中心 Modeling method for geophysical model
CN103116183B (en) * 2011-11-16 2015-05-27 中国石油天然气集团公司 Method of oil earthquake collection surface element covering degree property body slicing mapping
CN103777229A (en) * 2012-10-23 2014-05-07 中国石油化工股份有限公司 VSP observation system design method facing objective layer
CN103149590B (en) * 2013-02-26 2016-01-27 佟小龙 geophysical imaging method and device
CN104142518B (en) * 2013-10-29 2017-05-03 中国石油化工股份有限公司 Method for analyzing pre-stack time migration response of seismological observation system
CN104635259B (en) * 2013-11-06 2017-09-01 中国石油天然气集团公司 It is a kind of to determine the method and device that observation system arranges piece scope
CN103901465A (en) * 2014-03-21 2014-07-02 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Design method of holographic three-dimensional seismic prospecting and observing system
CN104050359B (en) * 2014-05-30 2016-10-05 中国石油大学(华东) A kind of the Forward Modeling based on stereo observing system arrangement sheet data segmentation
CN104090297A (en) * 2014-06-14 2014-10-08 吉林大学 Reverse illumination method for optimizing earthquake collection observing system
CN104536041B (en) * 2014-12-17 2017-02-22 中国石油天然气集团公司 Optimization method of seismological observation system parameters
EP3413092B1 (en) * 2017-06-08 2022-06-01 Total Se Method for evaluating a geophysical survey acquisition geometry over a region of interest, related process, system and computer program product
CN111723517B (en) * 2020-05-26 2022-09-02 河海大学 Geophysical model optimization method integrating machine learning principle and gradual deformation method
CN114722590B (en) * 2022-03-23 2024-03-22 中石化石油工程技术服务有限公司 Design optimization method of random acquisition observation system based on geophysical model

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102147480A (en) * 2011-01-12 2011-08-10 中国石油集团川庆钻探工程有限公司 Test method of three-dimensional observation system

Also Published As

Publication number Publication date
CN1337583A (en) 2002-02-27

Similar Documents

Publication Publication Date Title
CN1158539C (en) Optimizing design method of 3D seismic observation system based on geologic geophysical model
CN100557464C (en) Seismic prospecting horizon calibration method based on the prestack wave-field simulation
EP2250520B1 (en) Method for indexing a subsurface volume for the purpose of inferring geologic information
CN1170170C (en) Spectral decomposition for seismic interpretation
US8694261B1 (en) 3D-well log invention
CN102937720B (en) Well control improves the method for seismic data resolution
CN1189218A (en) Method and apparatus for seismic signal processing and exploration
CN104678434A (en) Method for predicting storage layer crack development parameters
CN106556869A (en) Method that tomography walk sliding displacement is portrayed under a kind of extention environment quantitatively
CN108445538A (en) The method and system of Depth Domain layer Q model is established based on reflected seismic information
CN104142518A (en) Method for analyzing pre-stack time migration response of seismological observation system
CN111722284B (en) Method for establishing speed depth model based on gather data
CN106777707B (en) Method for quantitatively identifying well logging lithology by using improved spider web diagram
CN101359056A (en) Method for generating longitudinal wave time domain high precision converted wave sectional plane
US20220221604A1 (en) Method of analysing seismic data
Martuganova et al. 3D deep geothermal reservoir imaging with wireline distributed acoustic sensing in two boreholes
CN107942379A (en) A kind of method for improving complex fault block rate pattern precision
CN104459798A (en) Velocity modeling method based on RTM imaging
Lafond et al. Migration of wide‐aperture onshore‐offshore seismic data, central California: Seismic images of late stage subduction
CN113031070B (en) Method for making depth domain synthetic seismic record
CN114609675A (en) Quantitative recovery method for carbonate rock stratum sedimentary micro-landform based on high-frequency cycle
CN104237939B (en) A kind of 3 D seismic observation system multi-frequency migration noise acquisition methods and device
AU2014394076B2 (en) Methods and systems for identifying and plugging subterranean conduits
CN110954956B (en) Method for evaluating acquisition trace of observation system and computer-readable storage medium
CN110531421B (en) Sedimentary microfacies identification method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C19 Lapse of patent right due to non-payment of the annual fee
CF01 Termination of patent right due to non-payment of annual fee