CN110579803A - earthquake interpretation method for identifying biological reef body - Google Patents

earthquake interpretation method for identifying biological reef body Download PDF

Info

Publication number
CN110579803A
CN110579803A CN201910999333.2A CN201910999333A CN110579803A CN 110579803 A CN110579803 A CN 110579803A CN 201910999333 A CN201910999333 A CN 201910999333A CN 110579803 A CN110579803 A CN 110579803A
Authority
CN
China
Prior art keywords
reef
abnormal
seismic
biological
amplitude
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
CN201910999333.2A
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 CN201910999333.2A priority Critical patent/CN110579803A/en
Publication of CN110579803A publication Critical patent/CN110579803A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30Analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30Analysis
    • G01V1/303Analysis for determining velocity profiles or travel times
    • G01V1/305Travel times
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30Analysis
    • G01V1/307Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/622Velocity, density or impedance
    • G01V2210/6222Velocity; travel time
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/63Seismic attributes, e.g. amplitude, polarity, instant phase
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/63Seismic attributes, e.g. amplitude, polarity, instant phase
    • G01V2210/632Amplitude variation versus offset or angle of incidence [AVA, AVO, AVI]

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

the invention relates to the technical field of oil exploration, in particular to an earthquake interpretation method for identifying a biological reef body, which comprises the following steps: the abnormal reflector depicting technology is characterized in that a one-dimensional plane point segment representing an abnormal reflector of a biological reef is precisely depicted by using various seismic attributes; an abnormal time difference identification technology, wherein the abnormal time difference identification technology utilizes biological reef top surface reflection and earthquake TC2marking the characteristic of low abnormal time difference between the reflecting layers, performing variable speed mapping, and preliminarily realizing the two-dimensional planar distribution of the reefs. On the basis of comprehensive well logging-geology-seismic data analysis, a main body utilizes a seismic profile to conclude and determine the seismic response characteristics of the reef bodies in Otao Cujiki period in the Erdos basin, and then a set of three-anomaly identification method is explored, so that a good effect is achieved in the application and well position deployment of actual data. The application has good effect.

Description

