CN101290356A - Auxiliary method for dynamite source exploration - Google Patents

Auxiliary method for dynamite source exploration Download PDF

Info

Publication number
CN101290356A
CN101290356A CNA2007100655785A CN200710065578A CN101290356A CN 101290356 A CN101290356 A CN 101290356A CN A2007100655785 A CNA2007100655785 A CN A2007100655785A CN 200710065578 A CN200710065578 A CN 200710065578A CN 101290356 A CN101290356 A CN 101290356A
Authority
CN
China
Prior art keywords
partiald
wave
reconstruction
field
tau
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.)
Pending
Application number
CNA2007100655785A
Other languages
Chinese (zh)
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.)
BGP Inc
Original Assignee
BGP Inc
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 BGP Inc filed Critical BGP Inc
Priority to CNA2007100655785A priority Critical patent/CN101290356A/en
Publication of CN101290356A publication Critical patent/CN101290356A/en
Pending legal-status Critical Current

Links

Images

Abstract

The invention relates to a geophysical prospecting technology which is an auxiliary method for prospecting explosive source. The concrete steps are as follows: reconstruction parameters during the reconstruction of a wave field are acquired according to an actual field observation system; velocity and density parameters of longitudinal and transverse wave layers corresponding to an underground geological structure are obtained through the early seismic data processing, a conventional method is used to establish a discrete strata property model; wave field output of the position where a detection point is positioned at each moment is taken as a wave field reconstruction record of the moment to obtain stacked and migrated sections; and the stacked and migrated sections of the wave field reconstruction are compared with stacked and migrated sections obtained by processing the acquired data. The technology reduces the interference of man-made transverse waves in the record obtained by the wave field reconstruction, ensures that the seismic record of the reconstruction more conforms to the actual situation, is applicable to the wave field reconstruction of any completed structure, has the characteristics of high accuracy, high efficiency and low frequency dispersion, and ensures that the result of the oil-gas reservoir exploration is more reliable.

Description

