Stratum boundary identification method and device, electronic equipment and storage medium
By acquiring and correcting phosphorus content data in historical well logging data, and combining cross plots and multi-dimensional data verification, the boundary between the Sinian and Cambrian strata was identified, solving the problems of low efficiency, poor accuracy and high cost in existing technologies, and achieving efficient and low-cost boundary identification.
Patent Information
- Application Number
- CN202311102524.7
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2023-08-29
- Publication Date
- 2026-01-23
- Estimated Expiration
- 2043-08-29
AI Technical Summary
Existing methods for identifying the stratigraphic boundary between the Sinian and Cambrian systems are inefficient, inaccurate, and costly, especially the analysis of small shell fossils, which is time-consuming and labor-intensive, and the analysis of seismic data has large errors.
By acquiring historical logging data of the target exploration area, corrected phosphorus content data is obtained after calibration. The target exploration values of the boundary depth range and candidate boundary depth are determined. The stratigraphic boundaries are identified using the target cross-plot, and verified by carbon isotope and small shell fossil data.
It improves the efficiency and accuracy of stratigraphic boundary identification, reduces identification costs, and fills the gap in using phosphorus for rapid stratigraphic delineation.
Smart Images

Figure CN119535577B_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of geological research technology, and in particular to a method, apparatus, electronic device, and storage medium for identifying stratigraphic boundaries. Background Technology
[0002] Oil and gas exploration has revealed that the Sinian System contains high-quality reservoirs and the Cambrian System contains high-quality source rocks, forming a favorable source-reservoir combination. This is an important oil and gas exploration system in my country, forming several gas-bearing areas with reserves of hundreds of billions and trillions of cubic meters. However, identifying the stratigraphic boundary between the Sinian and Cambrian Systems presents certain difficulties, and the delineation of this boundary remains controversial.
[0003] Currently, the main methods for identifying the stratigraphic boundary between the Sinian and Cambrian systems include: analysis of small shell fossils and high-precision carbon isotope characteristics, and seismic data analysis. The former requires a large number of samples and is time-consuming and expensive, making it unsuitable for well drilling for stratigraphic boundary identification. The latter often has large errors and cannot accurately identify the stratigraphic boundary. In summary, there is currently a lack of a highly efficient, accurate, and low-cost method for identifying the Sinian-Cambrian stratigraphic boundary. Summary of the Invention
[0004] This invention provides a method, apparatus, electronic device, and storage medium for identifying stratigraphic boundaries, in order to solve the technical problems of low efficiency, poor accuracy, and high cost in identifying the Sinian-Cambrian stratigraphic boundary.
[0005] According to one aspect of the present invention, a method for identifying stratigraphic boundaries is provided, wherein the method includes:
[0006] The target exploration area is determined, the historical logging data corresponding to the target exploration area is obtained, and the historical logging data is corrected to obtain corrected logging data, wherein the corrected logging data includes corrected phosphorus content data;
[0007] Based on the corrected phosphorus content data, a limit depth range is determined, and the target survey value corresponding to each candidate limit depth in the limit depth range is determined.
[0008] Based on the target survey values, a target intersection map is determined, and based on the target intersection map, the target stratigraphic boundary is determined, 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 stratigraphic boundary identification device is provided, wherein the device comprises:
[0010] The data correction module is used to determine the target exploration area, acquire the historical logging data corresponding to the target exploration area, correct the historical logging data, and obtain corrected logging data, wherein the corrected logging data includes corrected phosphorus content data.
[0011] The data calculation module is used to determine the boundary depth range based on the corrected phosphorus content data, and to determine the target survey value corresponding to each candidate boundary depth in the boundary depth range.
[0012] The wiring identification module is used to determine the target intersection map based on the target survey value, and to determine the target stratigraphic boundary based on the target intersection map, 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 that can be executed by the at least one processor, which enables the at least one processor to perform the stratigraphic boundary identification method according to any embodiment of the present invention.
[0017] According to another aspect of the present invention, a computer-readable storage medium is provided, the computer-readable storage medium storing computer instructions for causing a processor to execute and implement the stratigraphic boundary identification method according to any embodiment of the present invention.
[0018] The technical solution of this invention involves determining a target exploration area, acquiring historical logging data corresponding to the target exploration area, correcting the historical logging data to obtain corrected logging data, wherein the corrected logging data includes corrected phosphorus content data; determining a boundary depth range based on the corrected phosphorus content data, and determining a target exploration value corresponding to each candidate boundary depth within the boundary depth range; determining a target cross-plot based on the target exploration value, and determining a target stratigraphic boundary based on the target cross-plot, wherein the target stratigraphic boundary is the stratigraphic boundary between the Sinian and Cambrian systems. This solves the problems of cumbersome and costly experimental analysis using a large number of drill cores for small-shell fossils, carbon isotopes, etc., and the low accuracy of using seismic data for stratigraphic division. It fills the gap in rapidly dividing stratigraphic layers using phosphorus. Based on the technical solution of this invention, the efficiency and accuracy of stratigraphic boundary identification are improved, and the cost of stratigraphic boundary identification is reduced.
[0019] It should be understood that the description in this section is not intended to identify key or essential features of the embodiments of the present invention, nor is it intended to limit the scope of the invention. Other features of the invention will become readily apparent from the following description. Attached Figure Description
[0020] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.
[0021] Figure 1 This is a flowchart of a stratigraphic boundary identification method provided in Embodiment 1 of the present invention;
[0022] Figure 2 This is a flowchart of a stratigraphic boundary identification method provided in Embodiment 2 of the present invention;
[0023] Figure 3 This is a target intersection diagram provided according to an embodiment of the present invention;
[0024] Figure 4 This is an overall flowchart of a stratigraphic boundary identification method provided by an embodiment of the present invention;
[0025] Figure 5 This is a schematic diagram of the structure of a stratigraphic boundary identification device provided in Embodiment 3 of the present invention;
[0026] Figure 6 This is a schematic diagram of the structure of an electronic device that implements the stratigraphic boundary identification method of this invention. Detailed Implementation
[0027] To enable those skilled in the art to better understand the present invention, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort 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, claims, and accompanying drawings of this invention are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of the invention described herein can be implemented in orders other than those illustrated or described herein. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover a non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.
[0029] Example 1
[0030] Figure 1 This is a flowchart illustrating a stratigraphic boundary identification method according to Embodiment 1 of the present invention. This embodiment is applicable to geological research. The method can be executed by a stratigraphic boundary identification device, which can be implemented in hardware and / or software and can be configured in computer software. Figure 1 As shown, the method includes:
[0031] S110. Determine the target exploration area, obtain the historical logging data corresponding to the target exploration area, and correct the historical logging data to obtain corrected logging data, wherein the corrected logging data includes corrected phosphorus content data.
[0032] The target exploration area can be understood as the exploration area corresponding to the stratigraphic boundary between the Sinian and Cambrian systems to be identified. In this embodiment of the invention, the target exploration area can be set according to the needs of the scenario, and is not specifically limited here. The exploration area can be a region that may contain high-quality oil and gas. The target exploration area can be a portion of the exploration area.
[0033] The historical logging data can be understood as the historical logging data corresponding to the target exploration area. In this embodiment of the invention, the historical logging data can be preset according to scenario 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 surrounding exploration areas, etc.
[0034] The corrected logging data can be understood as logging data obtained after correcting the historical logging data. In this embodiment of the invention, the corrected logging data can be closer to the actual logging data than the historical logging data.
[0035] The corrected phosphorus content data can be understood as the corrected phosphorus content data corresponding to the target survey area. Optionally, the corrected phosphorus content data may include the corrected depth phosphorus content curve.
[0036] Optionally, the historical logging data includes historical lithological data and / or historical phosphorus content data, wherein the historical lithological data includes the first historical depth of carbonate rocks and / or the second historical depth of clastic rocks, and the step of correcting the historical logging data to obtain corrected logging data includes:
[0037] The actual lithological data of the target survey area are determined, wherein the actual lithological data includes the first actual depth of the carbonate rocks and / or the second actual depth of the clastic rocks, wherein the carbonate rocks are the main lithology of the upper Sinian system and the clastic rocks are the main lithology of the lower Cambrian system;
[0038] Based on the first historical depth, the second historical depth, the first actual depth, and the second actual depth, a correction coefficient is determined, 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] The historical lithological data can be understood as the lithological data of the target exploration area recorded in the historical well logging data. Optionally, the historical lithological data may include a first historical depth of carbonate rocks and / or a second historical depth of clastic rocks. The carbonate rocks are the main lithology of the upper Sinian system. The clastic rocks are the main lithology of the lower Cambrian system. The Cambrian system is closer to the Earth's surface than the Sinian system. The first historical depth is the depth data of the carbonate rocks recorded in the historical well logging data. The second historical depth is the depth data of the clastic rocks recorded in the historical well logging data.
[0041] The historical phosphorus content data can be understood as the phosphorus content data corresponding to the target exploration area recorded in the historical well logging data.
[0042] The actual lithological data can be understood as the lithological data corresponding to the target survey area that was actually measured. Optionally, the actual lithological 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 that was actually measured. The second actual depth can be understood as the depth data of the clastic rock that was actually measured.
[0043] The correction coefficient can be understood as a parameter used to correct the historical logging data. In this embodiment of the invention, the correction coefficient can be preset according to scenario requirements, and is not specifically limited here. Optionally, the correction coefficient may include a stretching coefficient and / or a depth correction amount.
[0044] Specifically, the correction coefficient is determined based on 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 H represents the first actual depth. l1 H represents the first historical depth. r2 H represents the second actual depth. l2 This indicates the second historical depth, where 'a' represents the stretching factor and 'b' represents the depth correction amount.
[0047] S120. Determine the boundary depth range based on the corrected phosphorus content data, and determine the target survey value corresponding to each candidate boundary depth in the boundary depth range.
[0048] The boundary depth range can be understood as the depth range within which the target stratigraphic boundary is located. Specifically, the target stratigraphic boundary can be understood as the stratigraphic boundary between the Sinian and Cambrian systems. In this embodiment of the invention, the boundary depth range is the approximate range within which the target stratigraphic boundary is located, determined based on the corrected phosphorus content data. For example, the boundary depth range can be 200m-220m, 200m-260m, or 500m-660m, etc.
[0049] The candidate limit depth can 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 the required value for identifying the target boundary depth among a plurality of candidate boundary depths. In this embodiment of the invention, the target survey value can be preset according to scenario requirements, and is not specifically limited here. Optionally, the target survey value may include the 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 boundary depth range based on the corrected phosphorus content data includes:
[0052] The morphological characteristics of the depth phosphorus content curve are determined, wherein the morphological characteristics include at least one of linear, sawtooth, step-shaped, and box-shaped.
[0053] The depth range corresponding to the curved stages of the stepped and box-shaped structures is determined as the limit depth range.
[0054] The depth phosphorus content curve can be understood as a curve characterizing the relationship between depth and phosphorus content. For example, the depth phosphorus content curve is a curve plotted in a depth-phosphorus content two-dimensional coordinate system, characterizing the relationship between depth and phosphorus content in the target survey area.
[0055] The morphological features can be understood as the morphological characteristics of the depth phosphorus content curve. In embodiments of the present invention, the depth 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 surrounding survey areas adjacent to the target survey area, and determining the target survey value corresponding to each candidate boundary depth in the boundary depth range includes:
[0057] For each candidate boundary depth within the boundary depth range, the target survey value is determined based on the amount of adjacent phosphorus content data and the corresponding data value of each adjacent phosphorus content data.
[0058] The surrounding exploration area can be understood as the exploration area adjacent to the target exploration area. Optionally, the surrounding exploration area can be a circle defined with the target exploration area as the center and a preset value as the radius. The preset value can be 10m. Optionally, the surrounding exploration area can include multiple historical exploration areas, and the historical exploration area can be a portion of the surrounding exploration area. The historical exploration area can be understood as the exploration area where the corresponding corrected phosphorus content data is recorded in the corrected logging data.
[0059] The adjacent phosphorus content data can be understood as the corrected phosphorus content data of the surrounding exploration area at the current candidate boundary depth, recorded in the corrected logging data.
[0060] In this embodiment of the invention, the amount of adjacent phosphorus content data is the amount of adjacent phosphorus content data recorded in the calibration logging data. Since in this embodiment, for each candidate boundary depth, one historical exploration area corresponds to one adjacent phosphorus content data point, the amount of adjacent phosphorus content data can also be understood as the number of historical exploration areas recorded in the calibration logging data. Typically, when the preset value is 10m, the calibration logging data records adjacent phosphorus content data corresponding to 15 historical exploration areas; that is, the amount of adjacent phosphorus content data is 15.
[0061] In this embodiment of the invention, the data value corresponding to the adjacent phosphorus content data can be understood as the data value corresponding to each of the adjacent phosphorus content data recorded in the calibration logging data, that is, the value of the phosphorus content corresponding to the current historical exploration area.
[0062] The target survey value can be understood as the survey value required to select the target boundary depth from the candidate boundary depths, corresponding to each candidate boundary depth. In this embodiment of the invention, the target survey value can be preset according to scenario requirements, and is not specifically limited here. Optionally, the target survey value may include the expected value and / or standard deviation.
[0063] Specifically, the formula for determining the expected value for each candidate boundary depth within the boundary depth range, based on the amount of adjacent phosphorus content data and the corresponding data value of each adjacent phosphorus content data, can be:
[0064]
[0065] in, Let x represent the expected value, n represent the number of adjacent phosphorus content data points, i represent the i-th adjacent phosphorus content data point, and x represent the expected value. i This represents the data value corresponding to the i-th adjacent phosphorus content data.
[0066] Specifically, the formula for determining the expected value for each candidate boundary depth within the boundary depth range, based on the amount of adjacent phosphorus content data and the corresponding data value of each adjacent phosphorus content data, can be:
[0067]
[0068] Where σ represents the expected value, n represents the number of adjacent phosphorus content data points, i represents the i-th adjacent phosphorus content data point, and x i This represents the data value corresponding to the i-th adjacent phosphorus content data. This represents the expected value.
[0069] S130. Determine the target intersection map based on the target survey values, and determine the target stratigraphic boundary based on the target intersection map, wherein the target stratigraphic boundary is the stratigraphic boundary between the Sinian System and the Cambrian System.
[0070] The target intersection map can be understood as an intersection map obtained based on the target survey value corresponding to the depth of each candidate boundary.
[0071] The target stratigraphic boundary can be understood as the stratigraphic boundary between the Sinian System and the Cambrian System.
[0072] Optionally, the stratigraphic boundary identification method further includes:
[0073] Based on the corrected logging data, verification data of the target formation boundary within the verification depth range is obtained. The verification data includes carbon isotope data and / or small shell fossil data. The verification depth range includes a first depth range in a first direction of the target formation boundary and / or a second depth range in a second direction of the target formation boundary. The first direction and the second direction are perpendicular to the ground surface, and the first direction is opposite to the second direction.
[0074] If the verification data meets the verification conditions, the target stratigraphic boundary is determined as the actual stratigraphic boundary. The verification conditions are that the data difference between the carbon isotope data corresponding to the first depth range and the second depth range 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 exceeds a second difference threshold.
[0075] The actual stratigraphic boundary can be understood as the actual stratigraphic boundary between the Sinian System and the Cambrian System.
[0076] In this embodiment of the invention, the accuracy of the determined actual stratigraphic boundary is verified using other dimensions of data, namely carbon isotope data and small shell fossil data, to ensure the accuracy of the determined actual stratigraphic boundary. Identifying the actual stratigraphic boundary based on multi-dimensional data ensures the accuracy of the stratigraphic boundary identification.
[0077] The technical solution of this invention involves determining a target exploration area, acquiring historical logging data corresponding to the target exploration area, correcting the historical logging data to obtain corrected logging data, wherein the corrected logging data includes corrected phosphorus content data; determining a boundary depth range based on the corrected phosphorus content data, and determining a target exploration value corresponding to each candidate boundary depth within the boundary depth range; determining a target cross-plot based on the target exploration value, and determining a target stratigraphic boundary based on the target cross-plot, wherein the target stratigraphic boundary is the stratigraphic boundary between the Sinian and Cambrian systems. This solves the problems of cumbersome and costly experimental analysis using a large number of drill cores for small-shell fossils, carbon isotopes, etc., and the low accuracy of using seismic data for stratigraphic division. It fills the gap in rapidly dividing stratigraphic layers using phosphorus. Based on the technical solution of this invention, the efficiency and accuracy of stratigraphic boundary identification are improved, and the cost of stratigraphic boundary identification is reduced.
[0078] Example 2
[0079] Figure 2 This is a flowchart of a stratigraphic boundary identification method provided in Embodiment 2 of the present invention. This embodiment refines the method described in the above embodiments for determining a target intersection map based on the target survey values and determining the target stratigraphic boundary based on the target intersection map. Figure 2 As shown, the method includes:
[0080] S210. Determine the target exploration area, obtain the historical logging data corresponding to the target exploration area, and correct the historical logging data to obtain corrected logging data, wherein the corrected logging data includes corrected phosphorus content data.
[0081] S220. Determine the boundary depth range based on the corrected phosphorus content data, and determine the target survey value corresponding to each candidate boundary depth in the boundary depth range.
[0082] S230. For each candidate boundary depth, coordinate points are determined according to the expected value and the standard deviation, and the coordinate points are marked on a two-dimensional coordinate system based on the target order to obtain the target intersection map, wherein the target order is sequentially from the candidate boundary depth with the larger depth value to the candidate boundary depth with the smaller depth value.
[0083] Here, the coordinate point can be understood as a coordinate point determined based on the expected value and the standard deviation. For example, the coordinate point is... Where A represents the coordinate point, Let σ represent the expected value and σ represent 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 map can be understood as an intersection map obtained by labeling the coordinate points corresponding to the depth of each candidate boundary on the two-dimensional coordinate system (reference). Figure 3 ).
[0086] S240. Determine the target stratigraphic boundary based on the target intersection map.
[0087] Optionally, determining the target stratigraphic boundary based on the target intersection map includes:
[0088] The coordinate points that meet the preset conditions and are initially marked are determined as target points, wherein the preset conditions are the coordinates that deviate from the origin of the two-dimensional coordinate system;
[0089] The target stratigraphic boundary is obtained by taking the candidate boundary depth corresponding to the target point as the target boundary depth.
[0090] The target point can be understood as the coordinate point that is first marked in the target intersection diagram and deviates from the origin of the two-dimensional coordinate system.
[0091] The target boundary depth can be understood as the boundary depth where the target stratum boundary is located.
[0092] The technical solution of this invention determines coordinate points for each candidate boundary depth based on the expected value and the standard deviation, and marks these coordinate points on a two-dimensional coordinate system according to the target order to obtain the target cross-plot map. The target order is sequentially from candidate boundary depths with larger depth values to candidate boundary depths with smaller depth values. This ensures the accuracy of the obtained target cross-plot map, further improving the accuracy of the target boundary depths and target stratigraphic boundaries determined based on the target cross-plot map.
[0093] Figure 4 This is an overall flowchart of a stratigraphic boundary identification method provided by an embodiment of the present invention. Figure 4 As shown, the overall process of the stratigraphic boundary identification method can be as follows:
[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 tectonic and sedimentary background of the target exploration area to understand the development characteristics of the main lithologies of the Sinian-Cambrian system. For example, the upper part of the Sinian system is mainly composed of carbonate rocks, while the lower part of the Cambrian system is mainly composed of clastic rocks. Collect well logging data from the system's drilling, and correlate the observed lithology with the lithology reflected by the natural gamma logging curves in the well logging data to establish the depth correspondence formula (1) H r =aHl +b, where H r H represents the actual depth. l Here, 'a' represents the historical depth, 'a' represents the stretching coefficient, and 'b' represents the depth correction amount. The depths of other logging curves are 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 map.
[0097] Based on the corrected phosphorus content data curve obtained from relational formula (1), determine the morphological characteristics (linear, sawtooth, step-shaped, box-shaped) and stability of the phosphorus element curve, and select the depth range (limit depth range) where the phosphorus content morphology is step-shaped or box-shaped. Then, use relational formula (2). Relationship formula (3) Calculate the expected value and standard deviation of a certain number of phosphorus content data points corresponding to each depth (candidate boundary depth). Plot the expected value and standard deviation data obtained for each depth into a two-dimensional scatter plot to obtain... Intersection diagram (target intersection diagram).
[0098] 3. Analyze the target intersection diagram to obtain the stratigraphic boundary between the Sinian System and the Cambrian System.
[0099] based on In the intersection diagram, find the target point that deviates significantly from the origin of the coordinate system from the intersection point (coordinate point), and record the corresponding candidate boundary depth, from deep to shallow. The candidate boundary depth corresponding to the first data that deviates from the origin of the coordinate system should be the Sinian-Cambrian boundary (target boundary depth).
[0100] 4. Verify the stratigraphic boundary between the Sinian and Cambrian systems by combining isotopic data and / or small shell fossil data.
[0101] By comparing the characteristics of other logging curves near the target boundary depth, there should be a significant difference between the logging curves above and below this depth. At the same time, the characteristics of small shell fossils and carbon isotopes in the vicinity of this depth can be analyzed. This depth should be the boundary between whether small shell fossils have developed and whether carbon isotope composition has changed significantly. If inconsistencies are found, it should be analyzed based on the geological background to determine whether tectonic activity has led to erosion of the strata, and whether some typical phenomena are missing.
[0102] This invention provides a method for rapidly and accurately determining the Sinian-Cambrian boundary based on phosphorus content in well logging. This method solves the problems of cumbersome and costly experimental analyses using large numbers of well cores for small-shell fossil and carbon isotope analysis, as well as the low accuracy of using seismic data for stratigraphic delineation. It also fills the gap in rapidly delineating stratigraphy using phosphorus, significantly improving work efficiency, saving costs, and providing technical and theoretical support for improving the benefits of oil and gas exploration.
[0103] Example 3
[0104] Figure 5 This is a schematic diagram of a stratigraphic boundary identification device provided in Embodiment 3 of the present invention. Figure 5 As shown, the device includes: a data correction module 310, a data calculation module 320, and a boundary recognition module 330; wherein,
[0105] The data correction module 310 is used to determine the target exploration area, acquire the historical logging data corresponding to the target exploration area, and correct the historical logging data to obtain corrected logging data, wherein the corrected logging data includes 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 map based on the target exploration value, and determine the target stratigraphic boundary based on the target intersection map, wherein the target stratigraphic boundary is the stratigraphic boundary between the Sinian System and the Cambrian System.
[0106] The technical solution of this invention involves determining a target exploration area, acquiring historical logging data corresponding to the target exploration area, correcting the historical logging data to obtain corrected logging data, wherein the corrected logging data includes corrected phosphorus content data; determining a boundary depth range based on the corrected phosphorus content data, and determining a target exploration value corresponding to each candidate boundary depth within the boundary depth range; determining a target cross-plot based on the target exploration value, and determining a target stratigraphic boundary based on the target cross-plot, wherein the target stratigraphic boundary is the stratigraphic boundary between the Sinian and Cambrian systems. This solves the problems of cumbersome and costly experimental analysis using a large number of drill cores for small-shell fossils, carbon isotopes, etc., and the low accuracy of using seismic data for stratigraphic division. It fills the gap in rapidly dividing stratigraphic layers using phosphorus. Based on the technical solution of this invention, the efficiency and accuracy of stratigraphic boundary identification are improved, and the cost of stratigraphic boundary identification is reduced.
[0107] Optionally, the historical logging data includes historical lithological data and / or historical phosphorus content data, wherein the historical lithological data includes the first historical depth of carbonate rocks and / or the second historical depth of clastic rocks, and the data correction module 310 is used for:
[0108] The actual lithological data of the target survey area are determined, wherein the actual lithological data includes the first actual depth of the carbonate rocks and / or the second actual depth of the clastic rocks, wherein the carbonate rocks are the main lithology of the upper Sinian system and the clastic rocks are the main lithology of the lower Cambrian system;
[0109] Based on the first historical depth, the second historical depth, the first actual depth, and the second actual depth, a correction coefficient is determined, 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 straight line, sawtooth, step, and box shape.
[0113] The boundary range determination unit is used to determine the depth range corresponding to the curve stages of the stepped type and the box type as the boundary depth range.
[0114] Optionally, the corrected phosphorus content data also includes adjacent phosphorus content data corresponding to surrounding survey areas 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 based on the data volume of the adjacent phosphorus content data and the data value corresponding to each of the 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: a cross-plot determination unit, used to determine coordinate points for each candidate boundary depth based on the expected value and the standard deviation, and to label the coordinate points on a two-dimensional coordinate system based on the target order to obtain the target cross-plot.
[0116] The target order is sequentially from the candidate boundary depth with the largest depth value to the candidate boundary depth with the smallest 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 points that meet the preset conditions and are initially marked as target points, wherein the preset conditions are the coordinate origins that are deviated from the two-dimensional coordinate system.
[0119] The boundary identification unit is used to take the candidate boundary depth corresponding to the target point as the target boundary depth to obtain the target stratum boundary.
[0120] Optionally, the stratigraphic 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 formation boundary within the verification depth range based on the calibration logging data. The verification data includes carbon isotope data and / or small shell fossil data. The verification depth range includes a first depth range in a first direction of the target formation boundary and / or a second depth range in a second direction of the target formation 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. The verification conditions are that the data difference between the carbon isotope data corresponding to the first depth range and the second depth range 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 exceeds a second difference threshold.
[0123] The stratigraphic boundary identification device provided in this embodiment of the invention can execute the stratigraphic boundary identification method provided in any embodiment of the invention, and has the corresponding functional modules and beneficial effects of the method.
[0124] Example 4
[0125] Figure 6 A schematic diagram of an electronic device 10 that can be used to implement embodiments 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 processors, cellular phones, smartphones, wearable devices (e.g., helmets, glasses, watches, etc.), and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely illustrative and are not intended to limit the implementation of the invention described and / or claimed herein.
[0126] like Figure 6As shown, the electronic device 10 includes at least one processor 11 and a memory, such as a read-only memory (ROM) 12 or a random access memory (RAM) 13, communicatively connected to the at least one processor 11. The memory stores computer programs executable by the at least one processor. The processor 11 can perform various appropriate actions and processes based on the computer program stored in the ROM 12 or loaded from storage unit 18 into the RAM 13. The RAM 13 may also store various programs and data required for the operation of the electronic device 10. The processor 11, ROM 12, and RAM 13 are interconnected via a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.
[0127] Multiple components in electronic device 10 are connected to I / O interface 15, including: input unit 16, such as keyboard, mouse, etc.; output unit 17, such as various types of displays, speakers, etc.; storage unit 18, such as disk, optical disk, etc.; and communication unit 19, such as network card, modem, wireless transceiver, etc. Communication unit 19 allows electronic device 10 to exchange information / data with other devices through computer networks such as the Internet and / or various telecommunications networks.
[0128] Processor 11 can be a variety of general-purpose and / or special-purpose processing components with processing and computing capabilities. Some examples of processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special-purpose artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. Processor 11 performs the various methods and processes described above, such as stratigraphic boundary identification methods.
[0129] In some embodiments, the stratigraphic boundary identification method may be implemented as a computer program tangibly contained in a computer-readable storage medium, such as storage unit 18. In some embodiments, part or all of the computer program may be loaded and / or installed on electronic device 10 via ROM 12 and / or communication unit 19. When the computer program is loaded into RAM 13 and executed by processor 11, one or more steps of the stratigraphic boundary identification method described above may be performed. Alternatively, in other embodiments, processor 11 may be configured to perform the stratigraphic boundary identification method by any other suitable means (e.g., by means of firmware).
[0130] Various embodiments 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-a-chip (SoCs), payload-programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various embodiments may include implementations in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which may be a dedicated or general-purpose programmable processor, capable of receiving data and instructions from a storage system, at least one input device, and at least one output device, and transmitting data and instructions to the storage system, the at least one input device, and the at least one output device.
[0131] Computer programs used to implement 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, such that when executed by the processor, the computer programs cause the functions / operations specified in the flowcharts and / or block diagrams to be performed. The computer programs may be executed entirely on a machine, partially on a machine, or as a standalone software package, partially on a machine and partially on a remote machine, or entirely on a remote machine or server.
[0132] In the context of this invention, a computer-readable storage medium can be a tangible medium that may contain or store a computer program for use by or in conjunction with an instruction execution system, apparatus, or device. A computer-readable storage medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, apparatus, or devices, or any suitable combination thereof. Alternatively, a computer-readable storage medium may be a machine-readable signal medium. More specific examples of machine-readable storage media include electrical connections based on one or more wires, portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fibers, portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof.
[0133] To provide interaction with a user, the systems and techniques described herein can 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 pointing device (e.g., a mouse or trackball) through which the user provides input to the electronic device. Other types of devices can also be used to provide interaction with the user; for example, feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including sound input, voice input, or tactile input).
[0134] The systems and technologies described herein can be implemented in computing systems that include backend components (e.g., as data servers), or computing systems that include middleware components (e.g., application servers), or computing systems that include frontend components (e.g., user computers with graphical user interfaces or web browsers through which users can interact with implementations of the systems and technologies described herein), or any combination of such backend, middleware, or frontend components. The components of the system can be interconnected via digital data communication of any form or medium (e.g., communication networks). Examples of communication networks include local area networks (LANs), wide area networks (WANs), blockchain networks, and the Internet.
[0135] A computing system can include clients and servers. Clients and servers are generally located far apart and typically interact through communication networks. The client-server relationship is created by computer programs running on the respective computers and having a client-server relationship with each other. The server can be a cloud server, also known as a cloud computing server or cloud host, which is a hosting product within the cloud computing service system to address the shortcomings of traditional physical hosts and VPS services, such as high management difficulty and weak business scalability.
[0136] It should be understood that the various forms of processes shown above can be used, with steps reordered, added, or deleted. For example, the steps described in this invention can be executed in parallel, sequentially, or in different orders, as long as the desired result of the technical solution of this invention can be achieved, and this is not limited herein.
[0137] The specific embodiments described above do not constitute a limitation on the scope of protection of this invention. Those skilled in the art should understand that various modifications, combinations, sub-combinations, and substitutions can be made according to design requirements and other factors. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of this invention should be included within the scope of protection of this invention.
Claims
1. A method of formation boundary identification, comprising: The method comprises the following steps: determining a target survey area, obtaining historical logging data corresponding to the target survey area, correcting the historical logging data to obtain corrected logging data, wherein the corrected logging data comprises 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; determining a target cross plot based on the target survey value, and determining a target stratigraphic boundary based on the target cross plot, wherein the target stratigraphic boundary is a stratigraphic boundary between Sinian and Cambrian systems; wherein the historical logging data comprises historical lithology data and / or historical phosphorus content data, the historical lithology data comprises a first historical depth of carbonate rock and / or a second historical depth of clastic rock, and the correction of the historical logging data to obtain corrected logging data comprises: determining actual lithology data of the target survey area, wherein the actual lithology data comprises a first actual depth of the carbonate rock and / or a second actual depth of the clastic rock, the carbonate rock being the main lithology at the top of the Sinian system, and the clastic rock being the main lithology at the bottom of the Cambrian system; determining 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 comprises a stretching coefficient and / or a depth correction amount; correcting the historical logging data based on the correction coefficient to obtain corrected logging data; the target survey value comprises an expected value and / or a standard deviation, and the determination of the target cross plot based on the target survey value comprises: for each candidate limit depth, determining a coordinate point according to the expected value and the standard deviation, and labeling the coordinate point on a two-dimensional coordinate system based on a target order to obtain the target cross plot, wherein the target order is along the candidate limit depths from a candidate limit depth with a large depth value to a candidate limit depth with a small depth value.
2. The method of claim 1, wherein, The corrected phosphorus content data comprises a depth phosphorus content curve corresponding to the target survey area, and the determination of the limit depth range based on the corrected phosphorus content data comprises: determining a morphological feature of the depth phosphorus content curve, wherein the morphological feature comprises at least one of a straight line type, a sawtooth shape, a step type, and a box type; determining a depth range corresponding to a curve stage of the step type and the box type as the limit depth range.
3. The method of claim 1, wherein, The corrected phosphorus content data further comprises adjacent phosphorus content data corresponding to a peripheral survey area adjacent to the target survey area, and the determination of the target survey value corresponding to each candidate limit depth in the limit depth range comprises: for each candidate limit depth in the limit depth range, determining a target survey value according to a data amount of the adjacent phosphorus content data and a data value corresponding to each adjacent phosphorus content data.
4. The method of claim 1, wherein, The determination of the target stratigraphic boundary based on the target cross plot comprises: determining a first labeled coordinate point meeting a preset condition as a target point, wherein the preset condition is deviating from a coordinate origin of the two-dimensional coordinate system. The candidate limit depth corresponding to the target point is taken as a target limit depth, and the target stratum limit is obtained.
5. The method of claim 1, wherein, Further comprising: obtaining verification data of the target stratum limit in 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 stratum limit and / or a second depth range in a second direction of the target stratum limit, the first direction and the second direction being perpendicular to a ground surface, the first direction being opposite to the second direction; in a case where the verification data meets a verification condition, determining the target stratum limit as an actual stratum limit, the verification condition being that a 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 a 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.
6. A formation boundary identification apparatus characterized by, Further comprising: a data correction module configured to determine a target survey area, obtain historical logging data corresponding to the target survey area, and correct the historical logging data to obtain corrected logging data, wherein the corrected logging data includes corrected phosphorus content data; a data calculation module configured to determine a limit depth range based on the corrected phosphorus content data, and determine a target survey value corresponding to each candidate limit depth in the limit depth range; a limit identification module configured to determine a target cross plot based on the target survey value, and determine a target stratum limit based on the target cross plot, wherein the target stratum limit is a stratum limit between Sinian and Cambrian systems; 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, the data correction module is specifically configured to 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 a main lithology at the top of the Sinian system, and the clastic rock is a main lithology at the bottom of the Cambrian system; determine a correction coefficient based on 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; correct the historical logging data based on the correction coefficient to obtain the corrected logging data; the target survey value includes an expected value and / or a standard deviation, and the limit identification module includes: a cross plot determination unit configured to, for each candidate limit depth, determine a coordinate point according to the expected value and the standard deviation, and label the coordinate point on a two-dimensional coordinate system based on a target order to obtain the target cross plot, wherein the target order is a sequence along the candidate limit depths from a candidate limit depth with a large depth value to a candidate limit depth with a small depth value.
7. An electronic device, comprising: The electronic device includes: at least one processor; and a memory connected with the at least one processor in communication; 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 to enable the at least one processor to perform the method for identifying formation boundaries according to any one of claims 1-5.
8. A computer-readable storage medium, characterized in that, The computer readable storage medium stores computer instructions for causing the processor to implement the method for identifying formation boundaries according to any one of claims 1-5 when executed.