CN112462417B - Inversion structure zone plane range identification method based on inversion intensity coefficient - Google Patents

Inversion structure zone plane range identification method based on inversion intensity coefficient Download PDF

Info

Publication number
CN112462417B
CN112462417B CN202010493517.4A CN202010493517A CN112462417B CN 112462417 B CN112462417 B CN 112462417B CN 202010493517 A CN202010493517 A CN 202010493517A CN 112462417 B CN112462417 B CN 112462417B
Authority
CN
China
Prior art keywords
inversion
thickness
stratum
intensity
formula
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.)
Active
Application number
CN202010493517.4A
Other languages
Chinese (zh)
Other versions
CN112462417A (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 National Offshore Oil Corp CNOOC
CNOOC China Ltd Tianjin Branch
Original Assignee
China National Offshore Oil Corp CNOOC
CNOOC China Ltd Tianjin Branch
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 National Offshore Oil Corp CNOOC, CNOOC China Ltd Tianjin Branch filed Critical China National Offshore Oil Corp CNOOC
Priority to CN202010493517.4A priority Critical patent/CN112462417B/en
Publication of CN112462417A publication Critical patent/CN112462417A/en
Application granted granted Critical
Publication of CN112462417B publication Critical patent/CN112462417B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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. for interpretation or for event detection
    • G01V1/282Application of seismic models, synthetic seismograms
    • 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. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • 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/624Reservoir parameters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/66Subsurface modeling
    • G01V2210/661Model from sedimentation process modeling, e.g. from first principles

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)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

The invention discloses a method for identifying a plane range of an inversion structural band based on an inversion intensity coefficient, which comprises the following steps of S1 establishing an inversion mode of the inversion structural band, S2 identifying whether the inversion structural type accords with the mode established in the step S1 from the section of a research area, and being used as the basis for whether the method is applicable, S3 analyzing according to the inversion mode of the step S1, establishing an inversion intensity formula, S4 obtaining the thickness of a stratum of the research area, wherein the thickness of the stratum in the inversion period is delta H, S5, the denudation amount of the research area is recovered, the thickness of the denudation amount of the stratum in the inversion period is E, S6, substituting the thickness of the stratum obtained in the step S4 and the denudation amount obtained in the step S5 into the inversion intensity calculation formula established in the step S2, obtaining the maximum value and the minimum value of the inversion intensity, and S7 establishing an inversion intensity coefficient formula, substituting the maximum value and the minimum value of the inversion intensity obtained in the step S6 into the formula, obtaining a final dimensionless inversion intensity coefficient R', thereby identifying the extent of the inverted build band on the plane.

Description

