CN113740930A - Method and device for automatically extracting profile information of single well evolution history - Google Patents
Method and device for automatically extracting profile information of single well evolution history Download PDFInfo
- Publication number
- CN113740930A CN113740930A CN202010482342.7A CN202010482342A CN113740930A CN 113740930 A CN113740930 A CN 113740930A CN 202010482342 A CN202010482342 A CN 202010482342A CN 113740930 A CN113740930 A CN 113740930A
- Authority
- CN
- China
- Prior art keywords
- data
- time
- well
- depth
- burial
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000009933 burial Methods 0.000 claims abstract description 124
- 238000000605 extraction Methods 0.000 claims abstract description 29
- IYLGZMTXKJYONK-ACLXAEORSA-N (12s,15r)-15-hydroxy-11,16-dioxo-15,20-dihydrosenecionan-12-yl acetate Chemical compound O1C(=O)[C@](CC)(O)C[C@@H](C)[C@](C)(OC(C)=O)C(=O)OCC2=CCN3[C@H]2[C@H]1CC3 IYLGZMTXKJYONK-ACLXAEORSA-N 0.000 claims abstract description 18
- IYLGZMTXKJYONK-UHFFFAOYSA-N ruwenine Natural products O1C(=O)C(CC)(O)CC(C)C(C)(OC(C)=O)C(=O)OCC2=CCN3C2C1CC3 IYLGZMTXKJYONK-UHFFFAOYSA-N 0.000 claims abstract description 18
- 238000004088 simulation Methods 0.000 claims description 48
- 239000011435 rock Substances 0.000 claims description 43
- 239000004215 Carbon black (E152) Substances 0.000 claims description 28
- 229930195733 hydrocarbon Natural products 0.000 claims description 28
- 150000002430 hydrocarbons Chemical class 0.000 claims description 28
- 238000010606 normalization Methods 0.000 claims description 22
- 230000015572 biosynthetic process Effects 0.000 claims description 16
- 230000008021 deposition Effects 0.000 claims description 16
- 238000012545 processing Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 3
- 239000003921 oil Substances 0.000 description 14
- 238000004519 manufacturing process Methods 0.000 description 9
- 239000003208 petroleum Substances 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000005553 drilling Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 239000010779 crude oil Substances 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 206010027336 Menstruation delayed Diseases 0.000 description 1
- 208000009989 Posterior Leukoencephalopathy Syndrome Diseases 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000013075 data extraction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000010198 maturation time Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling 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 profile information of single well evolution history. 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 isoline burial depth data based on the single well data; calculating the burial depth and time of an intersection point of the attribute contour value and the buried depth of the target layer based on the attribute value contour buried depth data to obtain buried depth and time data; and extracting the burial depth and time of set conditions from the burial depth and time data to finish the extraction of the profile information of the single well evolution history. The buried depth and time of an intersection point of the attribute contour value and the buried depth of the target layer are calculated by forming attribute value contour buried depth data based on the extracted single well data, and the buried depth and time meeting set conditions are extracted from the buried depth and time of the intersection point. Effectively improve work efficiency and precision.
Description
Technical Field
The invention belongs to the technical field of petroleum geology, and particularly relates to a method and a device for automatically extracting profile information of a single well evolution history.
Background
Basin simulation plays an important role in modern oil and gas exploration, petroleum geological comprehensive research and oil and gas resource evaluation, and one-dimensional simulation of a plurality of single wells in a basin is one of the simplest and most convenient methods for researching the basin burial history, thermal history, maturity history and hydrocarbon evolution history, and a plurality of scholars use single-well one-dimensional basin simulation to carry out a large amount of research work and obtain better effects. At present, the commercial basin simulation systems most widely applied in the industry mainly include temisipack of IFP, PetroMod of schlumberger and basemod of plate River company (PRA) in the united states, and TSM (mesopetro chemical), BASIMS (medium petroleum oil), PRES-MIGS (medium sea oil) mainly applied are developed by three oil companies in China in a crossed manner.
The dynamic evolution process of the oil-gas-containing basin can be displayed most intuitively and conveniently by one-dimensional single-well basin simulation (hereinafter referred to as single-well simulation), relevant information such as temperature, Ro and hydrocarbon generation amount is mainly displayed to a user through a visual graph, and the user needs to further analyze after obtaining the single-well evolution history profile to obtain more detailed information such as different evolution time, temperature and maturation time of a hydrocarbon source rock stratum. According to the existing basin simulation method, when the information is obtained, a user needs to use a mouse to position a target position and then display corresponding data, the user needs to obtain a set of key evolution time data of the hydrocarbon source rocks through a plurality of complex interactive operations and records, time and labor are wasted, the accuracy is inaccurate, and data interpretation texts need to be compiled after the data are obtained. And the prior basin simulation method has numerous data formats, thereby having the problems of inconvenient reading and low efficiency.
Disclosure of Invention
In view of this, the embodiment of the invention provides a method and a device for automatically extracting profile information of a single well evolution history, which at least solve the problem of low efficiency and low precision in a basin simulation method in the prior art.
In a first aspect, an embodiment of the present invention provides an automatic extraction method for single-well evolution history profile information, including:
extracting single well data;
forming attribute value isoline burial depth data based on the single well data;
calculating the burial depth and time of an intersection point of the attribute contour value and the buried depth of the target layer based on the attribute value contour buried depth data to obtain buried depth and time data;
and extracting the burial depth and time of set conditions from the burial depth and time data to finish the extraction of the profile information of the single well evolution history.
Optionally, the single-well data includes single-well simulation data and single-well four-history simulation result data;
the single well simulation data comprises a well name, a stratum bottom age, a horizon burial depth and/or a hydrocarbon source rock stratum;
the single well four-history simulation result data comprises the burial depth, the temperature, the maturity and/or the raw oil gas strength of each drill well horizon in the evolution stage.
Optionally, the forming attribute value contour buried depth data based on the single well data includes:
calculating the burial depth of maturity points of each geological evolution stage on a burial history profile based on the single-well data;
and connecting the maturity points to form a contour line.
Optionally, the forming attribute value contour buried depth data based on the single well data includes:
setting k evolution times for a single well simulation at a certain time, wherein the corresponding time is t1...tkWherein t iskFor the present moment, suppose a certain geological evolution stage t moment, the t moment is contained in k moments, the stratum participating in the evolution is m layers, and the corresponding stratum burial depth is h1...hmCorresponding maturity of Ro1...RomThe earth surface is 0, the buried depth is positive downwards, and h is setn+1＝0，Ron+1If the value is 0, calculating the corresponding burial depth Ht of the Ro value at the end of the time t as follows:
from tk-1To t1Judging the Ro value at each moment;
if Ro1≤Ro≤Rom+1Searching from bottom to top, if Ro is satisfiediRo and not less than Roi+1Ro is not more than 1, i is not less than 1 and n is not less than n, the buried depth H is calculated by adopting the following formulat，hiFor formation burial depth, RoiIn order to be of a degree of maturity,
if Ro > Ro1Then search for the first Ro from bottom to topjNot equal to Ro1J, the buried depth H is calculated by the following formulatJ is a natural number,
if Ro is greater than the bottom maturity of the lowest stratum of the well, the corresponding buried depth Ht1...Htk-1All invalid values are represented by-1.
Optionally, the calculating the burial depth and the time of the intersection point of the attribute contour value and the target layer burial depth line based on the attribute value contour burial depth data to obtain the burial depth and time data includes:
tRj＝ti+(Hti-HLi)/(HLi+1-HLi+Hti-HLi)×(ti-ti+1)，
HRj＝HLi+(Hti-HLi)/(HLi+1-HLi+Hti-HLi)×(HLi+1-HLi)，
Htito be buried deep, HLiAnd HLi+1Acquiring the burial depths of the ith and i +1 evolution stages corresponding to the L stratum through the evolution data of the burial history; ti is the geological evolution stage time, tRjTime of intersection, HRjIs the burial depth of the junction.
Optionally, the extracting the burial depth and the time of the set condition from the burial depth and time data includes:
and obtaining s intersection point burial depths and time for each intersection point, wherein one of the intersection point burial depths and time is geological time and stratum burial depth when the first maturity of the L stratum reaches the corresponding Ro value, extracting a first intersection point value of the isoline and the L stratum according to set conditions, and judging a stratum deposition period corresponding to the moment according to the input stratum age.
Optionally, the method further includes:
and establishing a data normalization template, and performing normalization processing on parameters in the single-well data, the attribute value isoline burial depth data and the burial depth and time data based on the data normalization template.
Optionally, in the data normalization template,
p1 is the well name;
TY is the name of the hydrocarbon source formation;
r1.. R4 is a maturity value of a key time node of the evolution of the source rock;
t1.. T4 is the time for a defined source rock to first reach a key maturity evolutionary node;
d4 is the time source rock bottom burial depth;
l1.. L4 is the time determined according to T1.. T4 compared with the age of the bottom of the drilled formation, if TDj≤Ti≤TDj-1Li is the name of the stratum of the j-1 th layer, Ti is the key time determined by the user, and TDj-1Is TDjLayer j-1 and stratum insulation input in single well simulation by userFor the age;
r, D is the bottom maturity and burial depth of the source rock today.
In a second aspect, an embodiment of the present invention further provides an automatic extraction apparatus for single-well evolution history profile information, including:
an extraction unit: for extracting single well data;
forming a unit: for forming attribute value isopleth burial depth data based on the single well data;
a calculation unit: the buried depth and time of the intersection point of the attribute contour value and the buried depth of the target layer are calculated based on the attribute value contour buried depth data, and buried depth and time data are obtained;
an information extraction unit: and the method is used for extracting the burial depth and time of set conditions from the burial depth and time data to finish the extraction of the profile information of the single well evolution history.
Optionally, the method further includes: a specification unit: the method is used for establishing a data normalization template, and carrying out normalization processing on parameters in single well data, attribute value isoline burial depth data and burial depth and time data based on the data normalization template.
According to the method, attribute value isoline burial depth data are formed on the basis of single well data extraction, so that the burial depth and time of an intersection point of the attribute isoline value and a target layer burial depth line are calculated, and the burial depth and time meeting set conditions are extracted from the burial depth and time of the intersection point. Effectively improve work efficiency and precision.
And the data is normalized by establishing a data normalization template, so that the working efficiency and the precision are further improved.
The method comprises the steps of taking single well simulation input data and simulation data of a buried history, a thermal history and a maturity history as bases, forming initialization information by reading the drilling name, the stratum bottom age and the horizon buried depth data which are input by single well simulation, extracting the data of the stratum buried depth, the temperature, the maturity and the like of each geological evolution stage of the drilling well by reading the simulation data of the single well buried history, the thermal history, the maturity history, the hydrocarbon generation history and the like, calculating the temperature and the maturity contour line on the section of the single well buried history by adopting a linear method, obtaining the intersection point of the contour line and the buried depth line according to different application requirements, extracting corresponding time and depth data according to the contour line, combining a standardized description text format to form a data explanation, and effectively improving the working efficiency and precision.
Additional features and advantages of the invention will be set forth in the detailed description which follows.
Drawings
The above and other objects, features and advantages of the present invention will become more apparent by describing in more detail exemplary embodiments thereof with reference to the attached drawings, in which like reference numerals generally represent like parts throughout.
FIG. 1 is a schematic diagram showing a one-dimensional single well evolution history superimposed maturity profile of one embodiment of the present invention;
FIG. 2 illustrates a flow chart of a method for automatic extraction of single well evolution history profile information according to an embodiment of the invention;
FIG. 3 is a schematic diagram illustrating the junction of a stratigraphic layer and a property contour in a single well evolution history profile according to an embodiment of the present invention.
Detailed Description
Preferred embodiments of the present invention will be described in more detail below. While the following describes preferred embodiments of the present invention, it should be understood that the present invention may be embodied in various forms and should not be limited by the embodiments set forth herein.
The invention relates to an automatic extraction method of profile information of single well evolution history, which is mainly used for a numerical simulation method of a hydrocarbon-containing basin in petroleum geological research.
An automatic extraction method for single-well evolution history profile information comprises the following steps:
extracting single well data;
forming attribute value isoline burial depth data based on the single well data;
calculating the burial depth and time of an intersection point of the attribute contour value and the buried depth of the target layer based on the attribute value contour buried depth data to obtain buried depth and time data;
and extracting the burial depth and time of set conditions from the burial depth and time data to finish the extraction of the profile information of the single well evolution history.
Optionally, the single-well data includes single-well simulation data and single-well four-history simulation result data;
the single well simulation data comprises a well name, a stratum bottom age, a horizon burial depth and/or a hydrocarbon source rock stratum;
the single well four-history simulation result data comprises the burial depth, the temperature, the maturity and/or the raw oil gas strength of each drill well horizon in the evolution stage.
Optionally, the forming attribute value contour buried depth data based on the single well data includes:
calculating the burial depth of maturity points of each geological evolution stage on a burial history profile based on the single-well data;
and connecting the maturity points to form a contour line.
Optionally, the forming attribute value contour buried depth data based on the single well data includes:
setting k evolution times for a single well simulation at a certain time, wherein the corresponding time is t1...tkWherein t iskFor the present moment, suppose a certain geological evolution stage t moment, the t moment is contained in k moments, the stratum participating in the evolution is m layers, and the corresponding stratum burial depth is h1...hmCorresponding maturity of Ro1...RomThe earth surface is 0, the buried depth is positive downwards, and h is setn+1＝0，Ron+1If equal to 0, calculating the Ro value corresponding to the buried depth H at the end of t timetThe method comprises the following specific steps:
from tk-1To t1Judging the Ro value at each moment;
if Ro1≤Ro≤Rom+1Searching from bottom to top, if Ro is satisfiediRo and not less than Roi+1Ro is not more than 1, i is not less than 1 and n is not less than n, the buried depth H is calculated by adopting the following formulat，hiFor formation burial depth, RoiIn order to be of a degree of maturity,
if Ro > Ro1Then search for the first Ro from bottom to topjNot equal to Ro1J, the buried depth H is calculated by the following formulatJ is a natural number,
if Ro is greater than the bottom maturity of the lowest stratum of the well, the corresponding buried depth Ht1...Htk-1All invalid values are represented by-1.
Optionally, the calculating the burial depth and the time of the intersection point of the attribute contour value and the target layer burial depth line based on the attribute value contour burial depth data to obtain the burial depth and time data includes:
tRj＝ti+(Hti-HLi)/(HLi+1-HLi+Hti-HLi)×(ti-ti+1)，
HRj＝HLi+(Hti-HLi)/(HLi+1-HLi+Hti-HLi)×(HLi+1-HLi)，
Htito be buried deep, HLiAnd HLi+1Acquiring the burial depths of the ith and i +1 evolution stages corresponding to the L stratum through the evolution data of the burial history; ti is the geological evolution stage time, tRjTime of intersection, HRjIs the burial depth of the junction.
Optionally, the extracting the burial depth and the time of the set condition from the burial depth and time data includes:
and obtaining s intersection point burial depths and time for each intersection point, wherein one of the intersection point burial depths and time is geological time and stratum burial depth when the first maturity of the L stratum reaches the corresponding Ro value, extracting a first intersection point value of the isoline and the L stratum according to set conditions, and judging a stratum deposition period corresponding to the moment according to the input stratum age.
Optionally, the method further includes:
and establishing a data normalization template, and performing normalization processing on parameters in the single-well data, the attribute value isoline burial depth data and the burial depth and time data based on the data normalization template.
Optionally, in the data normalization template,
p1 is the well name;
TY is the name of the hydrocarbon source formation;
r1.. R4 is a maturity value of a key time node of the evolution of the source rock;
t1.. T4 is the time for a defined source rock to first reach a key maturity evolutionary node;
d4 is the time source rock bottom burial depth;
l1.. L4 is the time determined according to T1.. T4 compared with the age of the bottom of the drilled formation, if TDj≤Ti≤TDj-1Li is the name of the stratum of the j-1 th layer, TiKey time, T, determined for the userDj-1Is TDjThe absolute age of the layer j-1 and the layer stratum input in the single well simulation of the user;
r, D is the bottom maturity and burial depth of the source rock today.
The embodiment of the invention also discloses an automatic extraction device for the profile information of the single well evolution history, which comprises the following steps:
an extraction unit: for extracting single well data;
forming a unit: for forming attribute value isopleth burial depth data based on the single well data;
a calculation unit: the buried depth and time of the intersection point of the attribute contour value and the buried depth of the target layer are calculated based on the attribute value contour buried depth data, and buried depth and time data are obtained;
an information extraction unit: and the method is used for extracting the burial depth and time of set conditions from the burial depth and time data to finish the extraction of the profile information of the single well evolution history.
Optionally, the method further includes: a specification unit: the method is used for establishing a data normalization template, and carrying out normalization processing on parameters in single well data, attribute value isoline burial depth data and burial depth and time data based on the data normalization template.
The first embodiment is as follows:
as shown in fig. 2, because the single-well simulation not only outputs the burial depth of each stratum at any evolution stage in the history, but also outputs a plurality of attribute information such as the temperature, maturity, crude oil gas strength, etc. of each stratum, the maturity simulation data indicates the most common flow for single-well simulation analysis, and therefore, the embodiment takes the extraction of the maturity data information as an example to explain the implementation manner without loss of generality, and other attributes are similar, and are not described herein.
(1) Extracting single-well data: the process relates to two parts of single-well simulation input data and single-well four-history simulation result data. The single well simulation data comprises well drilling names, stratum bottom ages, horizon burial depths, hydrocarbon source strata and the like, and the data are basic input data of the burial history simulation, are generally stored in an input data part of a single well simulation project and can be directly read. Secondly, the single-well four-history simulation result data comprises data of burial depth, temperature, maturity, crude oil gas strength and the like of each drill well layer 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, the data can be conveniently and directly read after the file storage format is obtained, and the maturity data is taken as an example in the embodiment.
(2) Attribute value contour (maturity) burial depth data are formed, as shown by the small triangles in fig. 3: calculating the burial depths of maturity points such as geological evolution stages on a burial history profile, and connecting the points to form an isoline, wherein the method comprises the following specific steps:
setting k evolution times for a single well simulation at a certain time, wherein the corresponding time is t1...tkWherein t iskFor the present moment, it is assumed that at time t (including the k times) of a certain geological evolution stage, the strata involved in the evolution are m layers, and the corresponding stratum burial depth is h1...hm(monotonically decreasing) corresponds to a maturity of Ro1…Rom(monotonically decreasing), surface 0, buried depth positive, set hn+1＝0，Ron+1If equal to 0, calculating the Ro value corresponding to the buried depth H at the end of t timetThe following procedure was used.
From tk-1To t1Step two is executed at each moment to obtain Ht1...Htk-1When t iskFor the initial time of evolution, no stratum deposition exists at this time, and then, H is orderedtk＝Htk-1。
② if Ro1≤Ro≤Rom+1If not, executing the step IV;
searching from bottom to top, if Ro is satisfiediRo and not less than Roi+1Ro is not more than 1, i is not less than 1 and n is not less than n, then formula 1 is adopted to calculate Ht；
Fourthly, if Ro is more than Ro1, searching the first Ro from bottom to topjNot equal to Ro1J, using equation 2 to calculate Ht
Fifthly, if Ro is larger than the bottom maturity of the lowest stratum of the well, Ht1...Htk-1All invalid values are represented by-1, and the junction points in the step (3) are meaningless replaced by-1.
(3) Calculating the burial depth and time (small square in fig. 3) of the intersection point of the attribute contour value and the buried depth line of the target layer: h calculated in step (2)t1...HtkBased on the above, further calculating the intersection point of the contour line and the target stratum, as the dotted line in the figure is a Ro value contour line R, the intersection points of the contour line and the stratum L are s, wherein s is less than k, and the burial depth of the intersection points is set as HR1...HRsCorresponding to geological time tR1...tRsSetting the j-th intersection point of the isoline and L between i geological evolution stages and i +1 geological evolution stages, and obtaining the buried depth of the i-th and i +1 evolution stages corresponding to the L stratum as H through the buried history evolution dataLiAnd HLi+1At this time, k-1 is not less than i not less than 1. Calculating t using equation 3RjCalculating H using equation 4Rj。
tRj＝ti+(Hti-HLi)/(HLi+1-HLi+Hti-HLi)×(ti-ti+1) (3)，
HRj＝HLi+(Hti-HLi)/(HLi+1-HLi+Hti-HLi)×(HLi+1-HLi) (4)，
Extracting intersection point time and depth values: the depth and time of the s junction points can be obtained by adopting the algorithm of the step 3 for each junction point, wherein the burial depth and time of one junction point are the geological time when the first maturity of the L stratum reaches the corresponding RO value and the burial depth of the stratum, the first junction point value of the isoline and the L stratum is generally extracted according to the requirement (an oil production threshold, a large oil production period, a gas production threshold, a gas production late period and the like), and meanwhile, the stratum deposition period of the moment is judged according to the input stratum age, and the depth and time of the junction points can be flexibly configured according to the research requirement.
(5) Forming a standard format report: 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 data extracted or calculated in the steps (1) to (4).
According to the actual application needs, the template can be flexibly configured, and the specific styles are as follows:
a [ P1 ] well [ TY1 ] hydrocarbon source rock, the maturity of the bottom of the hydrocarbon source rock is [ R1 ] when the [ L1 ] deposition period (about [ T1 ] Ma) enters an oil production threshold, and the stratum is buried deep [ D1 ] m; entering a large-amount oil production period in a deposition period (about [ T2 ] Ma) of [ L2 ], wherein the maturity of the bottom of a hydrocarbon source rock is [ R2 ], and the stratum is buried deep by [ D2 ] m; entering a gas generation threshold in a deposition period (about [ T3 ] Ma) of [ L3 ], wherein the maturity of the bottom of a hydrocarbon source rock is [ R3 ], and the stratum burial depth is [ D3 ] m; entering a late gas generation stage in a deposition stage (about [ T4 ] Ma) of [ L4 ], wherein the maturity of the bottom of a hydrocarbon source rock is [ R4 ], and the stratum burial depth is [ D4 ] m; until now, the bottom maturity of the hydrocarbon source rock is [ R ], and the burial depth is [ D ].
In the formula: p1 is the well name;
TY hydrocarbon source formation name;
r4 hydrocarbon source rock evolution key time node maturity values (user-defined hydrocarbon source rock key evolution time maturity values, defined according to user needs, may be more than 4);
t1.. T4 is the time for the user-defined source rock to reach the key maturity evolution node for the first time (through step 4, the Ro value of the key time node is used, and the first intersection point time value extracted (as the above-mentioned R1.. R4) is the same as the Ro selected by the user);
d1.. D4 is the source rock bottom burial depth at the time (same as the T1.. T4 acquisition method);
l1.. L4 comparing the time determined according to T1.. T4 with the age of the bottom of the drilled formation, if TDj≤Ti≤TDj-1(TiKey time, T, determined for the userDj-1Is TDjThe absolute age of the j-1 th stratum and the stratum layer input by the user in the single well simulation), Li is the name of the j-1 th stratum layer.
R, D is the bottom maturity and burial depth of the source rock today.
Example two:
practical application is carried out aiming at a certain well, the names and the time of layers used in simulation are shown in table 1, the source rock layers P2 and T3, the burial depths of all strata are shown in table 2, maturity history simulation is carried out in the earlier stage and is shown in figure 1, the steps (1) - (4) are applied to obtain data in tables 3 and 4, the step (5) is applied to explain that a template automatically outputs P2 source rock and T3 source rock evolution for information explanation, and the method specifically comprises the following steps:
when a certain well P2 hydrocarbon source rock enters an oil production threshold in the T1f sedimentary period (about 244.35Ma), the maturity of the bottom of the hydrocarbon source rock is 0.5%, and the stratum is buried by 1141 meters deeply; entering a large oil production period at the T2 deposition period (about 236.86Ma), wherein the maturity of the bottom of the source rock is 1%, and the buried depth of the stratum is 2766 m; entering a gas generation threshold in a T3 deposition period (about 230.57Ma), wherein the maturity of the bottom of the source rock is 1.3%, and the buried depth of the stratum is 3062 meters; entering a late gas generation stage at the deposition period of J12 (about 195.21Ma), when the maturity of the bottom of the source rock is 2.5%, and when the stratum is buried 4729 m deep; until now, the bottom maturity of P2 hydrocarbon source rock is 3.21%, and the buried depth is 6932 meters.
A certain well T3 source rock enters the oil producing threshold in the T3x4 deposition period (about 219.18Ma), the maturity of the bottom of the source rock is 0.5%, and the stratum is buried for 829 meters; entering a large oil production period in a J12 deposition period (about 187.87Ma), wherein the maturity of the bottom of the source rock is 1%, and the buried depth of the stratum is 2427 meters; the gas generation threshold is entered in the deposition period of J3 (about 164.65Ma), the maturity of the bottom of the source rock is 1.3%, and the buried depth of the stratum is 3300 m; entering a later gas generation stage in a-1 deposition period (about-1 Ma), wherein the maturity of the bottom of the source rock is-1%, and the stratum is buried by-1 m; until now, the bottom maturity of T3 hydrocarbon source rock is 1.74%, and the burial depth is 4405 meters.
The method for automatically extracting the profile information of the single-well evolution history can configure a data description template according to the requirements of users aiming at different geological properties (buried history, thermal history, maturity history and the like), realizes the automatic, rapid and accurate extraction of the single-well information through computer programming, improves the analysis efficiency and the data acquisition precision of the single-well simulation data, and can be used as an important aid for one-dimensional single-well simulation work.
Table 1, well formation name and absolute age correspondence table:
table 2, well formation name and bottom boundary burial depth correspondence table:
| numbering | Formation of earth | Buried depth of bottom boundary/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 table of key time information of formation maturity evolution of a certain well P2 hydrocarbon source:
| P1 | TY | |||
| a certain | P2 | |||
| R1 | L1 | D1 | T1 | |
| 0.5 | T1f | 1141 | 244.35 | |
| R2 | | 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 table of key time information of formation maturity evolution of a certain well T3 hydrocarbon source:
| P1 | TY | |||
| a certain | T3 | |||
| R1 | L1 | D1 | T1 | |
| 0.5 | T3x4 | 829 | 219.18 | |
| R2 | | 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 |
。
Having described embodiments of the present invention, the foregoing description is intended to be exemplary, not exhaustive, and not 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 described embodiments.
Claims (10)
1. An automatic extraction method for single-well evolution history profile information is characterized by comprising the following steps:
extracting single well data;
forming attribute value isoline burial depth data based on the single well data;
calculating the burial depth and time of an intersection point of the attribute contour value and the buried depth of the target layer based on the attribute value contour buried depth data to obtain buried depth and time data;
and extracting the burial depth and time of set conditions from the burial depth and time data to finish the extraction of the profile information of the single well evolution history.
2. The automatic extraction method for single-well evolution history profile information according to claim 1, wherein the single-well data comprises single-well simulation data and single-well four-history simulation result data;
the single well simulation data comprises a well name, a stratum bottom age, a horizon burial depth and/or a hydrocarbon source rock stratum;
the single well four-history simulation result data comprises the burial depth, the temperature, the maturity and/or the raw oil gas strength of each drill well horizon in the evolution stage.
3. The method for automatically extracting individual well evolution history profile information according to claim 1, wherein the forming attribute value contour buried depth data based on the individual well data comprises:
calculating the burial depth of maturity points of each geological evolution stage on a burial history profile based on the single-well data;
and connecting the maturity points to form a contour line.
4. The method for automatically extracting individual well evolution history profile information according to claim 3, wherein the forming attribute value contour buried depth data based on the individual well data comprises:
setting k evolution times for a single well simulation at a certain time, wherein the corresponding time is t1...tkWherein t iskFor the present moment, suppose a certain geological evolution stage t moment, the t moment is contained in k moments, the stratum participating in the evolution is m layers, and the corresponding stratum burial depth is h1...hmCorresponding maturity of Ro1...RomThe earth surface is 0, the buried depth is positive downwards, and h is setn+1＝0，Ron+1If equal to 0, calculating the Ro value corresponding to the buried depth H at the end of t timetThe method comprises the following specific steps:
from tk-1To t1Judging the Ro value at each moment;
if Ro1≤Ro≤Rom+1Searching from bottom to top, if Ro is satisfiediRo and not less than Roi+1Ro is not more than 1, i is not less than 1 and n is not less than n, the buried depth H is calculated by adopting the following formulat，hiFor formation burial depth, RoiIn order to be of a degree of maturity,
if Ro > Ro1Then search for the first Ro from bottom to topjNot equal to Ro1J, the buried depth H is calculated by the following formulatJ is a natural number,
if Ro is greater than the bottom maturity of the lowest stratum of the well, the corresponding buried depth Ht1...Htk-1All invalid values are represented by-1.
5. The automatic extraction method of single-well evolution history profile information according to claim 4, wherein the calculating of the burial depth and the time of the intersection of the attribute contour value and the target layer burial depth based on the attribute value contour burial depth data to obtain burial depth and time data comprises:
tRj＝ti+(Hti-HLi)/(HLi+1-HLi+Hti-HLi)×(ti-ti+1)，
HRj＝HLi+(Hti-HLi)/(HLi+1-HLi+Hti-HLi)×(HLi+1-HLi)，
Htito be buried deep, HLiAnd HLi+1Acquiring the burial depths of the ith and i +1 evolution stages corresponding to the L stratum through the evolution data of the burial history; t is tiFor the geological evolution stage time, tRjTime of intersection, HRjIs the burial depth of the junction.
6. The automatic extraction method of single-well evolution history profile information according to claim 5, wherein the extracting of the burial depth and time of the set condition in the burial depth and time data comprises:
and obtaining s intersection point burial depths and time for each intersection point, wherein one of the intersection point burial depths and time is geological time and stratum burial depth when the first maturity of the L stratum reaches the corresponding Ro value, extracting a first intersection point value of the isoline and the L stratum according to set conditions, and judging a stratum deposition period corresponding to the moment according to the input stratum age.
7. The automatic extraction method of single-well evolution history profile information according to claim 1, further comprising:
and establishing a data normalization template, and performing normalization processing on parameters in the single-well data, the attribute value isoline burial depth data and the burial depth and time data based on the data normalization template.
8. The automatic extraction method of single-well evolution history profile information according to claim 7, characterized in that, in the data normalization template,
p1 is the well name;
TY is the name of the hydrocarbon source formation;
r1.. R4 is a maturity value of a key time node of the evolution of the source rock;
t1.. T4 is the time for a defined source rock to first reach a key maturity evolutionary node;
d4 is the time source rock bottom burial depth;
l1.. L4 is the time determined according to T1.. T4 compared with the age of the bottom of the drilled formation, if TDj≤Ti≤TDj-1Then L isiIs the name of the stratum at the j-1 st layer, TiKey time, T, determined for the userDj-1Is TDjThe absolute age of the layer j-1 and the layer stratum input in the single well simulation of the user;
r, D is the bottom maturity and burial depth of the source rock today.
9. The utility model provides a single well evolution history section information automatic extraction element which characterized in that includes:
an extraction unit: for extracting single well data;
forming a unit: for forming attribute value isopleth burial depth data based on the single well data;
a calculation unit: the buried depth and time of the intersection point of the attribute contour value and the buried depth of the target layer are calculated based on the attribute value contour buried depth data, and buried depth and time data are obtained;
an information extraction unit: and the method is used for extracting the burial depth and time of set conditions from the burial depth and time data to finish the extraction of the profile information of the single well evolution history.
10. The automatic extraction method of single-well evolution history profile information according to claim 9, further comprising:
a specification unit: the method is used for establishing a data normalization template, and carrying out normalization processing on parameters in single well data, attribute value isoline burial depth data and burial depth and time data based on the data normalization template.
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 true CN113740930A (en) | 2021-12-03 |
| CN113740930B 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 (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140067351A1 (en) * | 2012-08-31 | 2014-03-06 | Gary G. Gray | Method of Estimating Rock Mechanical Properties |
| 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 |
-
2020
- 2020-05-29 CN CN202010482342.7A patent/CN113740930B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140067351A1 (en) * | 2012-08-31 | 2014-03-06 | Gary G. Gray | Method of Estimating Rock Mechanical Properties |
| 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)
| Title |
|---|
| PER AVSETH ETAL: "combining burial history and rock-physics modeling to constrain AVO analysis during exploration", vol. 35, no. 6, pages 474 - 560 * |
| 徐春华;肖永军;邱连贵;: "长岭凹陷断陷期生烃阶段性与成藏圈闭的选择性", 新疆石油地质, no. 03 * |
| 戴世立;杨智;李弘;杨姣;: "鄂尔多斯盆地塔巴庙地区上古生界源岩热演化特征", 海洋石油, no. 01 * |
| 韩杨等: "饶阳凹陷北部洼陷区古近系油气运聚机理与成藏模式", 天然气地球科学, vol. 31, no. 3, pages 358 - 369 * |
| 鲍园;韦重韬;王超勇;曹佳;: "贵州织纳煤田水公河向斜上二叠统8煤层三史模拟", 煤田地质与勘探, no. 06 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113740930B (en) | 2024-02-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20210350208A1 (en) | Method and device for predicting production performance of oil reservoir | |
| CN103370494B (en) | System and method for performing downhole stimulation operations | |
| US9070172B2 (en) | Method and system for data context service | |
| CN101221250A (en) | Three-dimensional stratum characteristic decryption method based on geologic information acquisition data | |
| CN104881435A (en) | Data mining based automatic research flow well logging evaluation expert system | |
| CN102089677A (en) | Well log normalization | |
| CN106547941A (en) | A kind of method that BIM Engineering Exploration Information Systems are built based on revit figure models | |
| CN111950046B (en) | Drilling data model construction method based on BIM | |
| CN108986213B (en) | Three-dimensional stratum modeling method based on stacking technology | |
| CN114542056A (en) | Rapid construction method and device for horizontal well stratum interpretation model | |
| NO20171058A1 (en) | Integrated a priori uncertainty parameter architecture in simulation model creation | |
| 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 | |
| NO20160254A1 (en) | Data analytics for oilfield data repositories | |
| CN109165312A (en) | A kind of survey data management system | |
| CN113534283B (en) | Quantitative evaluation method for ore-forming element characteristics of sandstone-type uranium ores | |
| Albrecht et al. | Laboratory measurements of fluid transport properties on tight gas sandstones and applications | |
| CN114742526A (en) | Standard design process based exploration full-process informatization method and system | |
| CN111767662B (en) | Method and device for correcting azimuth information of single sand body of oil reservoir | |
| Robertson et al. | Optimization of unconventional well-pad area using reservoir simulation and intelligent sequential sampling | |
| CN118016217B (en) | Method for calculating expansion deformation of expansive soil foundation of interbedded structure | |
| CN107060744A (en) | A kind of Logging Geology system | |
| CN117633658B (en) | Rock reservoir lithology identification method and system | |
| Detay et al. | An expert system for groundwater exploration and exploitation | |
| Yan et al. | Integrated Mature Field Management System from Rapid Production Update to History Matching |
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 |