A method, device, electronic device and storage medium for restoring basin boundaries

Abstract

The present invention relates to a method, device, electronic device and storage medium for restoring a basin boundary, and belongs to the field of geological technology, wherein the method comprises: obtaining the distribution information of the current stratum interface; according to the distribution information of the current stratum interface, using a preset balanced profile method to obtain the initial state of stratum deposition; constructing the correlation relationship between the stratum thickness of each point in the initial state of stratum deposition and the distribution information of the current stratum interface; based on the correlation relationship, splicing multiple continuous deposition areas across faults to obtain a target stratum thickness and stratum interface correlation model; according to the target stratum thickness and stratum interface correlation model, determining the basin boundary. The present invention solves the technical problem that the basin boundary and scale are difficult to accurately determine in the absence of well data.

Description

Technical Field

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

[0002] A prototype basin, also known as a basin prototype, refers to a basin with a specific sedimentary entity formed in a relatively single geodynamic system or a single-cycle tectonic evolution stage. In other words, in the entire geological history evolution process, a tectonic layer formed at different evolutionary stages, and corresponding to the evolutionary stage of each specific tectonic layer, there is a relatively stable geotectonic environment and a dominant sedimentation mechanism and sedimentary filling combination, as well as a definite basin boundary. Due to the incomplete preservation of strata development caused by tectonic activities, it is difficult to determine the boundaries and scale of the basin, which affects the analysis and understanding of the basin, especially the understanding of the structure and sediment distribution inside the basin, and thus affects the exploration of resources. In order to accurately predict the distribution of sedimentary bodies in the basin, it is necessary to analyze the basin filling process and restore the original state or basin boundary of the basin at the time of basin deposition based on the remaining strata. At present, the restoration of basin boundaries and basin scales is mainly based on stratigraphic methods, basin margin phases, clastic rock mineral maturity, clastic heavy mineral stability coefficient and Pan-Kriging method. The advantages of various methods are used for comprehensive analysis to restore the prototype basin boundaries as accurately as possible.

[0003] However, these methods still have certain limitations, because complex geological processes such as tectonic activities and sedimentation during the evolution of basins may cause dislocation and faulting of strata, making it difficult to accurately determine the boundaries and scale of prototype basins in the absence of well data. Among them, well data refers to data on underground rock structure, properties and stratigraphic information obtained through geological exploration work around the wellhead. These data include core samples, logging data, rock sample analysis results, wellbore stratigraphic records, ground temperature, ground stress and other observation data during the well drilling process. Well data is an important source of data for studying geological structure, oil and gas exploration, groundwater resource evaluation and other fields, and can provide detailed geological information and parameter references for geological science research and resource exploration. Summary of the invention

[0004] In view of this, it is necessary to provide a basin boundary recovery method, device, electronic device and storage medium to solve the technical problem in the prior art that the basin boundary and scale are difficult to accurately determine in the absence of well data.

[0005] In order to solve the above problems, the present invention provides a method for restoring a basin boundary, comprising:

[0006] Get the distribution information of the current stratum interface;

[0007] According to the distribution information of the current stratigraphic interface, the initial state of stratigraphic deposition is obtained by using the preset equilibrium profile method;

[0008] Constructing the correlation between the stratum thickness at each point in the initial state of stratum deposition and the distribution information of the current stratum interface;

[0009] Based on the correlation, multiple continuous sedimentary areas across the faults are spliced ​​to obtain a correlation model between target stratum thickness and stratum interface;

[0010] The basin boundary is determined according to the target stratum thickness and stratum interface correlation model.

[0011] In a possible implementation, the method of obtaining the initial state of stratum deposition by using a preset equilibrium profile method according to the current distribution information of stratum interfaces includes:

[0012] Establish a comprehensive velocity field in the target area, and based on the comprehensive velocity field, convert the time domain profile into a depth domain profile to obtain the stratigraphic structure profile;

[0013] Based on the formation structure profile, the preset analysis method is used to decompact the formation to be restored, and the formation to be restored in a balanced state is obtained;