Inversion structure band plane range identification method based on inversion intensity coefficient
Technical Field
The invention belongs to the technical field of petroleum seismic exploration, and particularly relates to a reverse structure zone plane range identification method based on a reverse intensity coefficient.
Background
The concept of inversion zone (abbreviated as inversion zone) is a composite structure with opposite properties to the former structure, which is generated by the reverse change of the tectonic action of the same geologic body in different geologic history periods (Williams, 1989; McClay, 1991; Mitra, 1993). Due to the good anticline shape, the oil-gas separation device is often a favorable place for oil-gas enrichment and is a hot spot area for oil exploration. Therefore, the identification of the planar range of the tectonic reversal zone in the exploration stage not only facilitates the analysis of the mechanism of cause of the tectonic reversal zone, but also facilitates the rapid selection of the exploration target area.
Conventionally, the inverted structure band can be recognized only on a cross section, and it is difficult to accurately recognize a planar range. The current plane identification method of the inversion structural belt mainly comprises a landform method, a thickness method and a hand-drawing method, and the principle is that a rough distribution range is qualitatively circled on a plane through section identification. The method has two defects that the boundary of a reverse band and a non-reverse band is difficult to distinguish on one hand, and the boundary of the reverse band is marked in a fuzzy way on the other hand. Therefore, how to effectively and accurately identify the planar range of the inversion structure belt is a problem which needs to be solved by those skilled in the art.
Disclosure of Invention
The invention aims to provide a method for identifying a planar range of an inversion structure belt based on an inversion intensity coefficient, so as to solve the problems of the background art.
In order to achieve the above object, a specific technical solution of the method for identifying a planar range of an inversion structure band based on an inversion intensity coefficient of the present invention is as follows:
a reverse structure belt plane range identification method based on a reverse strength coefficient comprises the following steps:
s1: establishing a reverse structure belt reverse mode;
s2: identifying from the section of the investigation region whether the inversion structure type corresponds to the pattern established in step S1 as a basis for whether the method is applicable;
s3: establishing a reversal strength formula according to the reversal pattern analysis of the step S1;
s4: the thickness of the stratum in the research area is obtained, and the thickness of the stratum in the reversal period is delta H;
s5: restoring the denudation amount of the research area, wherein the thickness of the stratum denudation amount in the reversal period is E;
s6: substituting the stratum thickness obtained in the step S4 and the ablation amount obtained in the step S5 into the inversion strength calculation formula established in the step S2 to obtain the maximum value and the minimum value of the inversion strength;
s7: and establishing a reversal intensity coefficient formula, normalizing the obtained reversal intensity, and substituting the maximum value and the minimum value of the reversal intensity obtained in the step S6 into the formula to obtain a final dimensionless reversal intensity coefficient R', so that the range of the reversal structural belt is identified on the plane.
Further, step S1 is to establish two inversion patterns of inversion structural bands, one is a single-break inversion structural band for single-fault control, and the other is a double-break inversion structural band for two-fault control, and both types of inversion structural bands are subjected to inversion action of early settlement and later lifting.
Further, a relational expression of the thickness of the stratum of the inversion structural zone and the buried depth of the bottom boundary of the stratum, namely an inversion strength calculation formula, is established:
Figure GDA0003604827200000021
wherein H Ta The geological significance of the buried depth of the stratum bottom boundary of the inversion structural zone is to reflect the degree of inversion deformation of the inversion structural zone.
The denudation amount recovery method adopts a stratum contrast method, the stratum deposition has inheritance and continuity in the same structural layer, according to the characteristic, the deposition thickness of an overlying layer can be estimated according to the product of the ratio of the saved complete adjacent layer thickness (overlying layer thickness/underlying layer thickness) and the underlying layer thickness, and if the estimated value is larger than the overlying residual layer thickness, the exceeding part is the denudation thickness; otherwise, the upper cladding layer is not denuded. According to the above analysis, the adjacent layer thickness ratio is:
Figure GDA0003604827200000022
wherein Ha is the thickness of the adjacent overlying layer that remains intact and Hb is the thickness of the adjacent underlying layer that remains intact. Further calculating the thickness of the stratum denudation as follows:
E=λ×H' b -H' a
Wherein E is the thickness of the denudation amount of the stratum, H 'b is the thickness of the overburden which is not completely stored, and H' a is the thickness of the underlayer which is completely stored.
Further, the formation thickness Δ H obtained in step S4 and the ablation thickness E obtained in step S5 are substituted into the inversion strength formula created in step S3, and the maximum value Rmax and the minimum value Rmin of the inversion strength are calculated.
Further, an inversion strength coefficient formula is established:
Figure GDA0003604827200000031
the obtained inversion intensity is normalized, the maximum value and the minimum value of the inversion intensity are substituted into a formula, and a final dimensionless inversion intensity coefficient R' is obtained.
Compared with the prior art, the invention has the following beneficial effects:
1. by summarizing the mode of the inversion structural band and analyzing the mechanism of the inversion structural band, formulas of inversion strength and inversion strength coefficients are established, so that the distribution range of the inversion structural band on a plane is rapidly and accurately identified, and the limitation that the prior art cannot accurately identify is broken through.
2. Meanwhile, parameters required in calculation are easy to obtain, universality is high, the plane distribution range of the inversion strip can be accurately identified, and the inversion strength of each area on the inversion structural strip can be further reflected.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a schematic view of the inversion structure with inversion operation according to the present invention;
FIG. 3 is the thickness of the stratums of the sand river block in the research area;
FIG. 4 is a thickness of denudation in the investigation region;
fig. 5 shows a planar distribution range of the inversion structure band calculated by the inversion intensity coefficient.
Detailed Description
For a better understanding of the objects, structure and function of the invention, reference should be made to the following detailed description taken in conjunction with the accompanying drawings.
Referring to fig. 1-5, a method for identifying a planar range of an inversion structural band based on an inversion intensity coefficient comprises the following steps:
s1: establishing inversion patterns of inversion structural bands, as shown in fig. 2, two inversion patterns of inversion structural bands, one being a single-break inversion structural band for single-fault control and the other being a double-break inversion structural band for two-fault control, both types of inversion structural bands undergoing inversion by early settlement and later lifting;
s2: identifying from the cross section of the investigation region whether the inversion structure type corresponds to the pattern established in step S1, identifying the type of inversion structure band in the target region on the cross section as a basis for whether the method is applicable, and it can be seen that the inversion structure band cross section characteristics correspond to the inversion patterns 1 and 2 established in fig. 1, so that the planar range identification of the inversion structure band can be performed by using the inversion intensity factor method;
S3: and (4) establishing a reverse strength formula according to the reverse mode analysis of the step S1, and solving the planar distribution range of the reverse structural belt. Through the mode and mechanism analysis established in the step 1, the stratum thickness on the inversion structural zone is large, the lifting is high, therefore, the stratum bottom interface buries shallowly, and the adjacent depression zone buries deeply although the stratum thickness is large. Wherein, the larger the thickness of the inversion structural belt is, the higher the lifting height is, and the larger the inversion strength is. In addition, although the buried depth of the bottom boundary of the stratum is shallow, the thickness of the stratum is small. Therefore, the thickness of the stratum with the inversion structure and the buried depth of the stratum bottom boundary have certain relations, and the two parameters are relatively easy to obtain. Based on the knowledge, a relational expression of the thickness of the stratum of the inversion structure zone and the buried depth of the bottom boundary of the stratum, namely an inversion strength calculation formula is established;
s4: and (3) calculating the thickness of the stratum in the research area, wherein the thickness of the stratum in the inversion period is delta H, the thickness of the stratum in the target area in the inversion action period is delta H, and the difference between the depth of the top surface and the depth of the bottom surface of the stratum in the inversion period is obtained. As shown in fig. 3, the main period of inversion in the research area is the sand-river street group period, so the depth of the top surface T3 of the sand-river street group stratum is different from the depth of the bottom surface T8 to obtain the thickness Δ H of the sand-river street group stratum;
S5: and recovering the denudation amount of the research area, wherein the thickness of the stratum denudation amount in the reversal period is E, the denudation amount of the research area is subjected to recovery calculation, and the denudation amount recovery adopts a stratum contrast method. As shown in fig. 4, a plurality of two-dimensional cross sections are selected to recover the ablation amount, and the recovered ablation amount is interpolated on a plane to form an interpolated plane ablation amount thickness distribution map;
s6: substituting the stratum thickness obtained in the step S4 and the ablation amount obtained in the step S5 into the inversion strength calculation formula established in the step S2 to obtain the maximum value and the minimum value of the inversion strength;
s7: establishing a reversal intensity coefficient formula, carrying out normalization processing on the obtained reversal intensity, and substituting the maximum value and the minimum value of the reversal intensity obtained in the step S6 into the formula to obtain a final dimensionless reversal intensity coefficient R', so as to identify the range of the reversal structural band on the plane. Fig. 5 is a plane distribution range of the inversion structural band after the inversion strength calculation, and compared with the thickness method of fig. 3, the plane distribution range of the inversion structural band can be more accurately depicted by the present invention. It can be seen from the inversion intensity contour lines in fig. 5 that the inversion intensity is large in the north and south regions of the inversion structure band, and gradually decreases toward the inversion structure band.
Establishing a relational expression of the thickness of the stratum of the inversion structural zone and the buried depth of the bottom boundary of the stratum, namely an inversion strength calculation formula:
Figure GDA0003604827200000051
wherein H Ta The geological significance of the buried depth of the stratum bottom boundary of the inversion structural zone is to reflect the degree of inversion deformation of the inversion structural zone. The denudation amount recovery method adopts a stratum contrast method, the stratum deposition has inheritance and continuity in the same structural layer, according to the characteristic, the deposition thickness of an overlying layer can be estimated according to the product of the ratio of the saved complete adjacent layer thickness (overlying layer thickness/underlying layer thickness) and the underlying layer thickness, and if the estimated value is larger than the overlying residual layer thickness, the exceeding part is the denudation thickness; otherwise, the upper cladding layer is not denuded. According to the above analysis, the adjacent layer thickness ratio is:
Figure GDA0003604827200000052
wherein Ha is the thickness of the adjacent overlying layer that remains intact and Hb is the thickness of the adjacent underlying layer that remains intact. Further calculating the thickness of the stratum denudation as follows:
E=λ×H' b -H' a
wherein E is the thickness of the denudation amount of the stratum, H 'b is the thickness of the overlying layer which is not completely preserved, and H' a is the thickness of the underlying layer which is completely preserved.
The formation thickness Δ H obtained in step S4 and the ablation thickness E obtained in step S5 are substituted into the inversion strength formula created in step S3, and the maximum value Rmax of the inversion strength and the minimum value Rmin of the inversion strength are calculated. Establishing a reversal intensity coefficient formula:
Figure GDA0003604827200000061
The obtained inversion intensity is normalized, the maximum value and the minimum value of the inversion intensity are substituted into a formula, and a final dimensionless inversion intensity coefficient R' is obtained.
It is to be understood that the present invention has been described with reference to certain embodiments, and that various changes in the features and embodiments, or equivalent substitutions may be made therein by those skilled in the art without departing from the spirit and scope of the invention. In addition, many modifications may be made to adapt a particular situation or material to the teachings of the invention without departing from the essential scope thereof. Therefore, it is intended that the invention not be limited to the particular embodiment disclosed, but that the invention will include all embodiments falling within the scope of the appended claims.

