Stratum boundary identification method and device, electronic equipment and storage medium

Abstract

The invention discloses a stratum boundary identification method and device, electronic equipment and a storage medium. The stratum boundary recognition method comprises the steps that a target exploration area is determined, historical logging data corresponding to the target exploration area are obtained, the historical logging data are corrected, corrected logging data are obtained, and the corrected logging data comprise corrected phosphorus content data; determining a limit depth range based on the corrected phosphorus content data, and determining a target survey value corresponding to each candidate limit depth in the limit depth range; a target cross-plot is determined based on the target survey value, and a target formation boundary is determined based on the target cross-plot. Based on the technical scheme of the embodiment of the invention, the stratum boundary identification efficiency can be improved and the stratum boundary identification cost can be reduced on the basis of accurately identifying the stratum interface. Based on the technical scheme of the embodiment of the invention, the efficiency and accuracy of stratum boundary identification are improved, and the cost of stratum boundary identification is reduced.

Description

Technical Field

[0001] The present invention relates to the field of geological research technology, and in particular to a stratum boundary identification method, device, electronic equipment and storage medium. Background Art

[0002] Oil and gas exploration has revealed that the Sinian System has high-quality reservoirs and the Cambrian System has high-quality source rocks, which together form a good source-reservoir combination. They are important oil and gas exploration strata in my country, forming multiple gas-bearing areas of hundreds of billions and trillions of cubic meters. However, it is difficult to identify the stratigraphic boundary between the Sinian System and the Cambrian System, and there are many controversies over the demarcation of the stratigraphic boundary between the Sinian System and the Cambrian System.

[0003] At present, the methods for identifying the stratigraphic boundary between the Sinian and Cambrian systems mainly include: small shell fossils and high-precision carbon isotope characteristic analysis and seismic data analysis. Among them, the former requires a large number of samples, takes a long time and is expensive, and is not suitable for drilling stratigraphic boundary identification. The latter often has large errors and cannot accurately identify stratigraphic boundaries. In summary, there is currently a lack of a high-efficiency, high-accuracy and low-cost method for identifying the stratigraphic boundary between the Sinian and Cambrian systems. Summary of the invention

[0004] The present invention provides a method, device, electronic equipment and storage medium for identifying stratum boundaries, so as to solve the technical problems of low efficiency, poor accuracy and high cost in identifying the stratum boundaries between the Sinian and Cambrian systems.

[0005] According to one aspect of the present invention, a method for identifying stratum boundaries is provided, wherein the method comprises:

[0006] Determine a target survey area, obtain historical well logging data corresponding to the target survey area, and correct the historical well logging data to obtain corrected well logging data, wherein the corrected well logging data includes corrected phosphorus content data;

[0007] Determine a cutoff depth range based on the corrected phosphorus content data, and determine a target survey value corresponding to each candidate cutoff depth in the cutoff depth range;

[0008] A target intersection diagram is determined based on the target survey value, and a target stratigraphic boundary is determined based on the target intersection diagram, wherein the target stratigraphic boundary is the stratigraphic boundary between the Sinian System and the Cambrian System.

[0009] According to another aspect of the present invention, a stratum boundary identification device is provided, wherein the device comprises:

[0010] A data correction module, used to determine a target survey area, obtain historical well logging data corresponding to the target survey area, and correct the historical well logging data to obtain corrected well logging data, wherein the corrected well logging data includes corrected phosphorus content data;

[0011] A data calculation module, used to determine a boundary depth range based on the corrected phosphorus content data, and determine a target survey value corresponding to each candidate boundary depth in the boundary depth range;

[0012] The wiring identification module is used to determine a target intersection diagram based on the target survey value, and determine a target stratigraphic boundary based on the target intersection diagram, wherein the target stratigraphic boundary is the stratigraphic boundary between the Sinian System and the Cambrian System.

[0013] According to another aspect of the present invention, an electronic device is provided, the electronic device comprising:

[0014] at least one processor; and

[0015] a memory communicatively connected to the at least one processor; wherein,

[0016] The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the method for identifying stratigraphic boundaries described in any embodiment of the present invention.

[0017] According to another aspect of the present invention, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the method for identifying stratum boundaries described in any embodiment of the present invention when executed.

[0018] The technical solution of the embodiment of the present invention determines the target survey area, obtains the historical logging data corresponding to the target survey area, corrects the historical logging data, and obtains the corrected logging data, wherein the corrected logging data includes the corrected phosphorus content data; determines the boundary depth range based on the corrected phosphorus content data, and determines the target survey value corresponding to each candidate boundary depth in the boundary depth range; determines the target intersection diagram based on the target survey value, and determines the target stratigraphic boundary based on the target intersection diagram, wherein the target stratigraphic boundary is the stratigraphic boundary between the Sinian System and the Cambrian System. It solves the problem of cumbersome and high cost of using a large number of drilling cores for experimental analysis of small shell fossils, carbon isotopes, etc., and the problem of low accuracy in dividing stratigraphic layers using seismic data, fills the gap in the rapid division of stratigraphic layers using phosphorus elements, and based on the technical solution of the embodiment of the present invention, improves the efficiency and accuracy of stratigraphic boundary identification, and reduces the cost of stratigraphic boundary identification.