Earthquake interpretation method for identifying biological reef body
Technical Field
The invention relates to the technical field of oil exploration, in particular to an earthquake interpretation method for identifying a biological reef body.
Background
Exploration practices show that the reef flat reservoir stratum has high abundance and high capacity, the reserve size accounts for more than 10% of the total reserve volume found by oil and gas in the world, and the reef flat reservoir stratum is always the target of key research in the field of oil and gas exploration. At present, the reef flat exploration in China is mainly concentrated in southern areas, the prediction of reef flat reservoir beds with sediments at the edge of the Clalmori stage of the Ordovician Otto system in the ancient kingdom under the Orldos basin is in the initial recognition stage, and no mature reference technical means exists.
In the period of Otaojicri Molii in the Ordors basin, the basin receives a large area of marine invasion. After experiencing sea level rising, not only the sea water cloth cover is expanded, the water body is deepened, but also the stratum has obvious difference, some places are raised and some places are sunken due to the structure movement. The water at the ridge is shallow, the illumination condition is good, and the method is suitable for the growth and the propagation of organisms and the development of the organism reefs. Referring to fig. 1, the controlled feature of positive landform of the reef determines the geophysical response characteristics of the reef to the upper marker layer with low time difference and reef-to-reef overlap. In addition, the biological reef lithology is mainly limestone, and the later period can receive dolomitic transformation, but the velocity difference between the limestone and the dolomitic is small, so that an obvious wave impedance interface is difficult to form, and the interior of the reef has the characteristic of weak amplitude.
In combination with the background of formation of the biosphere at the edge of the kreim stage and the corresponding characteristics of the biosphere on the seismic section, the present invention summarizes five characteristics of the biosphere, and please refer to fig. 2: the integral display of hill-shaped reflection characteristics; secondly, obvious overlapping contact relation can be seen on two sides of the reef body; thirdly, the time difference between the top of the reef body and the main marker layer of the basin is reduced by the earthquake TC2 reflection layer; fourthly, the interior of the reef body is blank-chaotic reflection; the interior of the reef body is in the characteristic of weak amplitude reflection;
Research and analysis recognize that the main technical means for current reef body reservoir prediction is to identify the reflection characteristics thereof by using seismic profiles, and generally comprises direct and indirect identification methods. The direct identification is mainly to distinguish and depict the reef body from the external form and the internal reflection characteristics of the reef body; the indirect identification mainly utilizes the phenomena of coating, overspread and the like existing above the reef body development position to further analyze the speed difference and indirectly predict the distribution and development of the reef body reservoir. In any identification method, the seismic profile is the comprehensive reflection characteristic of the underground rock stratum and the geologic body, and has certain multi-solution, so that the disposition of the well position in the oil and gas exploration process is obviously still insufficient. The identification of the visible reef bodies is a relatively difficult matter, and then it is more difficult to find the effective reservoir of the reef bodies rich in oil and gas reservoirs. On the basis of comprehensive well logging-geology-seismic data analysis, a main body utilizes a seismic profile to conclude and determine the seismic response characteristics of Oerdos basin Ordovic Kilmori period reef bodies, a set of three-anomaly identification method is further explored, and good effects are achieved in the application of actual data and well position deployment. The application has good effect.
Disclosure of Invention
In order to solve the above problems, the present invention provides an earthquake interpretation method for identifying a biological reef body.
The technical problem to be solved by the invention is realized by adopting the following technical scheme:
a seismic interpretation method of identifying a biological reef body, the method comprising the steps of:
Step 1, an abnormal reflector depicting technology, wherein the abnormal reflector depicting technology is used for precisely depicting a one-dimensional plane point segment representing an abnormal reflector of a biological reef by using various seismic attributes;
Step 2, an abnormal time difference identification technology, wherein the abnormal time difference identification technology is implemented by utilizing biological reef top surface reflection and earthquake TC2Marking the characteristic of low abnormal time difference between the reflecting layers, performing variable speed mapping, and preliminarily realizing two-dimensional planar distribution of the reefs;
and 3, an abnormal amplitude prediction technology, wherein the abnormal amplitude prediction technology is used for counting the earthquake amplitude change rule on the basis of the variable-speed imaging, finely predicting an abnormal low-amplitude area and finishing fine distribution prediction of a target body plane.
further, the seismic attributes refer to specific metrics derived from the seismic data regarding geometry, kinematics, dynamics and statistical properties. Aiming at the reef-shaped appearance of the living beings, geometric attributes and physical attributes can be used for describing; wherein the content of the first and second substances,
On a texture attribute section in geometric attributes, the hill-shaped form of the Ordovician biological reef is obvious in performance, reef cores and reef wings are easy to distinguish, meanwhile, instantaneous phase attributes in the geometric attributes are integrated, the external form of a hill-shaped body is further implemented, meanwhile, the reflection characteristic of the overlarge from two sides to the reef cores is also very obvious, and a good foundation is laid for implementing one-dimensional point sections of the abnormal reflector plane of the biological reef.
Further, step 2 comprises:
step 201, tracing and explaining an earthquake TC2 reflecting layer and a reef top reflecting layer on a seismic section;
step 202, changing the two reflecting layers into graphs and performing subtraction operation on the two layers of structural graphs;
and 203, obtaining a plane rule with the reduced thickness of the reflecting layer through grid operation on the basis of the planar one-dimensional point section of the implemented biological reef abnormal reflector, and determining the two-dimensional planar distribution and the thickness of the reef.
Further, step 3 comprises:
step 301, extracting an amplitude attribute profile of a seismic profile target layer;
Step 302, counting the abnormal amplitude;
Step 303, thinning the amplitude abnormal area of the reef two-dimensional plane distribution area,
And step 304, defining the reef plane distribution rule.
compared with the prior art, the invention has the beneficial effects that: on the basis of comprehensive well logging-geology-seismic data analysis, a main body utilizes a seismic profile to conclude and determine the seismic response characteristics of Oerdos basin Ordovic Kilmori period reef bodies, a set of three-anomaly identification method is further explored, and good effects are achieved in the application of actual data and well position deployment. The application has good effect.
drawings
FIG. 1 is a cross-section of pre-stack time migration of the survey line S097079;
FIG. 2 is a pattern diagram of a biological reef;
FIG. 3 is a cross section of the H135629 line instantaneous cosine and H157259 line instantaneous phase attribute;
FIG. 4 is a pre-stack time migration profile of H135629 and 07HL5229 lines;
Fig. 5 shows the pre-stack time, average amplitude attribute profile and average amplitude statistics for the S097079 survey line.
Detailed Description
the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The present invention provides a preferred embodiment:
On the basis of defining five response characteristics of the geophysical of the biological reef, three abnormal identification methods of the seismic of the biological reef, namely an abnormal reflector depicting technology, an abnormal time difference identification technology and an abnormal amplitude prediction technology, are summarized;
Specifically, one-dimensional plane point segments representing the abnormal reflector of the biological reef are finely depicted by using various seismic attributes; by using biological reef top surface reflection and earthquake TC2marking the characteristic of low abnormal time difference between the reflecting layers, performing variable speed mapping, and preliminarily realizing two-dimensional planar distribution of the reefs; and finally, on the basis of the variable-speed mapping, counting the earthquake amplitude change rule, and finely predicting an abnormal low-amplitude area to finish fine spread prediction of the target body plane.
The biological reef is in a hill-like reflection form on the seismic section and has obvious difference with the normal sedimentary stratum reflection, and by using the characteristic, the seismic attribute analysis technology is developed to implement the point section distribution rule representing the abnormal reflection of the biological reef.
Seismic attributes refer to special measures derived from seismic data regarding geometry, kinematics, dynamics and statistical properties, which can be characterized for a biosystock-like shape using geometric and physical attributes.
Referring to fig. 3, on the texture attribute section in the geometric attribute, the mound-like form of the aotao biological reef is obvious, and the reef core and the reef wing are easily distinguished. Meanwhile, instantaneous phase attributes in the physical attributes are integrated, the external form of the dune is further implemented, meanwhile, the reflection characteristics of the reef core from two sides are very obvious, and a good foundation is laid for implementing one-dimensional point sections of the plane of the abnormal reflector of the biological reef.
The reef has positive landform characteristics, top reflection and upper earthquake TC2the reflection layer is significantly reduced in cross section compared to both sides. Analysis shows that the LT1 belongs to the top reflecting layer of the reef and the earthquake TC2the reflection layer has a minimum time difference of only 50ms, and the reef wing reflection layer and the earthquake TC2The reflection layer time difference is substantially 150 ms.
Referring to FIG. 4, the top, wing and T of the reef are utilizedC2Predicting the plane position of the reef by the characteristic that the time difference of the reflecting layers is reduced, firstly, tracing and explaining an earthquake TC2 reflecting layer and a reef top reflecting layer on an earthquake section, secondly, performing variable-speed mapping on the two reflecting layers and performing subtraction operation on the two layers of structural diagrams, and finally, obtaining a plane rule with the reduced thickness of the reflecting layers and the two-dimensional plane distribution and the thickness of the reef on the basis of the one-dimensional point section of the plane of the implemented biological reef abnormal reflector through grid operation.
referring to fig. 5, the impedance difference between the limestone of the biological reef in the krimo mile period and the limestone and the dolomite is small, the amplitude on the seismic section is not obviously changed, the interior of the reef body is blank-disordered reflection, and the reef body is compared with surrounding rocks on two sides on the seismic section and has the characteristic of weak amplitude, so that the amplitude abnormal area of the two-dimensional plane distribution area of the reef body is refined by extracting the amplitude attribute section of the target layer of the seismic section and counting the abnormal amplitude, and the plane distribution rule of the reef body is determined.
in the present invention, unless otherwise expressly stated or limited, the terms "mounted," "connected," "secured," and the like are to be construed broadly and can, for example, be fixedly connected, detachably connected, or integrally formed; can be mechanically or electrically connected; they may be directly connected or indirectly connected through intervening media, or they may be connected internally or in any other suitable relationship, unless expressly stated otherwise. The specific meanings of the above terms in the present invention can be understood by those skilled in the art according to specific situations.
In the present invention, unless otherwise expressly stated or limited, the first feature "on" or "under" the second feature may be directly contacting the first and second features or indirectly contacting the first and second features through an intermediate. Also, a first feature "on," "over," and "above" a second feature may be directly or diagonally above the second feature, or may simply indicate that the first feature is at a higher level than the second feature. A first feature being "under," "below," and "beneath" a second feature may be directly under or obliquely under the first feature, or may simply mean that the first feature is at a lesser elevation than the second feature.
while the invention has been described above with reference to an embodiment, various modifications may be made and equivalents may be substituted for elements thereof without departing from the scope of the invention. In particular, the various features of the embodiments disclosed herein may be used in any combination, provided that there is no structural conflict, and the combinations are not exhaustively described in this specification merely for the sake of brevity and conservation of resources. Therefore, it is intended that the invention not be limited to the particular embodiments disclosed, but that the invention will include all embodiments falling within the scope of the appended claims.