Claims (3)

1. A reverse structure belt plane range identification method based on a reverse strength coefficient is characterized in that: the method comprises the following steps:
s1: establishing two inversion structure belt inversion modes, wherein one inversion structure belt is a single-break inversion structure belt controlled by a single fault, and the other inversion structure belt is a double-break inversion structure belt controlled by two faults, and the two types of inversion structure belts both undergo inversion action of early settlement and later lifting;
S2: identifying from the section of the investigation region whether the inversion structure type corresponds to the pattern established in step S1 as a basis for whether the method is applicable;
s3: according to the inversion pattern analysis of the step S1, an inversion strength formula, namely a relation between the thickness of the inversion structure zone stratum and the buried depth of the stratum bottom boundary is established:
Figure FDA0003604827190000011
wherein H Ta The buried depth of the bottom boundary of the stratum of the inversion structural zone has geological significance in reflecting the inversion deformation degree of the inversion structural zone;
s4: the thickness of the stratum in the research area is obtained, and the thickness of the stratum in the reversal period is delta H;
s5: restoring the denudation amount of the research area, wherein the thickness of the stratum denudation amount in the reversal period is E;
s6: substituting the stratum thickness obtained in the step S4 and the ablation amount obtained in the step S5 into the inversion strength calculation formula established in the step S2 to obtain the maximum value and the minimum value of the inversion strength;
s7: and establishing a reversal intensity coefficient formula, normalizing the obtained reversal intensity, and substituting the maximum value and the minimum value of the reversal intensity obtained in the step S6 into the formula to obtain a final dimensionless reversal intensity coefficient R', so that the range of the reversal structural belt is identified on the plane.
2. The inversion configuration band plane range recognition method based on the inversion intensity coefficient as claimed in claim 1, wherein: the formation thickness Δ H obtained in step S4 and the ablation thickness E obtained in step S5 are substituted into the inversion strength formula created in step S3, and the maximum value Rmax of the inversion strength and the minimum value Rmin of the inversion strength are calculated.
3. The inversion-formation-band planar range recognition method based on the inversion intensity coefficient according to claim 2, characterized in that: establishing a reversal intensity coefficient formula:
Figure FDA0003604827190000021
the obtained inversion intensity is normalized, the maximum value and the minimum value of the inversion intensity are substituted into a formula, and a final dimensionless inversion intensity coefficient R' is obtained.
CN202010493517.4A 2020-06-02 2020-06-02 Inversion structure zone plane range identification method based on inversion intensity coefficient Active CN112462417B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010493517.4A CN112462417B (en) 2020-06-02 2020-06-02 Inversion structure zone plane range identification method based on inversion intensity coefficient

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010493517.4A CN112462417B (en) 2020-06-02 2020-06-02 Inversion structure zone plane range identification method based on inversion intensity coefficient