[0019] It should be understood that the contents described in this section are not intended to identify the key or important features of the embodiments of the present invention, nor are they intended to limit the scope of the present invention. Other features of the present invention will become easily understood through the following description. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0021] Figure 1 is a flow chart of a method for identifying stratum boundaries provided according to the first embodiment of the present invention;

[0022] Figure 2 is a flow chart of a method for identifying stratum boundaries provided according to Embodiment 2 of the present invention;

[0023] Figure 3 is a target intersection diagram provided according to an embodiment of the present invention;

[0024] Figure 4 is an overall flow chart of a method for identifying stratum boundaries provided according to an embodiment of the present invention;

[0025] Figure 5 is a schematic structural diagram of a stratum boundary identification device provided according to Embodiment 3 of the present invention;

[0026] Figure 6 It is a schematic diagram of the structure of an electronic device for implementing the stratum boundary identification method according to an embodiment of the present invention. DETAILED DESCRIPTION

[0027] In order to enable those skilled in the art to better understand the scheme of the present invention, the technical scheme in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present invention.

[0028] It should be noted that the terms "first", "second", etc. in the specification and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units that are clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0029] Embodiment 1

[0030] Figure 1 A flow chart of a method for identifying stratum boundaries is provided for the first embodiment of the present invention. This embodiment is applicable to geological research. The method can be performed by a stratum boundary identification device. The stratum boundary identification device can be implemented in the form of hardware and / or software. The stratum boundary identification device can be configured in computer software. Figure 1 As shown, the method includes:

[0031] S110, determining a target survey area, obtaining historical well logging data corresponding to the target survey area, and correcting the historical well logging data to obtain corrected well logging data, wherein the corrected well logging data includes corrected phosphorus content data.

[0032] The target survey area may be understood as a survey area corresponding to a survey area where the stratigraphic boundary between the Sinian System and the Cambrian System is to be identified. In an embodiment of the present invention, the target survey area may be set according to scenario requirements and is not specifically limited herein. The survey area may be an area where high-quality oil and gas may be stored. The target survey area may be a partial area of ​​the survey area.

[0033] The historical logging data can be understood as the historical logging data corresponding to the target survey area. In the embodiment of the present invention, the historical logging data can be preset according to the scene requirements and is not specifically limited here. Optionally, the historical logging data may include historical phosphorus content data, historical lithology data, historical carbon isotope data, historical small shell fossil data, and historical logging data corresponding to the surrounding survey area.

[0034] The corrected logging data may be understood as logging data obtained after data correction is performed on the historical logging data. In an embodiment of the present invention, the corrected logging data may be closer to the real logging data than the historical logging data.

[0035] The corrected phosphorus content data may be understood as the corrected phosphorus content data corresponding to the target survey area. Optionally, the corrected phosphorus content data may include a corrected depth phosphorus content curve.

[0036] Optionally, the historical logging data includes historical lithology data and / or historical phosphorus content data, the historical lithology data includes a first historical depth of carbonate rock and / or a second historical depth of clastic rock, and the correcting the historical logging data to obtain the corrected logging data includes:

[0037] Determine actual lithology data of the target survey area, wherein the actual lithology data includes a first actual depth of the carbonate rock and / or a second actual depth of the clastic rock, the carbonate rock is the main lithology of the top of the Sinian System, and the clastic rock is the main lithology of the bottom of the Cambrian System;

[0038] Determine a correction coefficient according to the first historical depth, the second historical depth, the first actual depth, and the second actual depth, wherein the correction coefficient includes a stretching coefficient and / or a depth correction amount;

[0039] The historical logging data is corrected based on the correction coefficient to obtain corrected logging data.

[0040] Wherein, the historical lithology data can be understood as the lithology data of the target survey area recorded in the historical logging data. Optionally, the historical lithology data may include a first historical depth of carbonate rock and / or a second historical depth of clastic rock. Wherein, the carbonate rock is the main lithology at the top of the Sinian System. The clastic rock is the main lithology at the bottom of the Cambrian System. The Cambrian System is closer to the surface of the earth than the Sinian System. The first historical depth is the depth data of the carbonate rock recorded in the historical logging data. The second historical depth is the depth data of the clastic rock recorded in the historical logging data.

[0041] The historical phosphorus content data may be understood as the phosphorus content data corresponding to the target survey area recorded in the historical well logging data.

[0042] The actual lithology data can be understood as the lithology data corresponding to the target survey area actually measured. Optionally, the actual lithology data may include a first actual depth of the carbonate rock and / or a second actual depth of the clastic rock. The first actual depth can be understood as the depth data of the carbonate rock actually measured. The second actual depth can be understood as the depth data of the clastic rock actually measured.