[0014] The preset analysis method is used to eliminate the influence of faults on sedimentation, and the initial state of stratigraphic sedimentation of the stratum to be restored in a balanced state is obtained.

[0015] In a possible implementation, the constructing of the correlation between the stratum thickness at each point in the initial state of stratum deposition and the distribution information of the current stratum interface includes:

[0016] According to the stratigraphic structural section, thickness information at different positions of the stratigraphic formation to be restored is obtained;

[0017] A plane rectangular coordinate system is constructed by taking the position where the formation thickness is 0 as the coordinate origin, the thickness information at different positions as the horizontal coordinate, and the distribution information of the current formation interface as the vertical coordinate;

[0018] Based on the plane rectangular coordinate system, a functional correlation relationship between the formation thickness and the distribution information of the current formation interface is established.

[0019] In a possible implementation, the formation thickness can be expressed by the following formula:

[0020]

[0021] Among them, H is the formation thickness of the target point, L is the distance between the target point and the basin boundary, and a, b and c are all formation thickness adjustment coefficients.

[0022] In a possible implementation, the formation thickness can also be expressed by the following formula:

[0023]

[0024] Among them, C 顶 is the top depth of the target stratum; C 底 is the bottom depth of the target stratum, and a, b and c are the adjustment coefficients of stratum thickness.

[0025] In a possible implementation, after splicing a plurality of continuous sedimentary areas across the fault to obtain a correlation model between target stratum thickness and stratum interface, the method further includes:

[0026] Obtain the original position of the fault;

[0027] According to the correlation model between the target stratum thickness and the stratum interface, the splicing position of the fault is obtained;

[0028] The length of the fault is determined based on the difference between the splicing position of the fault and the original position.

[0029] In a possible implementation, the distribution information of the current stratum interface includes a stratum boundary, a top boundary, and a bottom boundary.

[0030] In a second aspect, the present invention further provides a basin boundary restoration device, comprising:

[0031] An acquisition module is used to obtain the distribution information of the current stratum interface;

[0032] The initial deposition acquisition module is used to obtain the initial state of the stratum deposition using a preset equilibrium profile method according to the distribution information of the current stratum interface;

[0033] A relationship building module is used to build the relationship between the stratum thickness at each point in the initial state of stratum deposition and the distribution information of the current stratum interface;

[0034] A target model determination module is used to splice a plurality of continuous sedimentary areas across the fault based on the association relationship to obtain a target formation thickness and formation interface association model;

[0035] The basin boundary determination module is used to determine the basin boundary according to the target stratum thickness and stratum interface association model.

[0036] In a third aspect, the present invention further provides an electronic device, comprising: a processor and a memory;

[0037] The memory stores a computer-readable program executable by the processor;

[0038] When the processor executes the computer-readable program, the processor implements the steps in the basin boundary restoration method described above.

[0039] In a fourth aspect, the present invention further provides a computer-readable storage medium storing one or more programs, wherein the one or more programs can be executed by one or more processors to implement the steps in the basin boundary restoration method as described above.

[0040] The beneficial effects of the present invention are: firstly, the distribution information of the current formation interface is obtained, and then the initial deposition state of the formation is restored by using the balanced profile method, and then the sediments in the basin change regularly from the edge of the basin to the center of the basin, so as to construct the correlation between the formation thickness of each point in the initial state of formation deposition and the distribution information of the current formation interface, and obtain the target formation thickness and formation interface correlation model by splicing multiple continuous deposition areas across faults; finally, the basin boundary is determined according to the target formation thickness and formation interface correlation model. The purpose of accurately restoring the formation boundary and scale in the absence of well data is achieved. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] FIG1 is a method flow chart of an embodiment of a method for restoring a basin boundary provided by the present invention;

[0042] Figure 2 A schematic diagram of a target seismic profile selected in the basin boundary recovery method provided by the present invention;

[0043] Figure 3 A method flow chart of an embodiment of step S102 in the basin boundary restoration method provided by the present invention;

[0044] Figure 4 A schematic diagram of a stratigraphic structural cross section in the basin boundary restoration method provided by the present invention;