Claims (4)

1. A seismic interpretation method for identifying a biological reef body is characterized by comprising the following steps:
Step 1, an abnormal reflector depicting technology, wherein the abnormal reflector depicting technology is used for precisely depicting a one-dimensional plane point segment representing an abnormal reflector of a biological reef by using various seismic attributes;
Step 2, an abnormal time difference identification technology, wherein the abnormal time difference identification technology is implemented by utilizing biological reef top surface reflection and earthquake TC2Marking the characteristic of low abnormal time difference between the reflecting layers, performing variable speed mapping, and preliminarily realizing two-dimensional planar distribution of the reefs;
and 3, an abnormal amplitude prediction technology, wherein the abnormal amplitude prediction technology is used for counting the earthquake amplitude change rule on the basis of the variable-speed imaging, finely predicting an abnormal low-amplitude area and finishing fine distribution prediction of a target body plane.
2. The method of claim 1, wherein the seismic attributes are specific metrics derived from seismic data regarding geometric, kinematic, kinetic and statistical properties. Aiming at the reef-shaped appearance of the living beings, geometric attributes and physical attributes can be used for describing; wherein the content of the first and second substances,
On a texture attribute section in geometric attributes, the hill-shaped form of the Ordovician biological reef is obvious in performance, reef cores and reef wings are easy to distinguish, meanwhile, instantaneous phase attributes in the geometric attributes are integrated, the external form of a hill-shaped body is further implemented, meanwhile, the reflection characteristic of the overlarge from two sides to the reef cores is also very obvious, and a good foundation is laid for implementing one-dimensional point sections of the abnormal reflector plane of the biological reef.
3. the seismic interpretation method of identifying the biological reef bar as recited in claim 2, wherein the step 2 comprises:
Step 201, tracing and explaining an earthquake TC2 reflecting layer and a reef top reflecting layer on a seismic section;
Step 202, changing the two reflecting layers into graphs and performing subtraction operation on the two layers of structural graphs;
and 203, obtaining a plane rule with the reduced thickness of the reflecting layer through grid operation on the basis of the planar one-dimensional point section of the implemented biological reef abnormal reflector, and determining the two-dimensional planar distribution and the thickness of the reef.
4. the seismic interpretation method of identifying the biological reef body of claim 2, wherein step 3 comprises:
step 301, extracting an amplitude attribute profile of a seismic profile target layer;
Step 302, counting the abnormal amplitude;
step 303, refining an amplitude abnormal area of the reef two-dimensional plane distribution area;
and step 304, defining the reef plane distribution rule.
CN201910999333.2A 2019-10-21 2019-10-21 earthquake interpretation method for identifying biological reef body Pending CN110579803A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910999333.2A CN110579803A (en) 2019-10-21 2019-10-21 earthquake interpretation method for identifying biological reef body

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910999333.2A CN110579803A (en) 2019-10-21 2019-10-21 earthquake interpretation method for identifying biological reef body