A kind of householder method of dynamite source exploration
Technical field
The present invention relates to seismic exploration technique, specifically be a kind of under staggered-mesh the explosive source of high precision, high-level efficiency and low dispersion excite the householder method of the dynamite source exploration of elasticity wave-field reconstruction.
Background technology
In seismic exploration, by the seismic event that adopts focus (explosive source, vibroseis, weight whacker) to excite on ground, the reflective information that utilizes wave detector to receive underground geologic bodies obtains seismologic record, by a series of seismic data data processing with explain the various attributes obtain Depth Domain, and then definite oil reservoir distributes.At present, rebuild the reflective information (wave-field reconstruction) of geologic body according to the open-air earthquake data acquisition process of known formation property parameters emulation, and then with wave-field reconstruction to reflective information and the method that compares of the reflective information of open-air actual acquisition, can be more accurately hydrocarbon reservoir structure under definitely.Underground geologic bodies reflective information wave-field reconstruction and the process of descending hydrocarbon reservoir structure definitely describes for geologic interpretation, hydrocarbon reservoir structure and the hydrocarbon-bearing pool exploitation plays a very important role.The process energy true reappearance seismic event of wave-field reconstruction has very important meaning in the overall process of underground propagation to architectonic understanding.When the seismologic record that utilizes various layer attribute parameters reconstruction explosive sources to excite, describing exciting of explosive source is the gordian technique of wave-field reconstruction.
At present, the effect that excites of describing explosive source is not very good, mainly is when exciting compressional wave, can not guarantee the stressed shear wave that evenly produces of all directions.The size of the energy of ripple on all directions is also inconsistent in addition, and the size of excitation energy is often relevant with sizing grid, and wavefront neither be just round.In this respect, Liu Yanqiang (geophysical prospecting for oil, 1994, supplementary issue 2:1-18) a kind of method is disclosed, the top of Fig. 2 has provided the result, but do not break away from original mode, still approach explosive source (as Fig. 1), be called p ripple focus, p-wave source or single source by on the symmetric points of a plurality of centers, applying outside velocity of displacement.Prior art has been selected the focus item has been added in the vector-velocity of displacement with directivity, produced focus in the uneven problem of all directions size, the appearance of shear wave has influenced staggered-mesh precision, efficient, the elasticity wave field method for reconstructing that this explosive source excites can not accurately descend hydrocarbon reservoir structure definitely.
Summary of the invention
The object of the invention is to provide high precision under a kind of staggered-mesh, high-level efficiency and accurately descends the householder method of the dynamite source exploration of hydrocarbon reservoir structure definitely.
Concrete implementation step of the present invention:
(1) according to the reconstruction parameter in the recording geometry collection wave-field reconstruction process of open-air reality;
(reconstruction parameter comprises shot interval, track pitch, big gun number, road number, sampling interval, road length, cautious X, the Y coordinate of big gun)
(2), set up discrete stratum attribute model with conventional method through the interval velocity of ripple in length and breadth and the density parameter that conform to subsurface geologic structures that early stage, the seismic data data processing obtained;
(data is handled and is comprised that the single shot record to collecting carries out a series of processing such as static correction, normal moveout correction, velocity analysis, stack and skew)
(3) according to recording geometry p-and s-wave velocity field and density field, adopt following formula (1) to rebuild each wave field value on each discrete grid block point constantly, and with each geophone station position wave field output constantly as the wave-field reconstruction record in this moment;
∂ v x ∂ t = B ( ∂ τ xx ∂ x + ∂ τ xz ∂ z ) ∂ v z ∂ t = B ( ∂ τ xz ∂ x + ∂ τ zz ∂ z ) ∂ τ xx ∂ t = ( λ + 2 μ ) ∂ v x ∂ x + λ ∂ v z ∂ z + δ * F ∂ τ zz ∂ t = ( λ + 2 μ ) ∂ v z ∂ z + λ ∂ v x ∂ x + δ * F ∂ τ xz ∂ t = μ ( ∂ v x ∂ z + ∂ v z ∂ x ) - - - ( 1 )
Wherein: B=1/ ρ, λ+2 μ=ρ v p 2, μ=ρ v s 2
F is a source function
ρ, v p, v sBe respectively Media density, velocity of longitudinal wave, shear wave velocity
τ Xx, τ Zz, τ XzBe respectively normal stress and shearing stress
v x, v zBe respectively the line direction and the particle displacement speed of line direction downward vertically
(4) record that will obtain is done and the identical processing of step (2), obtains stack and migrated section;
(5) wave-field reconstruction stack and migrated section and image data are handled stack and migrated section contrast, the hydrocarbon reservoir structure under definitely that obtains.
The effect of invention
The present invention has been added in the tensor-normal stress that does not have directivity with the focus item, has then fundamentally solved focus in the uneven problem of all directions size, has avoided following the appearance of shear wave.Because normal stress appears at the same node (center position of Fig. 1) in the discrete model, focus item and discrete grid block are irrelevant, and its size can not change along with the variation of sizing grid.
The present invention is an explosive source wave-field reconstruction and descend the very effective method of hydrocarbon reservoir structure definitely under the two-dimensional case, has wave-field reconstruction that is applicable to any complex structure and the ability of descending hydrocarbon reservoir structure definitely.Have precision height, efficient height and the low characteristics of frequency dispersion.
The present invention has reduced the interference of artificial shear wave in the record that wave-field reconstruction obtains, the seismologic record of reconstruction is more tallied with the actual situation, and makes the achievement of reservoir exploration more reliable.
The reconstruction effect of the wave field that explosive source excites in the even two-dimension elastic medium of the present invention is Fig. 2, Fig. 2 A is the reconstruction result of prior art, have the weak jamming shear wave to occur in compressional wave wavefront surface inside, Fig. 2 B is the reconstruction result of this method, follows the shear wave that excites compressional wave to produce to disappear.Because this method ripple is added in the explosive source item in the normal stress equation, no longer has directivity as can be seen from Figure 3, wavefront energy even and for just round.Because the explosive source item only appears on the point, and is more approaching with actual conditions, make that the reconstruction wave field is more accurate.Excitation energy size is irrelevant with the discrete unit sizing grid in addition, makes the easier control of excitation energy.
Description of drawings
Fig. 1 is a focal point distribution plan of the present invention;
Fig. 2 A is a velocity of displacement wave field comparison diagram of the present invention;
Fig. 2 B is a velocity of displacement wave field comparison diagram of the present invention;
Fig. 3 is normal stress wave field figure of the present invention;
Fig. 4 is an embodiment of the invention velocity of longitudinal wave field sectional view;
The seismic cross-section that Fig. 5 embodiment of the invention is synthetic;
The actual seismic cross-section of Fig. 6 embodiment of the invention.
Embodiment
The present invention finds the solution two-dimentional wave-field reconstruction formula (1) under the single order staggered-mesh, to obtain the vertical component of wave field record.Adopt the staggered-mesh method to find the solution in the process of wave-field reconstruction, promptly find the solution two velocity of displacements earlier, find the solution the partial derivative equation of three stress again the time to the partial derivative equation of time.This two solution procedure that goes on foot that iterates is also exported each velocity of displacement vertical component of position, geophone station position constantly.
The specific embodiment of the present invention is:
(1) the open-air actual recording geometry acquisition parameter of foundation is determined the reconstruction parameter in the wave-field reconstruction process, (comprising shot interval, track pitch, big gun number, road number, sampling interval, road length, cautious X, the Y coordinate of big gun).
(2) read the P-wave interval velocity, S-wave interval velocity and the density that conform to this area's subsurface geologic structures that obtains through seismic data data processing in early stage (comprising that the single shot record to collecting carries out a series of processing such as static correction, normal moveout correction, velocity analysis, stack and skew).
(3) carry out interpolation according to recording geometry parameter and layer attribute (p-and s-wave velocity field and density field) and set up the discrete layer attribute of adaptation wave-field reconstruction: velocity of longitudinal wave field (Fig. 4), S-wave interval velocity field and density field.
(4) according to formula (1) and discrete wave field value of rebuilding on each each net point of the moment: velocity of displacement field and stress.With each constantly the speed component at place, geophone station position and the components of stress as this wave-field reconstruction record of each component constantly.
(5) processing identical with real data (comprising a series of processing such as static correction, normal moveout correction, velocity analysis, stack and skew) done in the wave-field reconstruction record that previous step is obtained, and obtains corresponding section (comprising stack and migrated section) (Fig. 5).
(6) section that wave-field reconstruction is obtained (stack and migrated section) (Fig. 6) simply contrasts with the section (stack and migrated section) that the field data collection conventional processing obtains, and can clearly differentiate underground hydrocarbon reservoir structure (arrow indication position among Fig. 6).