[0043] The correction coefficient can be understood as a parameter used to correct the historical logging data. In the embodiment of the present invention, the correction coefficient can be preset according to the scene requirements and is not specifically limited here. Optionally, the correction coefficient can include a stretching coefficient and / or a depth correction amount.

[0044] Specifically, the correction coefficient is determined according to the first historical depth, the second historical depth, the first actual depth, and the second actual depth. The calculation formula may be:

[0045] H r1 =aH l1 +b,H r2 =aH l2 +b

[0046] Among them, H r1 Indicates the first actual depth, H l1 Indicates the first history depth, H r2 Indicates the second actual depth, H l2 represents the second historical depth, a represents the stretching coefficient, and b represents the depth correction amount.

[0047] S120. Determine a boundary depth range based on the corrected phosphorus content data, and determine a target survey value corresponding to each candidate boundary depth in the boundary depth range.

[0048] Wherein, the boundary depth range can be understood as the depth range where the target stratigraphic boundary is located. Wherein, the target stratigraphic boundary can be understood as the stratigraphic boundary between the Sinian System and the Cambrian System. In an embodiment of the present invention, the boundary depth range is the approximate range where the target stratigraphic boundary is located determined based on the corrected phosphorus content data. Exemplarily, the boundary depth range can be 200m-220m, 200m-260m or 500m-660m, etc.

[0049] The candidate limit depths may be understood as each limit depth included in the limit depth range. For example, when the limit depth range is 200m-220m, the candidate limit depths include 200m, 201m, 202m...219m and 220m.

[0050] The target survey value can be understood as a required value for identifying the target boundary depth among the plurality of candidate boundary depths. In the embodiment of the present invention, the target survey value can be preset according to the scene requirements and is not specifically limited here. Optionally, the target survey value can include an expected value and / or standard deviation corresponding to each candidate boundary depth.

[0051] Optionally, the corrected phosphorus content data includes a depth phosphorus content curve corresponding to the target survey area, and determining the limit depth range based on the corrected phosphorus content data includes:

[0052] Determining the morphological characteristics of the depth phosphorus content curve, wherein the morphological characteristics include at least one of a straight line, a sawtooth shape, a step shape, and a box shape;

[0053] The depth range corresponding to the curve stages of the step type and the box type is determined as the limit depth range.

[0054] The depth phosphorus content curve can be understood as a curve representing the relationship between depth and phosphorus content. Exemplarily, the depth phosphorus content curve is a curve marked in a depth-phosphorus content two-dimensional coordinate system, representing the relationship between the depth and phosphorus content of the target survey area.

[0055] The morphological feature may be understood as a feature of the morphology of the deep phosphorus content curve. In an embodiment of the present invention, the deep phosphorus content curve may include one or more morphological features.

[0056] Optionally, the corrected phosphorus content data further includes adjacent phosphorus content data corresponding to a peripheral survey area adjacent to the target survey area, and the determining of the target survey value corresponding to each candidate boundary depth in the boundary depth range includes:

[0057] For each of the candidate boundary depths in the boundary depth range, a target survey value is determined according to the data amount of the adjacent phosphorus content data and the data value corresponding to each of the adjacent phosphorus content data.

[0058] Wherein, the peripheral survey area can be understood as a survey area adjacent to the target survey area. Optionally, the peripheral survey area can be a circle determined with the target survey area as the center and a preset value as the radius. Wherein, the preset value can be 10m. Optionally, the peripheral survey area can include multiple historical survey areas, and the historical survey area can be a partial area in the peripheral survey area. The historical survey area can be understood as a survey area in which the corresponding corrected phosphorus content data is recorded in the corrected logging data.

[0059] The adjacent phosphorus content data may be understood as the corrected phosphorus content data of the surrounding survey area at the current candidate boundary depth recorded in the corrected logging data.

[0060] In the embodiment of the present invention, the data volume of the adjacent phosphorus content data is the data volume of the adjacent phosphorus content data recorded in the corrected logging data. Since in the embodiment of the present invention, for each candidate boundary depth, one historical survey area corresponds to one adjacent phosphorus content data, the data volume of the adjacent phosphorus content data can also be understood as the number of the historical survey areas recorded in the corrected logging data. Usually, when the preset value is 10m, the corrected logging data records the adjacent phosphorus content data corresponding to 15 historical survey areas, that is, the data volume of the adjacent phosphorus content data is 15.

[0061] In the embodiment of the present invention, the data value corresponding to the adjacent phosphorus content data can be understood as the data value corresponding to each adjacent phosphorus content data recorded in the corrected logging data, that is, the phosphorus content value corresponding to the current historical survey area.