Publications (2)

Publication Number Publication Date
CN112462417A CN112462417A (en) 2021-03-09
CN112462417B true CN112462417B (en) 2022-07-29

Family

ID=74833340

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010493517.4A Active CN112462417B (en) 2020-06-02 2020-06-02 Inversion structure zone plane range identification method based on inversion intensity coefficient

Country Status (1)

Country Link
CN (1) CN112462417B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105317427A (en) * 2014-06-05 2016-02-10 中国石油化工股份有限公司 Representation method for describing old fault top surface structure diagram
CN106324701A (en) * 2016-10-10 2017-01-11 中国海洋石油总公司 Petroliferous basin analysis method based on fracture-uplift linkage
CN107783189A (en) * 2017-02-17 2018-03-09 中国地质大学(北京) A kind of method based on stress types conversion Trend judgement fault strike
CN109903663A (en) * 2019-03-29 2019-06-18 沃肯仕能源科技(北京)有限公司 The Positive Inverted Structures physical simulation experiment device of controllable tomography angle
CN110593861A (en) * 2019-07-29 2019-12-20 中国石油天然气股份有限公司 Quantitative characterization method for structural deformation strength

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101157792B1 (en) * 2011-12-06 2012-06-19 한국지질자원연구원 A method for 3-d geological structure analysis by using structure index
CN105204069B (en) * 2014-06-27 2018-08-17 中国石油化工股份有限公司 A kind of Eroded Thickness restoration methods

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105317427A (en) * 2014-06-05 2016-02-10 中国石油化工股份有限公司 Representation method for describing old fault top surface structure diagram
CN106324701A (en) * 2016-10-10 2017-01-11 中国海洋石油总公司 Petroliferous basin analysis method based on fracture-uplift linkage
CN107783189A (en) * 2017-02-17 2018-03-09 中国地质大学(北京) A kind of method based on stress types conversion Trend judgement fault strike
CN109903663A (en) * 2019-03-29 2019-06-18 沃肯仕能源科技(北京)有限公司 The Positive Inverted Structures physical simulation experiment device of controllable tomography angle
CN110593861A (en) * 2019-07-29 2019-12-20 中国石油天然气股份有限公司 Quantitative characterization method for structural deformation strength

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
Architecture and evolution of the Cenozoic offshore Bohai Bay basin, eastern China;Yixin Yu et al.;《Journal of Asian Earth Sciences》;20200207;第1-15页 *
Cenozoic positive inversion tectonics and its migration in the East China Sea Shelf Basin;GUOHUA ZHANG et al.;《GEOLOGICAL JOURNAL》;20161231;第1-12页 *
Stress field modelling of the Late Oligocene tectonic inversion in the Liaodong Bay Subbasin, Bohai Bay Basin (northern China): Implications for geodynamics and petroleum accumulation;Panpan Hu et al.;《Journal of Geodynamics》;20190314;第32-45页 *
济阳坳陷东营末期构造反转特征分析;鲍倩倩;《中国优秀博硕士学位论文全文数据库(硕士) 基础科学辑》;20100315;第46-56页 *
西湖凹陷反转构造样式与迁移规律;连小翠等;《上海国土资源》;20161231;第37卷(第04期);第83-88页 *
负反转断裂主反转期和反转强度分析;汤良杰等;《石油大学学报(自然科学版)》;19991231;第23卷(第06期);第1-5页 *
辽东湾坳陷新生代构造变形特征及营潍断裂带的表现;康琳等;《石油学报》;20191231;第40卷;第79-90页 *