Claims (3)

1, a kind of householder method of dynamite source exploration is characterized in that concrete implementation step is as follows:
(1) according to the reconstruction parameter in the recording geometry collection wave-field reconstruction process of open-air reality;
(2), set up discrete stratum attribute model with conventional method through the interval velocity of ripple in length and breadth and the density parameter that conform to subsurface geologic structures that early stage, the seismic data data processing obtained;
(3) according to recording geometry p-and s-wave velocity field and density field, adopt following formula (1) to rebuild each wave field value on each discrete grid block point constantly, and with each geophone station position wave field output constantly as the wave-field reconstruction record in this moment;
∂ v x ∂ t = B ( ∂ τ xx ∂ x + ∂ τ xz ∂ z ) ∂ v z ∂ t = B ( ∂ τ xz ∂ x + ∂ τ zz ∂ z ) ∂ τ xx ∂ t = ( λ + 2 μ ) ∂ v x ∂ x + λ ∂ v z ∂ z + δ * F ∂ τ zz ∂ t = ( λ + 2 μ ) ∂ v z ∂ z + λ ∂ v x ∂ x + δ * F ∂ τ xz ∂ t = μ ( ∂ v x ∂ z + ∂ v z ∂ x ) - - - ( 1 )
Wherein: B=1/ ρ, λ+2 μ=ρ v p 2, μ=ρ v s 2
F is a source function
ρ, v p, v sBe respectively Media density, velocity of longitudinal wave, shear wave velocity
τ Xx, τ Zz, τ XzBe respectively normal stress and shearing stress
v x, v zBe respectively the line direction and the particle displacement speed of line direction downward vertically
(4) record that will obtain is done and the identical processing of step (2), obtains stack and migrated section;
(5) wave-field reconstruction stack and migrated section and image data are handled stack and migrated section contrast, the hydrocarbon reservoir structure under definitely that obtains.
2, the householder method of a kind of dynamite source exploration according to claim 1 is characterized in that the reconstruction parameter of step (1) comprises shot interval, track pitch, big gun number, road number, sampling interval, road length, cautious X, the Y coordinate of big gun.
3, the householder method of a kind of dynamite source exploration according to claim 1 is characterized in that the processing of step (2) data comprises that the single shot record to collecting carries out static correction, normal moveout correction, velocity analysis, stack and migration processing.
CNA2007100655785A 2007-04-17 2007-04-17 Auxiliary method for dynamite source exploration Pending CN101290356A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNA2007100655785A CN101290356A (en) 2007-04-17 2007-04-17 Auxiliary method for dynamite source exploration

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNA2007100655785A CN101290356A (en) 2007-04-17 2007-04-17 Auxiliary method for dynamite source exploration