[0062] The target survey value can be understood as the survey value corresponding to each candidate boundary depth required to screen out the target boundary depth from the candidate boundary depths. In the embodiment of the present invention, the target survey value can be preset according to the scene requirements and is not specifically limited here. Optionally, the target survey value may include an expected value and / or a standard deviation.

[0063] Specifically, for each candidate boundary depth in the boundary depth range, according to the data amount of the adjacent phosphorus content data and the data value corresponding to each adjacent phosphorus content data, the calculation formula for determining the expected value may be:

[0064]

[0065] in, represents the expected value, n represents the number of adjacent phosphorus content data, i represents the ith adjacent phosphorus content data, x i Represents the data value corresponding to the i-th adjacent phosphorus content data.

[0066] Specifically, for each candidate boundary depth in the boundary depth range, according to the data amount of the adjacent phosphorus content data and the data value corresponding to each adjacent phosphorus content data, the calculation formula for determining the expected value may be:

[0067]

[0068] Among them, σ represents the expected value, n represents the number of adjacent phosphorus content data, i represents the i-th adjacent phosphorus content data, x i represents the data value corresponding to the i-th adjacent phosphorus content data, Indicates expected value.

[0069] S130. Determine a target intersection diagram based on the target survey value, and determine a target stratigraphic boundary based on the target intersection diagram, wherein the target stratigraphic boundary is the stratigraphic boundary between the Sinian System and the Cambrian System.

[0070] The target intersection map may be understood as an intersection map obtained based on the target survey value corresponding to each candidate boundary depth.

[0071] The target stratigraphic boundary can be understood as the stratigraphic boundary between the Sinian System and the Cambrian System.

[0072] Optionally, the stratum boundary identification method further includes:

[0073] Acquire verification data of the target stratigraphic boundary within a verification depth range based on the corrected logging data, the verification data including carbon isotope data and / or small shell fossil data, the verification depth range including a first depth range in a first direction of the target stratigraphic boundary and / or a second depth range in a second direction of the target stratigraphic boundary, the first direction and the second direction being perpendicular to the ground surface, and the first direction being opposite to the second direction;

[0074] When the verification data meets the verification conditions, the target stratigraphic boundary is determined as the actual stratigraphic boundary, and the verification conditions are that the data difference between the carbon isotope data corresponding to the first depth range and the second depth range respectively exceeds a first difference threshold and the data difference between the small shell fossil data corresponding to the first depth range and the second depth range respectively exceeds a second difference threshold.

[0075] The actual stratigraphic boundary may be understood as the actual stratigraphic boundary between the Sinian System and the Cambrian System.

[0076] In the embodiment of the present invention, for the target stratigraphic boundary determined based on the phosphorus content data, the accuracy is verified by using data of other dimensions, namely, carbon isotope data and small shell fossil data, thereby ensuring the accuracy of the actual stratigraphic boundary determined. The actual stratigraphic boundary is identified based on multi-dimensional data, thereby ensuring the accuracy of stratigraphic boundary identification.

[0077] The technical solution of the embodiment of the present invention determines the target survey area, obtains the historical logging data corresponding to the target survey area, corrects the historical logging data, and obtains the corrected logging data, wherein the corrected logging data includes the corrected phosphorus content data; determines the boundary depth range based on the corrected phosphorus content data, and determines the target survey value corresponding to each candidate boundary depth in the boundary depth range; determines the target intersection diagram based on the target survey value, and determines the target stratigraphic boundary based on the target intersection diagram, wherein the target stratigraphic boundary is the stratigraphic boundary between the Sinian System and the Cambrian System. It solves the problem of cumbersome and high cost of using a large number of drilling cores for experimental analysis of small shell fossils, carbon isotopes, etc., and the problem of low accuracy in dividing stratigraphic layers using seismic data, fills the gap in the rapid division of stratigraphic layers using phosphorus elements, and based on the technical solution of the embodiment of the present invention, improves the efficiency and accuracy of stratigraphic boundary identification, and reduces the cost of stratigraphic boundary identification.

[0078] Embodiment 2

[0079] Figure 2 This is a flow chart of a method for identifying stratum boundaries provided in the second embodiment of the present invention. This embodiment is to refine the method of determining a target intersection diagram based on the target survey value and determining a target stratum boundary based on the target intersection diagram described in the above embodiment. Figure 2 As shown, the method includes:

[0080] S210, determining a target survey area, obtaining historical well logging data corresponding to the target survey area, and correcting the historical well logging data to obtain corrected well logging data, wherein the corrected well logging data includes corrected phosphorus content data.

[0081] S220. Determine a boundary depth range based on the corrected phosphorus content data, and determine a target survey value corresponding to each candidate boundary depth in the boundary depth range.

[0082] S230. For each candidate boundary depth, determine a coordinate point according to the expected value and the standard deviation, and mark the coordinate point in a two-dimensional coordinate system based on a target order to obtain the target intersection map, wherein the target order is from the candidate boundary depth with a large depth value to the candidate boundary depth with a small depth value.

