CN113740930B - Automatic extraction method and device for single well evolution history profile information - Google Patents

Automatic extraction method and device for single well evolution history profile information Download PDF

Info

Publication number
CN113740930B
CN113740930B CN202010482342.7A CN202010482342A CN113740930B CN 113740930 B CN113740930 B CN 113740930B CN 202010482342 A CN202010482342 A CN 202010482342A CN 113740930 B CN113740930 B CN 113740930B
Authority
CN
China
Prior art keywords
data
single well
burial
time
depth
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
CN202010482342.7A
Other languages
Chinese (zh)
Other versions
CN113740930A (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 Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
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 Petroleum and Chemical Corp, Sinopec Exploration and Production Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN202010482342.7A priority Critical patent/CN113740930B/en
Publication of CN113740930A publication Critical patent/CN113740930A/en
Application granted granted Critical
Publication of CN113740930B publication Critical patent/CN113740930B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V20/00Geomodelling in general

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The embodiment of the invention discloses a method and a device for automatically extracting single well evolution history profile information. The method for automatically extracting the profile information of the single well evolution history comprises the following steps: extracting single well data; forming attribute value contour line buried depth data based on the single well data; calculating the burial depth and time of the intersection point of the attribute contour line value and the target layer burial line based on the attribute value contour line burial depth data to obtain burial depth and time data; and extracting the burial depth and time of the set condition from the burial depth and time data to finish the extraction of the single well evolution history profile information. And forming attribute value contour line buried depth data based on the extracted single well data, so as to calculate the buried depth and time of the intersection point of the attribute contour line value and the destination layer buried depth line, and extracting the buried depth and time meeting the set conditions from the buried depth and time of the intersection point. The working efficiency and the precision are effectively improved.

Description

Automatic extraction method and device for single well evolution history profile information
Technical Field
The invention belongs to the technical field of petroleum geology, and particularly relates to an automatic extraction method and device for single well evolution history profile information.
Background
Basin simulation plays an important role in modern oil and gas exploration, petroleum and geological comprehensive research and oil and gas resource evaluation, and one of the simplest and most convenient methods for researching the embedding history, the heat history, the maturity history and the hydrocarbon evolution history of one basin is realized by carrying out one-dimensional simulation on a plurality of single wells in one basin, so that a great deal of research work is carried out by applying the single-well one-dimensional basin simulation by a plurality of students, and a better effect is obtained. Currently, commercial basin simulation systems most widely used in the industry mainly comprise TemisPack of IFP, petroMod of Schlenmez and Basingmod of Platte River Co (PRA) in the United states, three major oil companies in China are intersected and developed to mainly use TSM (medium petrochemical), BASIM (medium petroleum), PRES-MIGS (medium sea oil), and the systems incorporate one-dimensional single well simulation into the basin simulation system without exception, so that the one-dimensional single well basin simulation has great application potential.
The one-dimensional single-well basin simulation (hereinafter referred to as single-well simulation) can display the dynamic evolution process of the oil-gas-containing basin most intuitively and conveniently, and mainly superimposes relevant information such as temperature, ro, hydrocarbon production and the like through a buried history evolution profile, and displays the information to a user through an intuitive drawing, wherein the user needs to further analyze to obtain more detailed information such as different evolution time, temperature, maturation time and the like of a hydrocarbon source stratum after obtaining the single-well evolution history profile. In the existing basin simulation method, when the information is acquired, a user needs to position a target position by using a mouse, then corresponding data is displayed, the user needs a plurality of complicated interactive operations and records to acquire a set of hydrocarbon source rock key evolution time data, the time and the labor are wasted, the accuracy is low, and a data interpretation text needs to be written after the data is acquired. In addition, the data formats in the traditional basin simulation method are numerous, so that the problems of inconvenient reading and low efficiency exist.
Disclosure of Invention
In view of the above, the embodiment of the invention provides a method and a device for automatically extracting single well evolution history profile information, which at least solve the problems of low efficiency and low precision in basin simulation methods in the prior art.
In a first aspect, an embodiment of the present invention provides a method for automatically extracting profile information of single well evolution history, including:
extracting single well data;
forming attribute value contour line buried depth data based on the single well data;
calculating the burial depth and time of the intersection point of the attribute contour line value and the target layer burial line based on the attribute value contour line burial depth data to obtain burial depth and time data;
and extracting the burial depth and time of the set condition from the burial depth and time data to finish the extraction of the single well evolution history profile information.
Optionally, the single well data comprises single well simulation data and single well four history simulation result data;
the single well simulation data comprises a drilling name, a stratum name, stratum bottom age, horizon burial depth and/or hydrocarbon source rock stratum;
the single well four history simulation result data includes burial depths, temperatures, maturity, and/or hydrocarbon formation strength of various layers of the well during the evolution stage.
Optionally, the forming attribute value contour line burial depth data based on the single well data includes:
calculating the maturity point burial depth of each geological evolution stage on the burial history section based on the single well data;
and connecting the maturity points to form a contour line.
Optionally, the forming attribute value contour line burial depth data based on the single well data includes:
setting up a single well simulationSetting k evolution times, corresponding time t 1 ...t k Wherein t is k For the moment at present, it is assumed that the moment t of a geological evolution stage is included in k moments, the stratum involved in evolution is m layers, and the corresponding stratum burial depth is h 1 ...h m Corresponding maturity of Ro 1 ...Ro m The earth surface is 0, the buried depth is positive downwards, let h n+1 =0,Ro n+1 =0, then the end of t Ro value is calculated to correspond to the burial depth Ht as follows:
from t k-1 To t 1 Judging the Ro value at each moment;
if Ro 1 ≤Ro≤Ro m+1 Searching from bottom to top, if Ro is satisfied i Equal to or greater than Ro and Ro i+1 Ro is more than or equal to 1 and i is more than or equal to n, and the following formula is adopted to calculate the buried depth H t ,h i For the depth of burial of the stratum, ro i In order to achieve the degree of maturity,
if Ro > Ro 1 Then search the first Ro from bottom to top j Not equal to Ro 1 Recording j, calculating the burial depth H using the following formula t J is a natural number,
if Ro is greater than the bottom maturity of the lowest stratum of the well, the corresponding burial depth H t1 ...H tk-1 All invalid values are indicated by-1.
Optionally, the calculating the burial depth and time of the intersection point of the attribute contour line value and the destination layer burial depth line based on the attribute value contour line burial depth data, to obtain burial depth and time data, includes:
t Rj =t i +(H ti -H Li )/(H Li+1 -H Li +H ti -H Li )×(t i -t i+1 ),
H Rj =H Li +(H ti -H Li )/(H Li+1 -H Li +H ti -H Li )×(H Li+1 -H Li ),
H ti is buried deep H Li H and H Li+1 Obtaining the i-th and i+1-th evolution stage burial depths corresponding to the L stratum through burial history evolution data; ti is the time of the geological evolution stage, tR j Time of junction, H Rj Is the burial depth of the junction.
Optionally, the extracting the burial depth and time of the set condition from the burial depth and time data includes:
and obtaining s junction burial depths and times for each junction, wherein one junction burial depth and time are geological time and stratum burial depth when the first maturity of the L stratum reaches a corresponding Ro value, extracting the value of the first junction of the contour line and the L stratum according to set conditions, and judging the stratum deposition period corresponding to the moment according to the input stratum age.
Optionally, the method further comprises:
establishing a data normalization template, and performing normalization processing on parameters in single well data, attribute value contour line burial depth data and burial depth and time data based on the data normalization template.
Optionally, in the data normalization template,
p1 is the name of the well;
TY is the hydrocarbon source stratum name;
r1. R4 is the maturity value of a key time node of hydrocarbon source rock evolution;
t1. T4 is the time when the defined hydrocarbon source rock reaches the critical maturity evolution node for the first time;
d1. D4 is the time source rock bottom burial depth;
l1.l4 is the time determined according to T1.t4 compared to the age of the bottom of the drilling formation, if a certain T Dj ≤T i ≤T Dj-1 Li is the name of the stratum of the j-1 layer, ti is the key time determined by the user, T Dj-1 Is T Dj The absolute ages of the j-1 th layer and the stratum layer are input in single well simulation of a user;
r, D is the bottom maturity and burial depth of the hydrocarbon source rock nowadays.
In a second aspect, an embodiment of the present invention further provides an apparatus for automatically extracting profile information of single well evolution history, including:
extraction unit: the method comprises the steps of extracting single well data;
a forming unit: the method comprises the steps of forming attribute value contour line buried depth data based on single well data;
a calculation unit: the method comprises the steps of calculating the burial depth and time of an intersection point of an attribute contour line value and a target layer burial depth line based on the attribute value contour line burial depth data to obtain burial depth and time data;
an information extraction unit: and the method is used for extracting the burial depth and time of the set condition from the burial depth and time data to finish the extraction of the single well evolution history profile information.
Optionally, the method further comprises: specification unit: the method is used for establishing a data normalization template, and parameters in single well data, attribute value contour line burial depth data and burial depth and time data are normalized based on the data normalization template.
According to the method, the attribute value contour line buried depth data are formed based on the extracted single well data, so that the buried depth and time of an intersection point of the attribute contour line value and a destination layer buried depth line are calculated, and the buried depth and time meeting set conditions are extracted from the buried depth and time of the intersection point. The working efficiency and the precision are effectively improved.
And the data is normalized by establishing a data normalization template, so that the working efficiency and the precision are further improved.
Based on single well simulation input data and buried history, heat history and maturity history simulation data, initializing information is formed by reading single well simulation input drilling name, stratum bottom age and horizon buried depth data, data such as stratum buried depth, temperature, maturity and the like of each geological evolution stage of drilling are extracted by reading single well buried history, heat history, maturity history, hydrocarbon generation history and the like simulation data, temperature and maturity contour lines on a single well buried history section are calculated by adopting a linear method, contour line and buried depth line intersection points are obtained according to different application requirements, corresponding time and depth data are extracted according to the contour lines, and data interpretation is formed by combining a normalized description text format, so that the working efficiency and the working accuracy are effectively improved.
Additional features and advantages of the invention will be set forth in the detailed description which follows.
Drawings
The foregoing and other objects, features and advantages of the invention will be apparent from the following more particular descriptions of exemplary embodiments of the invention as illustrated in the accompanying drawings wherein like reference numbers generally represent like parts throughout the exemplary embodiments of the invention.
FIG. 1 shows a schematic diagram of a one-dimensional single well evolution history superimposed maturity profile of one embodiment of the present invention;
FIG. 2 shows a flow chart of a method for automatically extracting single well evolution history profile information according to one embodiment of the invention;
FIG. 3 illustrates a schematic view of a formation and attribute contour junction in a single well history section in accordance with one embodiment of the present invention.
Detailed Description
Preferred embodiments of the present invention will be described in more detail below. While the preferred embodiments of the present invention are described below, it should be understood that the present invention may be embodied in various forms and should not be limited to the embodiments set forth herein.
The invention relates to an automatic extraction method of single well evolution history profile information, which is mainly used for a numerical simulation method of a hydrocarbon basin in petroleum geology research.
An automatic extraction method of single well evolution history profile information comprises the following steps:
extracting single well data;
forming attribute value contour line buried depth data based on the single well data;
calculating the burial depth and time of the intersection point of the attribute contour line value and the target layer burial line based on the attribute value contour line burial depth data to obtain burial depth and time data;
and extracting the burial depth and time of the set condition from the burial depth and time data to finish the extraction of the single well evolution history profile information.
Optionally, the single well data comprises single well simulation data and single well four history simulation result data;
the single well simulation data comprises a drilling name, a stratum name, stratum bottom age, horizon burial depth and/or hydrocarbon source rock stratum;
the single well four history simulation result data includes burial depths, temperatures, maturity, and/or hydrocarbon formation strength of various layers of the well during the evolution stage.
Optionally, the forming attribute value contour line burial depth data based on the single well data includes:
calculating the maturity point burial depth of each geological evolution stage on the burial history section based on the single well data;
and connecting the maturity points to form a contour line.
Optionally, the forming attribute value contour line burial depth data based on the single well data includes:
setting k evolution times for a single well simulation, wherein the corresponding time is t 1 ...t k Wherein t is k For the moment at present, it is assumed that the moment t of a geological evolution stage is included in k moments, the stratum involved in evolution is m layers, and the corresponding stratum burial depth is h 1 ...h m Corresponding maturity of Ro 1 ...Ro m The earth surface is 0, the buried depth is positive downwards, let h n+1 =0,Ro n+1 =0, then calculate the end of t Ro value to the depth H t The method comprises the following steps:
from t k-1 To t 1 Judging the Ro value at each moment;
if Ro 1 ≤Ro≤Ro m+1 Searching from bottom to top, if Ro is satisfied i Equal to or greater than Ro and Ro i+1 Ro is more than or equal to 1 and i is more than or equal to n, and the following formula is adopted to calculate the buried depth H t ,h i For the depth of burial of the stratum, ro i In order to achieve the degree of maturity,
if Ro > Ro 1 Then search the first Ro from bottom to top j Not equal to Ro 1 Recording j, calculating the burial depth H using the following formula t J is a natural number,
if Ro is greater than the bottom maturity of the lowest stratum of the well, the corresponding burial depth H t1 ...H tk-1 All invalid values are indicated by-1.
Optionally, the calculating the burial depth and time of the intersection point of the attribute contour line value and the destination layer burial depth line based on the attribute value contour line burial depth data, to obtain burial depth and time data, includes:
t Rj =t i +(H ti -H Li )/(H Li+1 -H Li +H ti -H Li )×(t i -t i+1 ),
H Rj =H Li +(H ti -H Li )/(H Li+1 -H Li +H ti -H Li )×(H Li+1 -H Li ),
H ti is buried deep H Li H and H Li+1 Obtaining the i-th and i+1-th evolution stage burial depths corresponding to the L stratum through burial history evolution data; ti is the time of the geological evolution stage, tR j Time of junction, H Rj Is the burial depth of the junction.
Optionally, the extracting the burial depth and time of the set condition from the burial depth and time data includes:
and obtaining s junction burial depths and times for each junction, wherein one junction burial depth and time are geological time and stratum burial depth when the first maturity of the L stratum reaches a corresponding Ro value, extracting the value of the first junction of the contour line and the L stratum according to set conditions, and judging the stratum deposition period corresponding to the moment according to the input stratum age.
Optionally, the method further comprises:
establishing a data normalization template, and performing normalization processing on parameters in single well data, attribute value contour line burial depth data and burial depth and time data based on the data normalization template.
Optionally, in the data normalization template,
p1 is the name of the well;
TY is the hydrocarbon source stratum name;
r1. R4 is the maturity value of a key time node of hydrocarbon source rock evolution;
t1. T4 is the time when the defined hydrocarbon source rock reaches the critical maturity evolution node for the first time;
d1. D4 is the time source rock bottom burial depth;
l1.l4 is the time determined according to T1.t4 compared to the age of the bottom of the drilling formation, if a certain T Dj ≤T i ≤T Dj-1 Li is the name of the stratum of the j-1 layer, T i Key time determined for user, T Dj-1 Is T Dj The absolute ages of the j-1 th layer and the stratum layer are input in single well simulation of a user;
r, D is the bottom maturity and burial depth of the hydrocarbon source rock nowadays.
The embodiment of the invention also discloses an automatic extraction device of the single well evolution history profile information, which comprises the following steps:
extraction unit: the method comprises the steps of extracting single well data;
a forming unit: the method comprises the steps of forming attribute value contour line buried depth data based on single well data;
a calculation unit: the method comprises the steps of calculating the burial depth and time of an intersection point of an attribute contour line value and a target layer burial depth line based on the attribute value contour line burial depth data to obtain burial depth and time data;
an information extraction unit: and the method is used for extracting the burial depth and time of the set condition from the burial depth and time data to finish the extraction of the single well evolution history profile information.
Optionally, the method further comprises: specification unit: the method is used for establishing a data normalization template, and parameters in single well data, attribute value contour line burial depth data and burial depth and time data are normalized based on the data normalization template.
Embodiment one:
as shown in fig. 2, the single well simulation not only outputs the burial depth of each stratum at any evolution stage in the ground history, but also can output various attribute information such as the temperature, the maturity, the gas production intensity, etc. of each stratum, and the maturity simulation data description is the most common flow of single well simulation analysis, so that the embodiment uses maturity data information extraction as an example without losing generality, and the description of the implementation is omitted herein, and other attributes are similar.
(1) Single well data extraction: the process relates to two parts of single well simulation input data and single well four history simulation result data. (1) The single well simulation data comprise well drilling name, stratum bottom age, horizon burial depth, hydrocarbon source rock stratum and the like, are basic input data of burial history simulation, are generally stored in an input data part of single well simulation engineering, and can be directly read. (2) The single well four history simulation result data comprises data such as burial depth, temperature, maturity, gas production intensity and the like of each layer of the well in any evolution stage, the data are stored in a simulation result output file and are generally accessed in a plain code or binary format, and the data can be conveniently and directly read after the file storage format is obtained.
(2) Forming attribute value contour (maturity) burial depth data, such as small triangles in fig. 3: calculating the equivalent maturity point burial depth of each geological evolution stage on the burial history section, and connecting the points to form an contour line, wherein the specific steps are as follows:
setting k evolution times for a single well simulation, wherein the corresponding time is t 1 ...t k Wherein t is k For the present moment, assume that the moment t (including the k moments) of a geological evolution stage, the stratum involved in evolution is m layers, and the corresponding stratum burial depth is h 1 ...h m (monotonically decreasing) corresponding maturity of Ro 1 …Ro m (monotonically decreasing), the earth surface is 0, the burial depth is positive, and h is set n+1 =0,Ro n+1 =0, then calculate the end of t Ro value to the depth H t The following steps are adopted.
(1) From t k-1 To t 1 Executing the step (2) once at each moment to obtain H t1 ...H tk-1 When t k For the initial moment of evolution, no stratum is deposited at the moment, and H is led to be tk =H tk-1
(2) If Ro 1 ≤Ro≤Ro m+1 Executing the step (3), otherwise executing the step (4);
(3) searching from bottom to top, if Ro is satisfied i Equal to or greater than Ro and Ro i+1 Ro is more than or equal to 1 and i is more than or equal to n, and H is calculated by adopting a formula 1 t
(4) If Ro > Ro1, then search for the first Ro from bottom to top j Not equal to Ro 1 Recording j, calculating H using equation 2 t
(5) If Ro is greater than the bottom of the lowest formation of the well, H t1 ...H tk-1 All invalid values are indicated by-1, and all intersection points in step (3) are meaningless to be replaced by-1.
(3) Calculating the burial depth and time (small square in fig. 3) of the intersection point of the attribute contour line value and the burial depth line of the target layer: h calculated in step (2) t1 ...H tk Based on the above, further calculate the intersection point of the contour line and the destination stratum, as in FIG. 3, assume that the dotted line in the graph is a certain Ro value contour line R, the intersection points of the contour line and the stratum L are s, where s < k, and the burial depths of these intersection points are set as H R1 ...H Rs Corresponding to the geological time t R1 ...t Rs Setting the jth intersection point of the contour line and the L to be positioned between the i geological evolution stages and the i+1 geological evolution stages, and obtaining the i and i+1 evolution stages corresponding to the L stratum with the burial depth H through burial history evolution data Li H and H Li+1 At this time, k-1 is not less than i is not less than 1. Calculation of t using equation 3 Rj Calculate H using equation 4 Rj
t Rj =t i +(H ti -H Li )/(H Li+1 -H Li +H ti -H Li )×(t i -t i+1 ) (3),
H Rj =H Li +(H ti -H Li )/(H Li+1 -H Li +H ti -H Li )×(H Li+1 -H Li ) (4),
Extracting intersection time and depth values: the s junction depths and times can be obtained by adopting the algorithm of the step 3 for each junction, wherein the junction burial depth and time are geological time and stratum burial depth when the first maturity of the L stratum reaches the corresponding RO value, the value of the first junction of the contour line and the L stratum (oil-producing threshold, a large amount of oil-producing period, gas-producing threshold, gas-producing late period and the like) is generally extracted according to the needs, meanwhile, which set of stratum deposition period the moment is can be judged according to the input stratum age, and the method can be flexibly configured according to the research needs.
(5) Forming a canonical format report: and designing a data description template according to the data description requirement, realizing parameter replacement through programming, and outputting a description report. Wherein the relevant parameters are the data extracted or calculated in the steps (1) - (4).
According to the actual application requirement, the template can be flexibly configured, and the specific style is as follows:
the method comprises the steps that (P1) a well (TY 1) is a hydrocarbon source rock, the hydrocarbon source rock enters a raw oil threshold in a deposition period (about (T1) Ma) of (L1) at the moment, the bottom maturity of the hydrocarbon source rock is (R1), and the formation burial depth is (D1) meters; entering a large amount of oil producing periods in a deposition period (about T2 Ma) of L2, wherein the bottom maturity of the hydrocarbon source rock is R2, and the formation burial depth D2 is m; entering a gas production threshold in a deposition period (about (T3) Ma) of (L3), wherein the bottom maturity of the hydrocarbon source rock is (R3), and the formation burial depth is (D3) m; entering a late gas generating stage in a deposition period (about (T4) Ma) of (L4), wherein the bottom maturity of the hydrocarbon source rock is (R4), and the formation burial depth is (D4) m; the bottom maturity of the hydrocarbon source rock is [ R ] until the today [ TY ], and the burial depth is [ D ] meter.
Wherein: p1 is the name of the well;
TY hydrocarbon source stratum name;
r1. R4 source rock evolution key time node maturity values (user-defined source rock key evolution time maturity values, optionally defined by the user, or more than 4);
t4 is the first time the user-defined source rock reaches the critical maturity evolution node (by step 4, the first junction time value extracted by using the Ro value of the critical time node, (e.g., R1..r4) above, the number is the same as the Ro selected by the user);
d1. D4 is the time source rock bottom burial depth (same as T1. T4 acquisition method);
l1.l4 time determined according to T1.t4 is compared to the age of the bottom of the drilling formation, if T Dj ≤T i ≤T Dj-1 (T i Key time determined for user, T Dj-1 Is T Dj The absolute ages of the j-1 th layer and the stratum layers input in single well simulation of a user are provided, and Li is the name of the j-1 th layer stratum.
R, D is the bottom maturity and burial depth of the hydrocarbon source rock nowadays.
Embodiment two:
for a certain well drilling, practical application is carried out, the time of the simulated horizon is shown in table 1, the time of the hydrocarbon source rock horizons P2 and T3 are shown in table 2, the current buried depths of various strata are shown in table 2, the early developed maturity history simulation is shown in fig. 1, the data of table 3 and table 4 are obtained by applying steps (1) - (4), the information description is carried out by automatically outputting the evolution of the P2 hydrocarbon source rock and the T3 hydrocarbon source rock by applying a description template of step (5), and the concrete steps are as follows:
a P2 hydrocarbon source rock of a certain well enters a raw oil threshold in a T1f deposition period (about 244.35 Ma), the bottom maturity of the hydrocarbon source rock is 0.5%, and the formation burial depth is 1141 m; entering a large amount of oil producing period in a T2 deposition period (about 236.86 Ma), wherein the bottom maturity of the hydrocarbon source rock is 1%, and the stratum is buried at 2766 m; entering a gas generation threshold in a T3 deposition period (about 230.57 Ma), wherein the bottom maturity of the hydrocarbon source rock is 1.3%, and the stratum burial depth is 3062 meters; entering a late gas generation stage at a J12 deposition period (about 195.21 Ma) when the bottom maturity of the hydrocarbon source rock is 2.5%, and the formation burial depth is 4729 meters; the bottom maturity of the P2 source rock is 3.21% until the present day, and the burial depth is 6932 meters.
A certain well T3 hydrocarbon source rock enters a raw oil threshold in a T3x4 deposition period (about 219.18 Ma), the bottom maturity of the hydrocarbon source rock is 0.5%, and the formation burial depth is 829 m; entering a large amount of oil producing period in a J12 sedimentation period (about 187.87 Ma), wherein the bottom maturity of the hydrocarbon source rock is 1%, and the stratum burial depth is 2427 m; entering a gas generation threshold in a J3 deposition period (about 164.65 Ma), wherein the bottom maturity of the hydrocarbon source rock is 1.3%, and the stratum burial depth is 3300 m; entering a late gas generation stage in a-1 deposition period (about-1 Ma), wherein the bottom maturity of the hydrocarbon source rock is-1%, and the stratum is buried at a depth of-1 m; the maturity of the bottom of the T3 source rock is 1.74 percent, and the burial depth is 4405 meters.
The single well evolution history profile information automatic extraction method provided by the invention can be used for configuring a data description template according to the requirements of users aiming at different geological attributes (embedding history, heat history, maturity history and the like), realizes the automatic, rapid and accurate extraction of single well information through computer programming, improves the single well simulation data analysis efficiency and the data acquisition precision, and can be used as important auxiliary for one-dimensional single well simulation work.
Table 1, well formation name and absolute age correspondence table:
table 2, well formation name and bottom land burial depth correspondence table:
numbering device Stratum layer Bottom boundary burial depth/m
11 Q 15
10 K 454
9 J3 1719
8 J12 3100
7 T3x4 3986
6 T3 4405
5 T2 5332
4 T1j 5764
3 T1f 6224
2 P3 6651
1 P2 6932
Table 3, a key moment information table of maturity evolution of a P2 hydrocarbon source rock layer of a certain well:
P1 TY
certain kind P2
R1 L1 D1 T1
0.5 T1f 1141 244.35
R2 L2 D2 T2
1 T2 2767 236.86
R3 L3 D3 T3
1.3 T3 3062 230.57
R4 L4 D4 T4
2.5 J12 4729 195.21
R5 D5
3.21 6932
Table 4, a key time information table of maturity evolution of a hydrocarbon source rock layer of a certain well T3:
P1 TY
certain kind T3
R1 L1 D1 T1
0.5 T3x4 829 219.18
R2 L2 D2 T2
1 J12 2427 187.87
R3 L3 D3 T3
1.3 J3 3300 164.65
R4 L4 D4 T4
2.5 -1 -1 -1
R5 D5
1.74 4405
The foregoing description of embodiments of the invention has been presented for purposes of illustration and description, and is not intended to be exhaustive or limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the various embodiments described.

Claims (8)

1. The method for automatically extracting the profile information of the single well evolution history is characterized by comprising the following steps of:
extracting single well data;
forming attribute value contour line buried depth data based on the single well data;
calculating the burial depth and time of the intersection point of the attribute contour line value and the target layer burial line based on the attribute value contour line burial depth data to obtain burial depth and time data;
extracting the burial depth and time of the set condition from the burial depth and time data, and finishing the extraction of the single well evolution history profile information;
the single well data comprises single well simulation data and single well four-history simulation result data;
the single well simulation data comprises a drilling name, a stratum name, stratum bottom age, horizon burial depth and/or hydrocarbon source rock stratum;
the single well four history simulation result data comprise burial depths, temperatures, maturity and/or gas-oil strength of each layer of the well in an evolution stage;
wherein, the four histories comprise a single well buried history, a heat history, a maturity history and a hydrocarbon generation history;
the forming the attribute value contour line buried depth data based on the single well data comprises the following steps:
calculating the maturity point burial depth of each geological evolution stage on the burial history section based on the single well data;
and connecting the maturity points to form a contour line.
2. The method for automatically extracting profile information of single well evolution history according to claim 1, wherein the forming attribute value contour line buried depth data based on the single well data comprises:
setting k evolution times for a single well simulation, wherein the corresponding time is t 1 …t k Wherein t is k For the moment at present, it is assumed that the moment t of a geological evolution stage is included in k moments, the stratum involved in evolution is m layers, and the corresponding stratum burial depth is h 1 …h m Corresponding maturity of Ro 1 …Ro m The earth surface is 0, the buried depth is positive downwards, let h n+1 =0,Ro n+1 =0, then calculate the end of t Ro value to the depth H t The method comprises the following steps:
from t k-1 To t 1 Judging the Ro value at each moment;
if Ro 1 ≤Ro≤Ro m+1 Searching from bottom to top, if Ro is satisfied i Equal to or greater than Ro and Ro i+1 Ro is more than or equal to 1 and i is more than or equal to n, and the following formula is adopted to calculate the buried depth H t ,h i Is the depth of burial of stratum, rv i In order to achieve the degree of maturity,
if Ro>Ro 1 Then search the first Ro from bottom to top j Not equal to Ro 1 Recording j, calculating the burial depth H using the following formula t J is a natural number,
if Ro is greater than the bottom maturity of the lowest stratum of the well, the corresponding burial depth H t1 …H tk-1 All invalid values are indicated by-1.
3. The method for automatically extracting profile information of single well evolution history according to claim 2, wherein the calculating the burial depth and time of the intersection of the attribute contour value and the destination layer burial depth line based on the attribute value contour line burial depth data, to obtain the burial depth and time data, comprises:
t Rj =t i +(H ti -H Li )/(H Li+1 -H Li +H ti -H Li )×(t i -t i+1 ),
H Rj =H Li +(H ti -H Li )/(H Li+1 -H Li +H ti -H Li )×(H Li+1 -H Li ),
H ti is buried deep H Li H and H Li+1 Obtaining the i-th and i+1-th evolution stage burial depths corresponding to the L stratum through burial history evolution data; t is t i For the time of geological evolution stage, t Rj Time of junction, H Rj Is the burial depth of the junction.
4. The method for automatically extracting profile information of single well evolution history according to claim 3, wherein the extracting the burial depth and time of the set condition from the burial depth and time data comprises:
and obtaining s junction burial depths and times for each junction, wherein one junction burial depth and time are geological time and stratum burial depth when the first maturity of the L stratum reaches a corresponding RO value, extracting the value of the first junction of the contour line and the L stratum according to set conditions, and judging the stratum deposition period corresponding to the moment according to the input stratum age.
5. The method for automatically extracting profile information of single well evolution history according to claim 1, further comprising:
establishing a data normalization template, and performing normalization processing on parameters in single well data, attribute value contour line burial depth data and burial depth and time data based on the data normalization template.
6. The method for automatically extracting single well evolution history profile information according to claim 5, wherein in the data normalization template,
p1 is the name of the well;
TY is the hydrocarbon source stratum name;
r1 … R4 is a hydrocarbon source rock evolution key time node maturity value;
t1 … T4 is defined time for the hydrocarbon source rock to reach the critical maturity evolution node for the first time;
d1 … D4 is the time source rock bottom burial depth;
l1 … L4 is the time determined from T1 … T4 compared to the age of the bottom of the formation, if a certain T Dj ≤T i ≤T Dj-1 L is then i For the j-1 th layer stratum name, T i Key time determined for user, T Dj-1 Is T Dj The absolute ages of the j-1 th layer and the stratum layer are input in single well simulation of a user;
r, D is the bottom maturity and burial depth of the hydrocarbon source rock nowadays.
7. An automatic extraction device of single well evolution history profile information, which is characterized by comprising:
extraction unit: the method comprises the steps of extracting single well data;
a forming unit: the method comprises the steps of forming attribute value contour line buried depth data based on single well data;
a calculation unit: the method comprises the steps of calculating the burial depth and time of an intersection point of an attribute contour line value and a target layer burial depth line based on the attribute value contour line burial depth data to obtain burial depth and time data;
an information extraction unit: the method comprises the steps of extracting the burial depth and time of a set condition from the burial depth and time data, and completing the extraction of single well evolution history profile information;
the single well data comprises single well simulation data and single well four-history simulation result data;
the single well simulation data comprises a drilling name, a stratum name, stratum bottom age, horizon burial depth and/or hydrocarbon source rock stratum;
the single well four history simulation result data comprise burial depths, temperatures, maturity and/or gas-oil strength of each layer of the well in an evolution stage;
wherein, the four histories comprise a single well buried history, a heat history, a maturity history and a hydrocarbon generation history;
the forming the attribute value contour line buried depth data based on the single well data comprises the following steps:
calculating the maturity point burial depth of each geological evolution stage on the burial history section based on the single well data;
and connecting the maturity points to form a contour line.
8. The method for automatically extracting profile information of single well evolution history according to claim 7, further comprising:
specification unit: the method is used for establishing a data normalization template, and parameters in single well data, attribute value contour line burial depth data and burial depth and time data are normalized based on the data normalization template.
CN202010482342.7A 2020-05-29 2020-05-29 Automatic extraction method and device for single well evolution history profile information Active CN113740930B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010482342.7A CN113740930B (en) 2020-05-29 2020-05-29 Automatic extraction method and device for single well evolution history profile information

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010482342.7A CN113740930B (en) 2020-05-29 2020-05-29 Automatic extraction method and device for single well evolution history profile information

Publications (2)

Publication Number Publication Date
CN113740930A CN113740930A (en) 2021-12-03
CN113740930B true CN113740930B (en) 2024-02-27

Family

ID=78727982

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010482342.7A Active CN113740930B (en) 2020-05-29 2020-05-29 Automatic extraction method and device for single well evolution history profile information

Country Status (1)

Country Link
CN (1) CN113740930B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107861152A (en) * 2017-10-01 2018-03-30 中国石油大港油田勘探开发研究院 A kind of Mesozoic Strata fine correlation method
US10190998B1 (en) * 2018-08-29 2019-01-29 Research Institute Of Petroleum Exploration & Development, Dagang Oil Field Of Cnpc Method and device for evaluating and predicting a shale oil enrichment areas of fault lacustrine basins
CN110517794A (en) * 2019-08-23 2019-11-29 长安大学 A method of that establishes shale gas reservoir buries-thermal evolution history figure
CN110555221A (en) * 2018-06-01 2019-12-10 中国石油化工股份有限公司 Method and device for calculating lifting amplitude of regional stratum

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9411071B2 (en) * 2012-08-31 2016-08-09 Exxonmobil Upstream Research Company Method of estimating rock mechanical properties

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107861152A (en) * 2017-10-01 2018-03-30 中国石油大港油田勘探开发研究院 A kind of Mesozoic Strata fine correlation method
CN110555221A (en) * 2018-06-01 2019-12-10 中国石油化工股份有限公司 Method and device for calculating lifting amplitude of regional stratum
US10190998B1 (en) * 2018-08-29 2019-01-29 Research Institute Of Petroleum Exploration & Development, Dagang Oil Field Of Cnpc Method and device for evaluating and predicting a shale oil enrichment areas of fault lacustrine basins
CN110517794A (en) * 2019-08-23 2019-11-29 长安大学 A method of that establishes shale gas reservoir buries-thermal evolution history figure

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Per Avseth etal.combining burial history and rock-physics modeling to constrain AVO analysis during exploration.2016,第35卷(第6期),474-560. *
贵州织纳煤田水公河向斜上二叠统8煤层三史模拟;鲍园;韦重韬;王超勇;曹佳;;煤田地质与勘探(第06期);全文 *
鄂尔多斯盆地塔巴庙地区上古生界源岩热演化特征;戴世立;杨智;李弘;杨姣;;海洋石油(第01期);全文 *
长岭凹陷断陷期生烃阶段性与成藏圈闭的选择性;徐春华;肖永军;邱连贵;;新疆石油地质(第03期);全文 *
饶阳凹陷北部洼陷区古近系油气运聚机理与成藏模式;韩杨等;天然气地球科学;第31卷(第3期);358-369 *

Also Published As

Publication number Publication date
CN113740930A (en) 2021-12-03

Similar Documents

Publication Publication Date Title
CN111950051B (en) BIM-based three-dimensional geological modeling and geological model-based construction application method
CN104380144B (en) Three-dimensional multimode formula core and Geologic modeling for optimal oil field development
CN101221250A (en) Three-dimensional stratum characteristic decryption method based on geologic information acquisition data
CN111950046B (en) Drilling data model construction method based on BIM
Gherabati et al. Assessment of hydrocarbon in place and recovery factors in the Eagle Ford Shale play
Drovosekova et al. Peculiarities of modelling hydro-lithospheric processes in the region of Kavkazskiye Mineralnye Vody (caucasus mineral springs)
CN108986213B (en) Three-dimensional stratum modeling method based on stacking technology
CN106547941A (en) A kind of method that BIM Engineering Exploration Information Systems are built based on revit figure models
CN107180452A (en) The stratified geological objects modeling method of the multi-source data coupling on the stratum containing standard
CN110197013B (en) Riverbed bedrock surface modeling method based on Morphing
Taipova et al. Verifying reserves opportunities with multi-well pressure pulse-code testing
CN113740930B (en) Automatic extraction method and device for single well evolution history profile information
CN111472763B (en) Stratum thickness prediction method and device
O'Sullivan et al. An integrated, mesh-independent geothermal modelling framework
Meyer et al. Pecos Valley aquifer, west Texas: structure and brackish groundwater
CN108074282B (en) Thickness correction method of physical property parameters based on three-dimensional geological modeling software
CN112991531B (en) Dynamic construction method of refined three-dimensional hydrogeologic model
Korus et al. Integrating AEM and borehole data for regional hydrogeologic synthesis: Tools and examples from Nebraska, USA
Robertson et al. Optimization of unconventional well-pad area using reservoir simulation and intelligent sequential sampling
Rafi et al. Preliminary prioritization on steam flood injection in oil field using random forest regression method
Detay et al. An expert system for groundwater exploration and exploitation
CN117251739B (en) Construction method of sample set for sandstone uranium resource quantitative prediction
CN112037329B (en) Irregular curved surface construction method
Cannia et al. Cooperative hydrology study COHYST hydrostratigraphic units and Aquifer characterization report
Wood et al. Model creation based on digital borehole records and interpreted geological cross‐sections

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