[0045] Figure 5 A flow chart of the inversion simulation technology in the basin boundary restoration method provided by the present invention;

[0046] Figure 6 A schematic diagram of a correlation model between target stratum thickness and stratum interface in a basin boundary restoration method provided by the present invention;

[0047] Figure 7 A method flow chart of an embodiment of step S103 in the basin boundary restoration method provided by the present invention;

[0048] Figure 8 is a schematic diagram of an embodiment of a basin boundary restoration device provided by the present invention;

[0049] Fig. 9It is a schematic diagram of the operating environment of an electronic device provided by the present invention. DETAILED DESCRIPTION

[0050] The preferred embodiments of the present invention are described in detail below in conjunction with the accompanying drawings, wherein the accompanying drawings constitute a part of this application and are used together with the embodiments of the present invention to illustrate the principles of the present invention, but are not used to limit the scope of the present invention.

[0051] A specific embodiment of the present invention discloses a method for restoring basin boundaries. Figure 1 ,include:

[0052] S101, obtaining distribution information of current stratum interface;

[0053] S102, according to the distribution information of the current stratum interface, using the preset equilibrium profile method to obtain the initial state of stratum deposition;

[0054] S103, constructing a correlation relationship between the stratum thickness at each point in the initial state of stratum deposition and the distribution information of the current stratum interface;

[0055] S104, based on the association relationship, splicing multiple continuous sedimentary areas across the fault to obtain a target formation thickness and formation interface association model;

[0056] S105. Determine the basin boundary according to the target stratum thickness and stratum interface correlation model.

[0057] In this embodiment, the distribution information of the current formation interface is first obtained, and then the initial deposition state of the formation is restored by using the balanced profile method. Then, according to the regular change of sediments in the basin from the edge of the basin to the center of the basin, the relationship between the formation thickness of each point in the initial state of formation deposition and the distribution information of the current formation interface is constructed, and the target formation thickness and formation interface association model is obtained by splicing multiple continuous deposition areas across faults; finally, the basin boundary is determined according to the target formation thickness and formation interface association model. The purpose of accurately restoring the formation boundary and scale in the absence of well data is achieved.

[0058] It should be noted that the inversion effect of the balanced profile method is used, that is, on the unit survey line, the balanced profile method is applied to restore its initial sedimentary state, and then restore the prototype basin boundary and scale. The balanced profile method is a commonly used sedimentary basin analysis method. By analyzing the thickness, dip, sequence and other information of the strata, the structural evolution during the basin deposition process and the changes in the sedimentary basin boundary are inverted.

[0059] Among them, when applying the balanced profile method to restore the basin boundary, first select a suitable survey line and obtain sufficient geological data, including stratigraphic structure, structural characteristics, lithology and other information. Then, the initial sedimentary state on the survey line is calculated and restored by the balanced profile method, that is, the original distribution of the strata in the absence of structural deformation. By comparing the restored original sedimentary state with the actual observed stratigraphic distribution, the changes in the basin boundary position, sedimentary system and structural characteristics can be identified. By using the balanced profile method to restore the basin boundary, we can better understand the evolutionary history and structural characteristics of the basin, and provide important geological background and reference information for oil and gas exploration, geological disaster risk assessment and other work. At the same time, combined with seismic, electromagnetic and other exploration technologies, the restoration results of the basin boundary can be further verified and improved to improve its reliability and accuracy. This comprehensive application method can more comprehensively reveal the geological characteristics of the basin and provide strong support for geological research and resource exploration.

[0060] In step S101, refer to Figure 2 In the three-dimensional seismic data volume, a section with a relatively clear stratigraphic distribution is selected based on the seismic section, the phase axis of the sedimentary interface is identified, the strata are interpreted, the positions of the top and bottom interfaces of the strata are found, and the distribution of the current stratigraphic interface is obtained.

[0061] In some embodiments, the initial state of stratum deposition is obtained by using a preset equilibrium profile method according to the distribution information of the current stratum interface. Figure 3 ,include:

[0062] S301, establishing a comprehensive velocity field in the target area, and based on the comprehensive velocity field, converting the time domain profile into a depth domain profile to obtain a stratigraphic structural profile;