[0083] The coordinate point can be understood as a coordinate point determined based on the expected value and the standard deviation. Among them, A represents the coordinate point, represents the expected value, and σ represents the standard deviation.

[0084] The two-dimensional coordinate system can be understood as A two-dimensional coordinate system with σ as the horizontal axis and σ as the vertical axis.

[0085] The target intersection diagram can be understood as an intersection diagram obtained by marking the coordinate points corresponding to each candidate boundary depth in the two-dimensional coordinate system (refer to Figure 3 ).

[0086] S240: Determine a target stratum boundary based on the target intersection map.

[0087] Optionally, determining the target stratum boundary based on the target intersection diagram includes:

[0088] Determine the coordinate point that meets the preset condition and is marked for the first time as the target point, wherein the preset condition is a coordinate origin that deviates from the two-dimensional coordinate system;

[0089] The candidate boundary depth corresponding to the target point is used as the target boundary depth to obtain the target stratum boundary.

[0090] The target point may be understood as the coordinate point that is first marked in the target intersection diagram and deviates from the coordinate origin of the two-dimensional coordinate system.

[0091] The target boundary depth may be understood as the boundary depth at which the target formation boundary is located.

[0092] The technical solution of the embodiment of the present invention determines the coordinate points according to the expected value and the standard deviation for each candidate boundary depth, and marks the coordinate points in a two-dimensional coordinate system based on the target order to obtain the target intersection map, wherein the target order is from the candidate boundary depth with a large depth value to the candidate boundary depth with a small depth value. The accuracy of the obtained target intersection map is guaranteed, and the accuracy of the target boundary depth and the target stratum boundary determined based on the target intersection map is further improved.

[0093] Figure 4 FIG. 1 is an overall flow chart of a method for identifying stratum boundaries provided according to an embodiment of the present invention. Figure 4 As shown, the overall process of the stratigraphic boundary identification method can be:

[0094] 1. Investigate the geological background, obtain historical logging data, and use natural gamma logging curves to correct the logging data.

[0095] First, investigate the structure and sedimentary background of the target survey area to understand the development characteristics of the main lithologies of the Sinian-Cambrian systems. For example, the top of the Sinian system is mainly carbonate rock, and the bottom of the Cambrian system is mainly clastic rock. Collect the logging data of the drilling wells, and compare the actual observed lithology with the lithology reflected by the natural gamma logging curve in the logging data, and establish the depth correspondence formula (1) H r =aHl +b, where H r Indicates the actual depth, H l represents the historical depth, a represents the stretching coefficient, and b represents the depth correction. The depth of other logging curves is adjusted using the current relationship formula to perform depth correction.

[0096] 2. Analyze the phosphorus content data curve, calculate the target survey value, and obtain the target intersection diagram.

[0097] According to the corrected phosphorus content data curve obtained in the relationship formula (1), the morphological characteristics (straight line, sawtooth, step, box) and stability of the phosphorus element curve are determined, and the depth range (limit depth range) where the phosphorus content morphology is step-shaped or box-shaped is selected. And the relationship formula (3) Calculate the expected value and standard deviation of a certain number of phosphorus content data corresponding to each depth (candidate boundary depth) respectively. Put the expected value and standard deviation data obtained at each depth into a two-dimensional scatter plot to obtain Intersection diagram (target intersection diagram).

[0098] 3. Analyze the target intersection map to obtain the stratigraphic boundary between the Sinian and Cambrian systems.

[0099] based on Intersection diagram, find the target point that is obviously deviated from the coordinate origin from the intersection point (coordinate point), and record the candidate boundary depth corresponding to it, from deep to shallow. The candidate boundary depth corresponding to the first data deviating from the coordinate origin should be the Sinian-Cambrian boundary (target boundary depth).

[0100] 4. Verify the stratigraphic boundary between the Sinian and Cambrian systems by combining isotope data and / or small shell fossil data.

[0101] By comparing the characteristics of other logging curves near the target boundary depth, there should be obvious differences in the logging curve characteristics above and below the depth; at the same time, the characteristics of small shell fossils and carbon isotopes in the range near the depth can be analyzed. This depth should be the boundary for the development of small shell fossils. The carbon isotope composition also changes significantly. If inconsistencies occur, it is analyzed based on the geological background to determine whether the erosion of the strata is caused by tectonic activities and some typical phenomena are missing.

[0102] The present invention provides a method for quickly and accurately determining the boundary between the Sinian and Cambrian systems based on the phosphorus content of element logging, which solves the problem of cumbersome and high cost of using a large number of drilling cores for experimental analysis of small shell fossils, carbon isotopes, etc., and the problem of low accuracy in using seismic data to divide strata. At the same time, it fills the gap in using phosphorus to quickly divide strata, significantly improves work efficiency, saves costs, and provides technical and theoretical support for improving the benefits of oil and gas exploration.