Publications (1)

Publication Number Publication Date
CN110579803A true CN110579803A (en) 2019-12-17

Family

ID=68815160

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910999333.2A Pending CN110579803A (en) 2019-10-21 2019-10-21 earthquake interpretation method for identifying biological reef body

Country Status (1)

Country Link
CN (1) CN110579803A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112462423A (en) * 2020-11-03 2021-03-09 中国石油天然气集团有限公司 Method and device for predicting thickness of biological reef

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104977611A (en) * 2014-04-08 2015-10-14 中国石油化工股份有限公司 Reef reservoir engraving method
CN107831554A (en) * 2017-10-25 2018-03-23 中国石油化工股份有限公司 A kind of complex biological reef microfacies identification and decision method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104977611A (en) * 2014-04-08 2015-10-14 中国石油化工股份有限公司 Reef reservoir engraving method
CN107831554A (en) * 2017-10-25 2018-03-23 中国石油化工股份有限公司 A kind of complex biological reef microfacies identification and decision method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
尚振广: ""开江—梁平海槽东侧长兴组生物礁沉积特征及有利区分析"", 《中国优秀硕士学位论文全文数据库•基础科学辑》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112462423A (en) * 2020-11-03 2021-03-09 中国石油天然气集团有限公司 Method and device for predicting thickness of biological reef