[0063] S302, based on the formation structure profile, using a preset analysis method to decompact the formation to be restored, to obtain the formation to be restored in a balanced state;

[0064] S303, using a preset analysis method to eliminate the influence of the fault on the sedimentation, and obtaining the initial state of the sedimentation of the stratum to be restored in a balanced state.

[0065] In step S301, establishing a comprehensive velocity field means conducting seismic exploration in the target area, and establishing a comprehensive velocity field by measuring the velocity distribution of different strata, which includes the velocity of seismic waves propagating in different strata. Converting the selected time domain profile into a depth domain profile means converting the time information in the seismic data into depth information. By referring to the velocity field, the time domain data can be accurately converted into the depth domain data. By obtaining the structural profile, such as Figure 4As shown, the structural characteristics of the underground structure in the target area can be more intuitively displayed to help determine the nature of the fault in the target area. By analyzing the shape and changes of the structural section, the strike, dip, slip direction and other information of the fault can be inferred.

[0066] In both step S302 and step S303, 2D move software is used to remove compaction or eliminate the influence of faults. It is understandable that other processing methods may also be used.

[0067] Furthermore, in step S302, the balanced profile technology is divided into forward method and inversion method. The inversion method is to restore a real, actual deformed profile to its original, undeformed state under the principle of relevant balance. The inversion method is used to restore the balanced profile, and its basic process is as follows: Figure 5 As shown. Decompaction correction is actually the restoration correction of stratum porosity, and porosity is mainly controlled by stratum burial depth and gravity load of overlying strata. In order to accurately reflect the ancient structural morphology, the rock strata must be decompacted before the profile is balanced; after stripping off the upper sedimentary strata layer by layer, the true thickness of the lower strata at the end of deposition is restored. If the result is unbalanced, the time-depth conversion map and the restoration map must be recalibrated. It has been verified that the restoration map after decompaction is balanced.

[0068] In step S303, the influence of faults on sedimentation is eliminated, and each layer is flattened to obtain the initial sedimentation state of the stratum and determine the sedimentation center. Figure 6 The image shown is the image after flattening.

[0069] In some embodiments, the relationship between the thickness of each point at the initial state of the constructed stratum deposition and the distribution information of the current stratum interface can be found in Figure 7 ,include:

[0070] S701, obtaining thickness information of different positions of the stratum to be restored according to the stratum structural section;

[0071] S702, constructing a plane rectangular coordinate system with the position where the formation thickness is 0 as the coordinate origin, the thickness information at different positions as the abscissa, and the distribution information of the current formation interface as the ordinate;

[0072] S703: Based on a plane rectangular coordinate system, a functional correlation relationship between the stratum thickness and the distribution information of the current stratum interface is established.

[0073] In this embodiment, by establishing a functional relationship to characterize the relationship between the stratum thickness and the distribution information of the current stratum interface, the sedimentary process and structural evolution can be better understood.

[0074] Specifically, based on the continuous sedimentary stratigraphic segments, the software is used to establish the equation:

[0075]

[0076] Where: H is the formation thickness of the target point, in meters; L is the distance between the target point and the basin boundary, in kilometers; a, b, and c are all formation thickness adjustment coefficients. This equation can be used to calculate the original formation thickness at the target point.

[0077] Furthermore, the thickness of the denuded strata can be analyzed by collecting data to establish a regression equation. First, the top and bottom depth data of the un-eroded area of ​​the target stratum of the survey line are read and statistically analyzed to ensure that the sedimentation model applicable to the sedimentation simulation method is established. Regression analysis is performed on the top and bottom depth data of the un-eroded area to establish a regression equation:

[0078]

[0079] In the formula: Ctop is the top depth of the target stratum; Cbottom is the bottom depth of the target stratum, and a, b and c are all stratum thickness adjustment coefficients. Using this equation, the top boundary of the original stratum in the unconformity area can be obtained, and the difference between the top boundary of the original stratum and the current unconformity surface depth can be used to calculate the erosion amount of the unconformity surface.