[0103] Embodiment 3

[0104] Figure 5 This is a schematic diagram of the structure of a stratum boundary identification device provided in the third embodiment of the present invention. Figure 5 As shown, the device includes: a data correction module 310, a data calculation module 320 and a boundary identification module 330; wherein,

[0105] The data correction module 310 is used to determine the target exploration area, obtain the historical logging data corresponding to the target exploration area, and correct the historical logging data to obtain the corrected logging data, wherein the corrected logging data includes the corrected phosphorus content data; the data calculation module 320 is used to determine the boundary depth range based on the corrected phosphorus content data, and determine the target exploration value corresponding to each candidate boundary depth in the boundary depth range; the boundary identification module 330 is used to determine the target intersection diagram based on the target exploration value, and determine the target stratigraphic boundary based on the target intersection diagram, wherein the target stratigraphic boundary is the stratigraphic boundary between the Sinian System and the Cambrian System.

[0106] The technical solution of the embodiment of the present invention determines the target survey area, obtains the historical logging data corresponding to the target survey area, corrects the historical logging data, and obtains the corrected logging data, wherein the corrected logging data includes the corrected phosphorus content data; determines the boundary depth range based on the corrected phosphorus content data, and determines the target survey value corresponding to each candidate boundary depth in the boundary depth range; determines the target intersection diagram based on the target survey value, and determines the target stratigraphic boundary based on the target intersection diagram, wherein the target stratigraphic boundary is the stratigraphic boundary between the Sinian System and the Cambrian System. It solves the problem of cumbersome and high cost of using a large number of drilling cores for experimental analysis of small shell fossils, carbon isotopes, etc., and the problem of low accuracy in dividing stratigraphic layers using seismic data, fills the gap in the rapid division of stratigraphic layers using phosphorus elements, and based on the technical solution of the embodiment of the present invention, improves the efficiency and accuracy of stratigraphic boundary identification, and reduces the cost of stratigraphic boundary identification.

[0107] Optionally, the historical logging data includes historical lithology data and / or historical phosphorus content data, and the historical lithology data includes a first historical depth of carbonate rock and / or a second historical depth of clastic rock. The data correction module 310 is used to:

[0108] Determine actual lithology data of the target survey area, wherein the actual lithology data includes a first actual depth of the carbonate rock and / or a second actual depth of the clastic rock, the carbonate rock is the main lithology of the top of the Sinian System, and the clastic rock is the main lithology of the bottom of the Cambrian System;

[0109] Determine a correction coefficient according to the first historical depth, the second historical depth, the first actual depth, and the second actual depth, wherein the correction coefficient includes a stretching coefficient and / or a depth correction amount;

[0110] The historical logging data is corrected based on the correction coefficient to obtain corrected logging data.

[0111] Optionally, the data calculation module 320 includes: a morphological feature recognition unit and a boundary range determination unit; wherein,

[0112] The morphological feature recognition unit is used to determine the morphological features of the depth phosphorus content curve, wherein the morphological features include at least one of a straight line, a sawtooth, a step, and a box;

[0113] The limit range determination unit is used to determine the depth range corresponding to the curve stages of the step type and the box type as the limit depth range.

[0114] Optionally, the corrected phosphorus content data also includes adjacent phosphorus content data corresponding to a peripheral survey area adjacent to the target survey area. The data calculation module 320 includes: a data calculation unit, used to determine the target survey value for each candidate boundary depth in the boundary depth range according to the data amount of the adjacent phosphorus content data and the data value corresponding to each adjacent phosphorus content data.

[0115] Optionally, the target survey value includes an expected value and / or a standard deviation. The boundary identification module 330 includes: an intersection diagram determination unit, configured to determine a coordinate point according to the expected value and the standard deviation for each candidate boundary depth, and mark the coordinate point in a two-dimensional coordinate system based on a target order to obtain the target intersection diagram.

[0116] The target sequence is from the candidate boundary depth with a larger depth value to the candidate boundary depth with a smaller depth value.

[0117] Optionally, the boundary recognition module 330 includes: a coordinate determination unit and a boundary recognition unit; wherein,

[0118] The coordinate determination unit is used to determine the coordinate point that meets the preset condition and is marked for the first time as the target point, wherein the preset condition is a coordinate origin that deviates from the two-dimensional coordinate system;

[0119] The boundary recognition unit is used to use the candidate boundary depth corresponding to the target point as the target boundary depth to obtain the target stratum boundary.

[0120] Optionally, the stratum boundary identification device further includes: a verification preparation module and a boundary verification module; wherein,

[0121] The verification preparation module is used to obtain verification data of the target stratigraphic boundary within a verification depth range based on the corrected logging data, the verification data including carbon isotope data and / or small shell fossil data, the verification depth range including a first depth range in a first direction of the target stratigraphic boundary and / or a second depth range in a second direction of the target stratigraphic boundary, the first direction and the second direction are perpendicular to the ground surface, and the first direction is opposite to the second direction;