Similar Documents

Publication Publication Date Title
CN101408624B (en) Forecasting and evaluating technologies of three-dimensional earthquake optimum time window river course sand body storage layer
CN103529475B (en) A kind of method identified and explain carbonate paleokarst reservoir three dimensional structure
CN101980053B (en) Complicated reef flat reservoir predicting method
CN105259572B (en) The seismic facies computational methods classified automatically based on seismic multi-attribute parametrical nonlinearity
CN109655896B (en) Method and system for identifying and depicting karst ancient river channel
CN103675913B (en) Detection method for fractured reservoirs
CN105093307B (en) The oblique true formation thickness acquiring method of Lower Paleozoic strata
CN105319585B (en) A kind of method hidden using thin-layers interference amplitude recovery identification of hydrocarbon
CN104698496A (en) Small sand body boundary identification method and small sand body space quantitative description method
CN106443781A (en) Method for predicting distribution of thin sand body
CN109541685A (en) A kind of river channel sand recognition methods
CN104965225A (en) Wavelet sidelobe information-based oil reservoir identification method
CN104977611B (en) A kind of reef reservoir engraving process
CN110927817B (en) Method for optimally designing and controlling thin-layer superposed riverway horizontal well track
CN110579803A (en) earthquake interpretation method for identifying biological reef body
CN106483562A (en) A kind of method based on seismic amplitude predictably layer thickness variation
CN107703548B (en) Shallow stratum stratum boundary division methods based on the deposit qualities factor and RL return loss grade curve peak valley
CN114265116B (en) Quantitative statistical method and device for channel and ridge amplitude of granite buried hill weathering crust
CN104991271A (en) Prestack depth migration method based on VTI medium bidirectional wave field superposition
CN109164486A (en) A kind of identification and evaluation method and system of favorable oil/gas target area
Sheasby Depositional patterns of the Upper Devonian Waterways Formation, Swan Hills area, Alberta
CN114545498A (en) Comprehensive depicting method for ancient karst fracture-cave structure
CN104318106B (en) Method for predicting gas content and productivity by attenuation of frequency division energy
CN111562615B (en) Method for determining sand body development area of retrograde beach dam, and method and device for identifying oil reservoir
CN111948716A (en) Method for calculating depth of submarine water channel by using seismic data

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20191217