[0080] In some embodiments, after splicing a plurality of continuous sedimentary areas across the faults to obtain a correlation model between target stratum thickness and stratum interface, the method further includes:

[0081] Obtain the original position of the fault;

[0082] According to the correlation model between the target stratum thickness and the stratum interface, the splicing position of the fault is obtained;

[0083] The length of the fault is determined based on the difference between the splicing position of the fault and the original position.

[0084] In this embodiment, the coordinates of other continuous sedimentary areas affected by the fault are fitted with the correlation model, and the original horizontal coordinates at the breakpoint are marked as , the horizontal coordinate after the fitting equation is shifted is marked as , the difference between the two coordinates , which is the length of the missing deposition area.

[0085] Furthermore, when calculating the initial sedimentary boundary of the target horizon in the area controlled by the tensional fault, the original horizontal coordinate at the breakpoint is marked as , the horizontal coordinate after the fitting equation is shifted is marked as , the difference between the two coordinates This is the distance the basin expanded during its evolution. Subtracting this distance from the current residual boundary gives the corrected residual boundary length.

[0086] Based on the above basin boundary restoration method, the present invention also provides a basin boundary restoration device, see Figure 8 , including: an acquisition module 810, an initial deposition acquisition module 820, a relationship construction module 830, a target model determination module 840 and a basin boundary determination module 850.

[0087] An acquisition module 810 is used to acquire the distribution information of the current stratum interface;

[0088] The initial deposition acquisition module 820 is used to obtain the initial state of the formation deposition using a preset equilibrium profile method according to the distribution information of the current formation interface;

[0089] A relationship building module 830 is used to build a correlation relationship between the stratum thickness at each point in the initial state of stratum deposition and the distribution information of the current stratum interface;

[0090] A target model determination module 840 is used to splice a plurality of continuous sedimentary areas across the fault based on the association relationship to obtain a target formation thickness and formation interface association model;

[0091] The basin boundary determination module 850 is used to determine the basin boundary according to the target formation thickness and formation interface association model.

[0092] like Fig. 9 As shown, based on the above basin boundary restoration method, the present invention also provides an electronic device, which can be a computing electronic device such as a mobile terminal, a desktop computer, a notebook, a palm computer, and a server. The electronic device includes a processor 910, a memory 920, and a display 930. Fig. 9 Only some components of the electronic device are shown, but it should be understood that it is not required to implement all of the components shown, and more or fewer components may be implemented instead.

[0093] In some embodiments, the memory 920 may be an internal storage unit of the electronic device, such as a hard disk or memory of the electronic device. In other embodiments, the memory 920 may also be an external storage electronic device of the electronic device, such as a plug-in hard disk, a smart memory card (Smart Media Card, SMC), a secure digital (SecureDigital, SD) card, a flash card (Flash Card), etc. equipped on the electronic device. Further, the memory 920 may also include both an internal storage unit of the electronic device and an external storage electronic device. The memory 920 is used to store application software and various types of data installed in the electronic device, such as program codes installed in the electronic device. The memory 920 may also be used to temporarily store data that has been output or is to be output. In one embodiment, a basin boundary recovery program 940 is stored on the memory 920, and the basin boundary recovery program 940 can be executed by the processor 910, thereby realizing the basin boundary recovery method of each embodiment of the present application.

[0094] In some embodiments, the processor 910 may be a central processing unit (CPU), a microprocessor or other data processing chip, and is used to run program codes or process data stored in the memory 920, such as executing a basin boundary recovery method.

[0095] In some embodiments, the display 930 may be an LED display, a liquid crystal display, a touch-sensitive liquid crystal display, an OLED (Organic Light-Emitting Diode) touch device, etc. The display 930 is used to display information of the recovery electronic device at the basin boundary and to display a visual user interface. The components 910-930 of the electronic device communicate with each other through a system bus.

[0096] Those skilled in the art will appreciate that all or part of the processes of the above-mentioned embodiments can be implemented by instructing related hardware through a computer program, and the program can be stored in a computer-readable storage medium, wherein the computer-readable storage medium is a disk, an optical disk, a read-only storage memory, or a random access memory, etc.