[0122] The boundary verification module is used to determine the target stratigraphic boundary as the actual stratigraphic boundary when the verification data meets the verification conditions, and the verification conditions are that the data difference between the carbon isotope data corresponding to the first depth range and the second depth range respectively exceeds a first difference threshold and the data difference between the small shell fossil data corresponding to the first depth range and the second depth range respectively exceeds a second difference threshold.

[0123] The stratum boundary identification device provided in the embodiment of the present invention can execute the stratum boundary identification method provided in any embodiment of the present invention, and has the corresponding functional modules and beneficial effects of the execution method.

[0124] Embodiment 4

[0125] Figure 6 A schematic diagram of the structure of an electronic device 10 that can be used to implement an embodiment of the present invention is shown. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as personal digital processing, cellular phones, smart phones, wearable devices (such as helmets, glasses, watches, etc.) and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the present invention described and / or required herein.

[0126] like Figure 6As shown, the electronic device 10 includes at least one processor 11, and a memory connected to the at least one processor 11, such as a read-only memory (ROM) 12, a random access memory (RAM) 13, etc., wherein the memory stores a computer program that can be executed by at least one processor, and the processor 11 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 12 or the computer program loaded from the storage unit 18 to the random access memory (RAM) 13. In the RAM 13, various programs and data required for the operation of the electronic device 10 can also be stored. The processor 11, the ROM 12, and the RAM 13 are connected to each other through a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

[0127] A number of components in the electronic device 10 are connected to the I / O interface 15, including: an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, speakers, etc.; a storage unit 18, such as a disk, an optical disk, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices through a computer network such as the Internet and / or various telecommunication networks.

[0128] The processor 11 may be a variety of general and / or special processing components with processing and computing capabilities. Some examples of the processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any appropriate processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as a formation boundary identification method.

[0129] In some embodiments, the formation boundary identification method may be implemented as a computer program, which is tangibly contained in a computer-readable storage medium, such as a storage unit 18. In some embodiments, part or all of the computer program may be loaded and / or installed on the electronic device 10 via the ROM 12 and / or the communication unit 19. When the computer program is loaded into the RAM 13 and executed by the processor 11, one or more steps of the formation boundary identification method described above may be performed. Alternatively, in other embodiments, the processor 11 may be configured to perform the formation boundary identification method in any other appropriate manner (e.g., by means of firmware).

[0130] Various implementations of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), systems on chips (SOCs), load programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various implementations can include: being implemented in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which can be a special purpose or general purpose programmable processor that can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit data and instructions to the storage system, the at least one input device, and the at least one output device.

[0131] Computer programs for implementing the methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, so that when the computer program is executed by the processor, the functions / operations specified in the flow chart and / or block diagram are implemented. The computer program may be executed entirely on the machine, partially on the machine, partially on the machine and partially on a remote machine as a stand-alone software package, or entirely on a remote machine or server.

[0132] In the context of the present invention, a computer-readable storage medium may be a tangible medium that may contain or store a computer program for use by or in combination with an instruction execution system, device or equipment. A computer-readable storage medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, devices or equipment, or any suitable combination of the foregoing. Alternatively, a computer-readable storage medium may be a machine-readable signal medium. A more specific example of a machine-readable storage medium may include an electrical connection based on one or more lines, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0133] To provide interaction with a user, the systems and techniques described herein may be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and a pointing device (e.g., a mouse or trackball) through which the user can provide input to the electronic device. Other types of devices may also be used to provide interaction with the user; for example, the feedback provided to the user may be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user may be received in any form (including acoustic input, voice input, or tactile input).

[0134] The systems and techniques described herein may be implemented in a computing system that includes backend components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes frontend components (e.g., a user computer with a graphical user interface or a web browser through which a user can interact with implementations of the systems and techniques described herein), or a computing system that includes any combination of such backend components, middleware components, or frontend components. The components of the system may be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: a local area network (LAN), a wide area network (WAN), a blockchain network, and the Internet.

[0135] A computing system may include a client and a server. The client and the server are generally remote from each other and usually interact through a communication network. The client and server relationship is generated by computer programs running on the corresponding computers and having a client-server relationship with each other. The server may be a cloud server, also known as a cloud computing server or cloud host, which is a host product in the cloud computing service system to solve the defects of difficult management and weak business scalability in traditional physical hosts and VPS services.

[0136] It should be understood that the various forms of processes shown above can be used to reorder, add or delete steps. For example, the steps described in the present invention can be executed in parallel, sequentially or in different orders, as long as the desired results of the technical solution of the present invention can be achieved, and this document does not limit this.

[0137] The above specific implementations do not constitute a limitation on the protection scope of the present invention. It should be understood by those skilled in the art that various modifications, combinations, sub-combinations and substitutions can be made according to design requirements and other factors. Any modification, equivalent substitution and improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims

1. A method for identifying stratum boundaries, characterized in that: include: Determine a target survey area, obtain historical well logging data corresponding to the target survey area, and correct the historical well logging data to obtain corrected well logging data, wherein the corrected well logging data includes corrected phosphorus content data; Determine a cutoff depth range based on the corrected phosphorus content data, and determine a target survey value corresponding to each candidate cutoff depth in the cutoff depth range; A target intersection diagram is determined based on the target survey value, and a target stratigraphic boundary is determined based on the target intersection diagram, wherein the target stratigraphic boundary is the stratigraphic boundary between the Sinian System and the Cambrian System.

2. The method according to claim 1, characterized in that The historical well logging data includes historical lithology data and / or historical phosphorus content data, the historical lithology data includes a first historical depth of carbonate rock and / or a second historical depth of clastic rock, and the correcting of the historical well logging data to obtain the corrected well logging data includes: Determine actual lithology data of the target survey area, wherein the actual lithology data includes a first actual depth of the carbonate rock and / or a second actual depth of the clastic rock, the carbonate rock is the main lithology of the top of the Sinian System, and the clastic rock is the main lithology of the bottom of the Cambrian System; Determine a correction coefficient according to the first historical depth, the second historical depth, the first actual depth, and the second actual depth, wherein the correction coefficient includes a stretching coefficient and / or a depth correction amount; The historical logging data is corrected based on the correction coefficient to obtain corrected logging data.

3. The method according to claim 1, characterized in that The corrected phosphorus content data includes a depth phosphorus content curve corresponding to the target survey area, and determining the limit depth range based on the corrected phosphorus content data includes: Determining the morphological characteristics of the depth phosphorus content curve, wherein the morphological characteristics include at least one of a straight line, a sawtooth shape, a step shape, and a box shape; The depth range corresponding to the curve stages of the step type and the box type is determined as the limit depth range.

4. The method according to claim 1, characterized in that: The corrected phosphorus content data also includes adjacent phosphorus content data corresponding to a peripheral survey area adjacent to the target survey area, and the determining of a target survey value corresponding to each candidate boundary depth in the boundary depth range includes: For each of the candidate boundary depths in the boundary depth range, a target survey value is determined according to the data amount of the adjacent phosphorus content data and the data value corresponding to each of the adjacent phosphorus content data.

5. The method according to claim 1, characterized in that The target survey value includes an expected value and / or a standard deviation, and determining a target intersection diagram based on the target survey value includes: For each candidate boundary depth, a coordinate point is determined according to the expected value and the standard deviation, and the coordinate point is marked in a two-dimensional coordinate system based on the target order to obtain the target intersection map, The target sequence is from the candidate boundary depth with a larger depth value to the candidate boundary depth with a smaller depth value.

6. The method according to claim 5, characterized in that Determining the target stratum boundary based on the target cross-plot includes: Determine the coordinate point that meets the preset condition and is marked for the first time as the target point, wherein the preset condition is a coordinate origin that deviates from the two-dimensional coordinate system; The candidate boundary depth corresponding to the target point is used as the target boundary depth to obtain the target stratum boundary.

7. The method according to claim 1, characterized in that Also includes: Acquire verification data of the target stratigraphic boundary within a verification depth range based on the corrected logging data, the verification data including carbon isotope data and / or small shell fossil data, the verification depth range including a first depth range in a first direction of the target stratigraphic boundary and / or a second depth range in a second direction of the target stratigraphic boundary, the first direction and the second direction being perpendicular to the ground surface, and the first direction being opposite to the second direction; When the verification data meets the verification conditions, the target stratigraphic boundary is determined as the actual stratigraphic boundary, and the verification conditions are that the data difference between the carbon isotope data corresponding to the first depth range and the second depth range respectively exceeds a first difference threshold and the data difference between the small shell fossil data corresponding to the first depth range and the second depth range respectively exceeds a second difference threshold.

8. A stratum boundary identification device, characterized in that: include: A data correction module, used to determine a target survey area, obtain historical well logging data corresponding to the target survey area, and correct the historical well logging data to obtain corrected well logging data, wherein the corrected well logging data includes corrected phosphorus content data; A data calculation module, used to determine a boundary depth range based on the corrected phosphorus content data, and determine a target survey value corresponding to each candidate boundary depth in the boundary depth range; The boundary identification module is used to determine a target intersection diagram based on the target survey value, and determine a target stratigraphic boundary based on the target intersection diagram, wherein the target stratigraphic boundary is the stratigraphic boundary between the Sinian System and the Cambrian System.

9. An electronic device, characterized in that: The electronic device comprises: at least one processor; and a memory communicatively connected to the at least one processor; wherein, The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can perform the method for identifying stratum boundaries according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the stratigraphic boundary identification method according to any one of claims 1 to 7 when executed.

Images