Also Published As

Publication number Publication date
CN112462417A (en) 2021-03-09

Similar Documents

Publication Publication Date Title
CN110488353B (en) Fault interpretation method based on combination of profile interaction and construction style guidance
CN111458767B (en) Method and system for identifying lithology based on intersection graph method
CN105866832B (en) The recognition methods for the low-grade fault that Deep Tight Sandstone is hidden and device
CN101266299A (en) Method for forecasting oil gas utilizing earthquake data object constructional features
CN103559368B (en) Chemically mechanical polishing analogy method and clearance computational methods thereof
CN107153770A (en) A kind of Forecasting Methodology of Unevenly-pressured shallow tunnel ground settlement deformation
CN106960105A (en) Local scour near pier depth computing method on a kind of stratification alluvial riverbed
CN112462417B (en) Inversion structure zone plane range identification method based on inversion intensity coefficient
CN105426620A (en) Quantitative analysis method and device of main controlling factor of oil layer remaining oil
CN102454400A (en) Method for recognizing carbonate rock crevice cave-shaped reservoir
CN106569272B (en) A kind of seismic attribute fusion method that dimension is risen based on data character space
CN103218750B (en) A kind of geology evaluation method for offshore sedimentary rock oil field
CN110208872A (en) A kind of alluvial fan is every interlayer integrated recognition method
CN109948291A (en) A kind of discontinuous boundary line direction-adaptive recognition methods of sand body
CN109919298A (en) A kind of airfield runway cutting automatic identification and measurement method based on long memory network and Naive Bayes Classifier in short-term
CN113126156B (en) Method and device for extracting high-angle fracture in radon region, storage medium and equipment
CN112711069B (en) Sand point dam prediction method based on plain area curved stream river channel evolution
CN113803047A (en) Method for judging whether fracturing is effective or not based on oil well reservoir physical property data
Huntington et al. The Hurricane fault in southwestern Utah
CN112068201B (en) Exploration method for ancient ridge edge unconformity V-belt
CN112668926A (en) Dam break influence evaluation method
CN111997604B (en) Single-stage river boundary identification method
CN110703325A (en) Hydrocarbon source identification method based on structure difference activity
Shang Research on fast identification model of water-flooded layer in old oilfield--Taking Xingbei area of Daqing Oilfield as an example
CN111665561B (en) Multi-phase superposed river channel tight sandstone oil reservoir dominant reservoir identification method and device

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
GR01 Patent grant
GR01 Patent grant