Publications (1)

Publication Number Publication Date
CN101290356A true CN101290356A (en) 2008-10-22

Family

ID=40034703

Family Applications (1)

Application Number Title Priority Date Filing Date
CNA2007100655785A Pending CN101290356A (en) 2007-04-17 2007-04-17 Auxiliary method for dynamite source exploration

Country Status (1)

Country Link
CN (1) CN101290356A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104516014A (en) * 2013-09-27 2015-04-15 中国石油天然气集团公司 Fitting landform based wave field reconstruction method
CN106772597A (en) * 2016-12-12 2017-05-31 中国矿业大学 Take man-made explosion as the method for vibration signal inverting coal mine work area stress distribution
CN112505763A (en) * 2020-10-30 2021-03-16 中国石油天然气集团有限公司 Transverse wave seismic data quality detection method and system

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104516014A (en) * 2013-09-27 2015-04-15 中国石油天然气集团公司 Fitting landform based wave field reconstruction method
CN104516014B (en) * 2013-09-27 2018-04-03 中国石油天然气集团公司 A kind of wave field reconstructing method based on fitting landform
CN106772597A (en) * 2016-12-12 2017-05-31 中国矿业大学 Take man-made explosion as the method for vibration signal inverting coal mine work area stress distribution
CN112505763A (en) * 2020-10-30 2021-03-16 中国石油天然气集团有限公司 Transverse wave seismic data quality detection method and system

Similar Documents

Publication Publication Date Title
CN104977618B (en) A kind of method evaluated shale gas reservoir and find dessert area
Regone et al. Geologic model building in SEAM Phase II—Land seismic challenges
CN104570125B (en) A kind of method utilizing well data to improve image taking speed model accuracy
CN102053263B (en) Method for inspecting surface structure
CN105629325B (en) Foreland basin alluvial fan meticulous depiction and Forecasting Methodology
Yang et al. Locating and monitoring microseismicity, hydraulic fracture and earthquake rupture using elastic time-reversal imaging
CN101290355A (en) Method for determining carbonate rock hydrocarbon reservoir structure
Rovelli et al. Assessment of potential strong ground motions in the city of Rome
Tselentis et al. High-resolution passive seismic tomography for 3D velocity, Poisson’s ratio ν, and P-wave quality QP in the Delvina hydrocarbon field, southern Albania
CN102269823A (en) Wave field reconstruction method based on model segmentation
CN101173988B (en) Method for confirming underground oil-gas reservoir construction
CN101290356A (en) Auxiliary method for dynamite source exploration
Giustiniani et al. Crustal structure of central sicily
McBride et al. Visualizing Precambrian basement tectonics beneath a carbon capture and storage site, Illinois Basin
CN104345337B (en) A kind of time control reservoir parameter modeling method for seismic inversion
Abd El-Aal et al. Near-surface seismic refraction applied to exploring subsurface clay layer at a new mining area in southeast Cairo, Egypt
Mandal et al. Detecting sub-basalt Mesozoic sediments and active tectonics in the Narmada-Tapti Rift Zone, central India
Li et al. Time-Lapse Seismic Inversion for Predicting Reservoir Parameters Based on a Two-Stage Dual Network
Burschil et al. Finite-difference modelling to evaluate seismic P-wave and shear-wave field data
Wang et al. Geotechnical site characterization with 3D ambient noise tomography
Chinwuko et al. Integration of seismic attributes in delineation of channel features in Rence field of Niger Delta, Nigeria
Strayer et al. Upper crustal seismic velocity structure of the Hayward fault zone, San Francisco Bay, California, USA: Results from the 2016 East Bay Seismic Experiment (EBSI-16)
Cheruiyot Analysis of oil exploration plays in Anza Basin from a trap-type perspective
Velasquez-Espejo 3D multicomponent seismic characterization of a clastic reservoir in the Middle Magdalena Valley Basin, Colombia
Yeluru Enhancing Usability of the Multi-channel Analysis of Surface Wave (MASW) Technique for Subsurface Physical Property Mapping by Incorporating Random-Array Seismic Acquisition

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Open date: 20081022