[0097] The above description is only a preferred specific implementation manner of the present invention, but the protection scope of the present invention is not limited thereto. Any changes or substitutions that can be easily conceived by any technician familiar with the technical field within the technical scope disclosed by the present invention should be covered within the protection scope of the present invention.

Claims

1. A method for restoring a basin boundary, characterized in that: include: Get the distribution information of the current stratum interface; A comprehensive velocity field is established in the target area, and based on the comprehensive velocity field, the time domain profile is converted into a depth domain profile to obtain a stratigraphic structural profile; based on the stratigraphic structural profile, a preset analysis method is used to decompact the stratum to be restored to obtain a stratum to be restored in a balanced state; a preset analysis method is used to eliminate the influence of faults on sedimentation to obtain the initial state of stratigraphic sedimentation of the stratum to be restored in a balanced state; According to the stratigraphic structural section, thickness information at different positions of the stratigraphic formation to be restored is obtained; A plane rectangular coordinate system is constructed with the position where the formation thickness is 0 as the coordinate origin, the thickness information at different positions as the abscissa, and the distribution information of the current formation interface as the ordinate. Based on the plane rectangular coordinate system, a functional correlation relationship between the formation thickness and the distribution information of the current formation interface is established; wherein the formation thickness is expressed by the following formula: , where H is the formation thickness of the target point, L is the distance between the target point and the basin boundary, and a, b and c are all formation thickness adjustment coefficients; Based on the correlation, multiple continuous sedimentary areas across the faults are spliced ​​to obtain a correlation model between target stratum thickness and stratum interface; The basin boundary is determined according to the target stratum thickness and stratum interface correlation model.

2. The method for restoring the basin boundary according to claim 1, characterized in that: After splicing multiple continuous sedimentary areas across faults and obtaining the correlation model between target stratum thickness and stratum interface, it also includes: Obtain the original position of the fault; According to the correlation model between the target stratum thickness and the stratum interface, the splicing position of the fault is obtained; The length of the fault is determined based on the difference between the splicing position of the fault and the original position.

3. The method for restoring the basin boundary according to claim 2, characterized in that: The distribution information of the current stratum interface includes stratum boundary, top layer boundary and bottom layer boundary.

4. A basin boundary restoration device, characterized in that: include: An acquisition module is used to obtain the distribution information of the current stratum interface; The initial deposition acquisition module is used to establish a comprehensive velocity field in the target area, and based on the comprehensive velocity field, convert the time domain profile into a depth domain profile to obtain a stratigraphic structural profile; based on the stratigraphic structural profile, a preset analysis method is used to decompact the stratum to be restored to obtain a stratum to be restored in a balanced state; a preset analysis method is used to eliminate the influence of faults on deposition to obtain the initial state of stratigraphic deposition of the stratum to be restored in a balanced state; A relationship building module, used to obtain thickness information at different positions of the stratum to be restored according to the stratum structural section; A plane rectangular coordinate system is constructed with the position where the formation thickness is 0 as the coordinate origin, the thickness information at different positions as the abscissa, and the distribution information of the current formation interface as the ordinate. Based on the plane rectangular coordinate system, a functional correlation relationship between the formation thickness and the distribution information of the current formation interface is established; wherein the formation thickness is expressed by the following formula: , where H is the formation thickness of the target point, L is the distance between the target point and the basin boundary, and a, b and c are all formation thickness adjustment coefficients; A target model determination module is used to splice a plurality of continuous sedimentary areas across the fault based on the association relationship to obtain a target formation thickness and formation interface association model; The basin boundary determination module is used to determine the basin boundary according to the target stratum thickness and stratum interface association model.

5. An electronic device, characterized in that: include: Processor and memory; The memory stores a computer-readable program executable by the processor; When the processor executes the computer-readable program, the processor implements the steps in the basin boundary restoration method according to claims 1-3.

6. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores one or more programs, and the one or more programs can be executed by one or more processors to implement the steps in the basin boundary restoration method according to claims 1-3.

Images