Method and device for predicting thickness of salt layer in salt structure recovery, equipment and medium

Abstract

The invention relates to a method and device for predicting the thickness of a salt layer in salt structure recovery, equipment and a medium, and relates to the field of structure evolution of a salt-containing basin, and the method comprises the steps: calculating the average change rate S of the salt layer according to the initial salt layer thickness and the current salt layer thickness of a salt structure in a research area; defining special events in the salt basin evolution process based on the basin evolution background, and obtaining change coefficients of the salt layer thickness under different special events; and based on the average change rate of the salt layer, predicting the salt layer thickness corresponding to the special event period by taking the change coefficient as a constraint. According to the prediction method provided by the invention, the influence of special events in the salt structure is considered, the salt bed thickness is restrained through geologic data such as geodynamics background and lithologic change aiming at the thickness change caused by nonconservation and viscoplastic deformation of salt rock in the salt structure recovery process, and the evaluation of the salt bed thickness in different periods is realized; and efficient prediction of the salt bed thickness change during salt structure recovery is realized.

Description

Technical Field

[0001] The present invention relates to restoration of structural evolution of salt-bearing basins, and in particular to a method, device, equipment and medium for predicting salt layer thickness in salt structure restoration. Background Art

[0002] Currently, traditional profile restoration emphasizes material conservation, enabling quantitative and semi-quantitative analysis of stratigraphic evolution. However, the fluidity of salt rock causes the salt rock area in current profiles to differ from that in different periods. Although it can be argued that plastic flow of salt rock can "flow in" or "out" of a two-dimensional profile, in the direction of regional transport (principal stress), salt bodies generally deform within the two-dimensional plane, resulting in minimal changes in salt rock area. Therefore, this is considered to conform to the principle of "two-dimensional area balance" to a certain extent.

[0003] However, this assumption fails to accurately reflect the evolution of salt structures, particularly the thickness variations of salt rocks. Furthermore, some salt structure reconstruction methods emphasize stratigraphic changes above the salt, estimating the initial length of the salt layer based on the deformation of the sedimentary layer. However, this method also fails to accurately determine the thickness variations of the salt layer. Therefore, calculating the thickness variations of the salt layer during salt structure profile reconstruction is a key issue currently facing researchers.

[0004] Profile restoration technology is a core technology in structural geology. It can intuitively reflect the evolution and details of a basin, the development of different structural styles, and the differential evolution under different stresses. It can also summarize the laws of fault activity, clarify the oil and gas migration pathways in the basin, and provide advantageous reservoirs. In saline basins, balanced profile technology is equally important and can provide technical support in many aspects, including providing a complete salt structure development process, quantifying the impact of structure and sedimentation on salt structures, demarcating key periods and nodes in salt basin development, and providing prototypes for other studies. More importantly, salt structure profile restoration can provide relevant researchers or oil and gas practitioners in non-salt structural fields with a rapid understanding of regional salt basins, and accurate assessment of salt layer thickness can lay the foundation for the next step of oil and gas resource exploration and development.

[0005] However, due to the unique rheological properties and complexity of structural styles of salt structures, salt-related structures cannot be restored by simply removing fault distances and flattening layers, which is different from other types of profile restoration. This is because basin dynamics, sedimentary processes, and evolutionary patterns of different salt structural styles in different periods have a certain impact on salt-related structures. Therefore, there is an urgent need for a relatively quantitative and standardized method to control the changes in salt layer thickness during salt structure restoration, so as to improve the prediction accuracy of changes in salt layer thickness during salt structure restoration. Summary of the Invention

[0006] In view of the problems existing in the prior art, the purpose of the present invention is to provide a method and device, equipment, and medium for predicting salt layer thickness in salt structure restoration, so as to solve the defect of low accuracy in predicting changes in salt layer thickness during salt-related structure restoration due to basin dynamics, sedimentary processes, and evolutionary patterns of different salt structure styles in different periods.

[0007] To achieve this object, the present invention adopts the following technical solutions:

[0008] In a first aspect, the present invention provides a method for predicting salt layer thickness in salt structure restoration, the prediction method comprising:

[0009] Based on the initial salt layer thickness and the current salt layer thickness of the salt structure in the study area, the average salt layer change rate S is calculated;

[0010] Based on the basin evolution background, special events in the salt basin evolution process are defined, and the variation coefficient of salt layer thickness under different special events is obtained;

[0011] Based on the average change rate of the salt layer and the change coefficient as a constraint, the salt layer thickness corresponding to the special event period is predicted.

[0012] The prediction method provided by the present invention considers the impact of special events in salt structures, targets the thickness changes caused by the non-conservation and viscoplastic deformation of salt rocks during the restoration of salt structures, constrains the thickness of salt layers through geological data such as geodynamic background and lithologic changes, realizes the evaluation of salt layer thickness at different periods, and achieves efficient prediction of salt layer thickness changes during salt structure restoration.

[0013] As a preferred technical solution of the present invention, the calculation formula of the average change rate S of the salt layer is as follows:

[0014] S=|Nn| / T

[0015] Where N is the initial salt layer thickness, m; n is the current salt layer thickness, m; and T is the evolution time to date, myr.

[0016] As a preferred technical solution of the present invention, the definition of special events in the evolution of salt basins based on the basin evolution background includes: establishing a complete basin dynamics framework and sedimentary stratigraphic model based on the geophysical data of the study area, and determining the key influencing events of salt basin development in different periods.

[0017] Preferably, the geophysical data includes: well logging data and seismic data.

[0018] As a preferred technical solution of the present invention, the variation coefficient includes: an initial coefficient a, a heat sink coefficient b, or a welding coefficient c, or a combination of at least two of them.

[0019] As a preferred technical solution of the present invention, the initial coefficient a includes: initial clastic rock layer coefficient a 碎 Crumbs.

[0020] Preferably, the initial clastic rock layer coefficient a 碎屑 =β×S, where, if the thickness of the clastic rock is less than the initial salt layer thickness N, β=k×clastic rock thickness, k=0.00025, β<0.5; if the thickness of the clastic rock is greater than or equal to the initial salt layer thickness N, β is equal to 0.5.

[0021] As a preferred technical solution of the present invention, the heat sink coefficient b includes: based on the substrate lifting angle during the heat sink period, when the substrate lifting angle is 2°, b is 3; based on the substrate lifting angle of 2°, b increases by 2 for every 1° increase in the substrate lifting angle; when the substrate lifting angle is greater than 4°, b is 7; when the substrate lifting angle is less than 2°, b decreases by 1 for every 1° decrease in the substrate lifting angle.

[0022] Preferably, the welding coefficient c includes: based on the welding residual thickness, if the welding residual thickness is ≤20m, then c is 0.01; if the welding residual thickness is >20m, then c=0.05.

[0023] As a preferred technical solution of the present invention, the relationship used to predict the salt layer thickness corresponding to the special event period includes:

[0024] S 初期 ×t1+S 热沉降 ×t2+S 焊接后 ×t3+S * ×t n =H

[0025] In the formula, when the rock layer is a carbonate rock layer, S 初期 =0.1×(h 初期碳酸 / 400m)×H / t 初期碳酸 , H is the total salt layer thickness change, m; h 初期碳酸 is the thickness of the carbonate formation determined by actual core sampling or seismic profile, m; t 初期碳酸 It is determined by actual core sampling or seismic profile, which is obtained by subtracting the age of the bottom from the age of the top of the carbonate rock. If the rock layer is a clastic rock layer, then S 初期 = initial coefficient a × average rate of change S; S 热沉降 = heat sink coefficient b×average rate of change S; S 焊接后 = welding coefficient c×average rate of change S; S * = other coefficients of change × average rate of change S; t1 + t2 + t3 + ... t n=T, t1 is the sedimentation time corresponding to the post-salt carbonate stratum or the sedimentation time corresponding to the thickness of the clastic sedimentary rock reaching the critical point of density inversion, t2 is the duration of basin thermal subsidence, t3 is the duration from the formation of salt welding to the present, t n is the duration corresponding to other variation coefficients.

[0026] In a second aspect, the present invention provides a device for predicting salt layer thickness in salt structure restoration, the device comprising:

[0027] The change rate acquisition module is used to calculate the average change rate S of the salt layer based on the initial salt layer thickness and the current salt layer thickness of the salt structure in the study area;

[0028] The variation coefficient acquisition module is used to define special events in the evolution of salt basins based on the basin evolution background and obtain the variation coefficient of salt layer thickness under different special events;

[0029] The prediction module is used to predict the salt layer thickness corresponding to the special event period based on the average change rate of the salt layer and the change coefficient as a constraint.

[0030] In a third aspect, the present invention provides an electronic device, comprising:

[0031] at least one processor; and a memory communicatively coupled to the at least one processor;

[0032] Wherein, the memory stores a computer program that can be executed by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the method for predicting salt layer thickness in salt structure restoration described in the first aspect.

[0033] In a fourth aspect, an embodiment of the present invention provides a computer storage medium storing computer executable instructions, which, when executed by a processor, implements the method for predicting salt layer thickness in salt structure restoration described in the first aspect.

[0034] Compared with the existing technical solutions, the present invention has the following beneficial effects:

[0035] (1) The prediction method provided by the present invention establishes a complete basin dynamics framework based on the acquisition of geophysical data such as well logging data and seismic data in the study area. The key special events in the development of the salt basin at different stages are determined, and each event is set to control the main changes in the thickness of the salt layer at that stage. There is no specific limit on the number of special events, which are usually set according to the basin dynamics background and sedimentary lithology changes. Each special event corresponds to a period of time, but setting more special events can establish a more accurate change process for the salt layer thickness. For time periods without obvious special events, it can be judged based on the accumulation rate of sediments, or the salt layer thickness change rate during that period can be directly equal to the overall average rate.

[0036] (2) The prediction method provided by the present invention targets the thickness changes caused by the non-conservation and viscoplastic deformation of salt rocks during the restoration of salt structures. The thickness of the salt layer is constrained by geological data such as the geodynamic background and lithologic changes, thereby realizing the evaluation of the thickness of the salt layer at different periods. The prediction process can objectively reduce the influence of subjective consciousness, thereby obtaining a more realistic profile restoration result. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] Figure 1 This is a flow chart of a method for predicting salt layer thickness in salt structure restoration provided by an embodiment of the present invention;

[0038] Figure 2 Schematic diagram of a device for predicting salt layer thickness in salt structure restoration provided by an embodiment of the present invention;

[0039] Figure 3 is a schematic diagram of an electronic device provided by an embodiment of the present invention;

[0040] Figure 4 Schematic diagram of the salt layer thickness calculation process in Example 1 of the present invention;

[0041] Figure 5 This is an example of salt structure restoration and thickness change pattern diagram in Example 1 of the present invention.

[0042] In the figure: 10 - electronic device, 11 - processor, 12 - ROM, 13 - RAM, 14 - bus, 15 - I / O interface, 16 - input unit, 17 - output unit, 18 - storage unit, 19 - communication unit;

[0043] 100-change rate acquisition module, 200-change coefficient acquisition module, 300-prediction module.

[0044] The present invention is further described in detail below. However, the following examples are merely simplified examples of the present invention and do not represent or limit the scope of protection of the present invention. The scope of protection of the present invention shall be subject to the claims. DETAILED DESCRIPTION

[0045] To better illustrate the present invention and facilitate understanding of the technical solutions of the present invention, typical but non-limiting embodiments of the present invention are as follows:

[0046] This embodiment provides a method for predicting the thickness of salt layers in salt structure restoration. The process is as follows: Figure 1 As shown, the prediction method includes:

[0047] Based on the initial salt layer thickness and the current salt layer thickness of the salt structure in the study area, the average salt layer change rate S is calculated;

[0048] Based on the basin evolution background, special events in the salt basin evolution process are defined, and the variation coefficient of salt layer thickness under different special events is obtained;

[0049] Based on the average change rate of the salt layer and the change coefficient as a constraint, the salt layer thickness corresponding to the special event period is predicted.

[0050] The calculation formula of the average change rate S of the salt layer is as follows:

[0051] S=|Nn| / T

[0052] Where N is the initial salt layer thickness, m; n is the current salt layer thickness, m; and T is the evolution time to date, myr.

[0053] Among them, the definition of special events in the evolution of salt basins based on the basin evolution background includes: establishing a complete basin dynamics framework and sedimentary stratigraphic model based on the geophysical data of the study area, and determining the key influencing events of salt basin development in different periods.

[0054] In the present invention, special events of salt basin development are defined, such as carbonate development and thermal deposition, which can be regarded as special events. Each special factor dominates the change in salt layer thickness in the corresponding period and corresponds to a different salt thickness variation coefficient. In addition, more special events can be set according to actual conditions.

[0055] The geophysical data include: well logging data and seismic data.

[0056] The variation coefficient includes: an initial coefficient a, a heat sink coefficient b, or a welding coefficient c, or a combination of at least two of them.

[0057] The initial coefficient a includes: initial clastic rock layer coefficient a 碎屑 .

[0058] Among them, the initial clastic rock layer coefficient a 碎屑=β×S, where, if the thickness of the clastic rock is less than the initial salt layer thickness N, β=k×clastic rock thickness, k=0.00025, β<0.5; if the thickness of the clastic rock is greater than or equal to the initial salt layer thickness N, β is equal to 0.5.

[0059] Among them, the heat sink coefficient b includes: based on the substrate lifting angle during the heat sinking period, when the substrate lifting angle is 2°, b is 3; based on the substrate lifting angle of 2°, b increases by 2 for every 1° increase in the substrate lifting angle; when the substrate lifting angle is greater than 4°, b is 7; when the substrate lifting angle is less than 2°, b decreases by 1 for every 1° decrease in the substrate lifting angle.

[0060] The welding coefficient c includes: based on the welding residual thickness, if the welding residual thickness is ≤20m, then c is 0.01; if the welding residual thickness is >20m, then c=0.05.

[0061] The relationship used to predict the salt layer thickness corresponding to the special event period includes:

[0062] S 初期 ×t1+S 热沉降 ×t2+S 焊接后 ×t3+S * ×t n =H, where, when the rock layer is a carbonate rock layer, S 初期 =0.1×(h 初期碳酸 / 400m)×H / t 初期碳酸 , H is the total salt layer thickness change, m; h 初期碳酸 is the thickness of the carbonate formation determined by actual core sampling or seismic profile, m; t 初期碳酸 It is determined by actual core sampling or seismic profile, which is obtained by subtracting the age of the bottom from the age of the top of the carbonate rock. If the rock layer is a clastic rock layer, then S 初期 = initial coefficient a × average rate of change S; S 热沉降 = heat sink coefficient b×average rate of change S; S 焊接后 = welding coefficient c×average rate of change S; S * = other coefficients of change × average rate of change S; t1 + t2 + t3 + ... t n =T, t1 is the sedimentation time corresponding to the post-salt carbonate stratum or the sedimentation time corresponding to the thickness of the clastic sedimentary rock reaching the critical point of density inversion, t2 is the duration of basin thermal subsidence, t3 is the duration from the formation of salt welding to the present, t n is the era corresponding to other coefficients of variation.

[0063] In the present invention, the specific value of time t1 must depend on actual geological data such as core sampling seismic profiles because the thickness of sedimentary strata in different basins varies greatly.

[0064] In the present invention, time t2 depends on the duration of the rift and the intensity of the magma intrusion activity. Since thermal subsidence lasts for a long time and is relatively slow in the later period, the value of t2 only considers the rapid subsidence stage in the early and middle stages, and the time span is generally 10-50 million years. Rapid rift extension or strong magma intrusion usually corresponds to a smaller time value, and vice versa.

[0065] In the present invention, time t3 is used to determine the formation time of salt welding, which is usually determined based on the characteristics of the top sedimentary strata. In areas without external tectonic forces, welding will cause the salt flow to stop and the strata above the salt to remain stable. When external stress exists, the bottom salt layer cannot flow, and the stratum deformation cannot be adjusted by the salt rock, forming structural features similar to those of a non-salt basin. Therefore, t3 is usually determined based on the time span of the corresponding sedimentary strata based on the actual seismic profile.

[0066] Furthermore, this embodiment provides a device for predicting the thickness of salt layers in salt structure restoration, such as Figure 2 As shown, the prediction device includes:

[0067] The change rate acquisition module 100 is used to calculate the average change rate S of the salt layer based on the initial salt layer thickness and the current salt layer thickness of the salt structure in the study area;

[0068] The variation coefficient acquisition module 200 is used to define special events in the salt basin evolution process based on the basin evolution background and obtain the variation coefficient of the salt layer thickness under different special events;

[0069] The prediction module 300 is used to predict the thickness of the salt layer corresponding to the special event period based on the average change rate of the salt layer and the change coefficient as a constraint.

[0070] Regarding the apparatus in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.

[0071] Further, the present invention provides an electronic device, which 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 may also represent various forms of mobile devices, such as personal digital assistants, cellular phones, smart phones, wearable devices (such as helmets, glasses, watches, etc.) and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the invention described and / or claimed herein.

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

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

[0074] The processor 11 can be any general-purpose and / or specialized processing component with processing and computing capabilities. Some examples of the processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various specialized artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as the aforementioned method for predicting salt layer thickness in salt structure restoration.

[0075] In some embodiments, the aforementioned method for predicting the thickness of the salt layer in salt structure restoration can be implemented as a computer program, which is tangibly contained in a computer-readable storage medium, such as a storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed on the electronic device 10 via the ROM 12 and / or the communication unit 19. When the computer program is loaded into the RAM 13 and executed by the processor 11, one or more steps of the aforementioned method for predicting the thickness of the salt layer in salt structure restoration described above can be performed. Alternatively, in other embodiments, the processor 11 can be configured to execute the aforementioned method for predicting the thickness of the salt layer in salt structure restoration in any other appropriate manner (for example, by means of firmware).

[0076] Various embodiments of the systems and techniques described 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), system-on-chip systems (SOCs), programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various embodiments can include being implemented in one or more computer programs that are executable and / or interpreted on a programmable system that includes at least one programmable processor, which can be a special purpose or general purpose programmable processor that can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit data and instructions to the storage system, the at least one input device, and the at least one output device.

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

[0078] In the context of the present invention, computer-readable storage media can be tangible media that can contain or store a computer program for use with an instruction execution system, device or equipment or used in combination with an instruction execution system, device or equipment. Computer-readable storage media can include but are not limited to electronic, magnetic, optical, electromagnetic, infrared or semiconductor systems, devices or equipment, or any suitable combination of the foregoing. Alternatively, computer-readable storage media can be machine-readable signal media. More specific examples of machine-readable storage media can include electrical connections based on one or more lines, portable computer disks, hard disks, random access memories (RAM), read-only memories (ROM), erasable programmable read-only memories (EPROM or flash memory), optical fibers, portable compact disk read-only memories (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination of the foregoing.

[0079] 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 can provide input to the electronic device. Other types of devices can also be used to provide interaction with the user; for example, the 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 acoustic input, voice input, or tactile input).

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

[0081] A computing system may include clients and servers. The clients and servers are typically remote from each other and typically interact via a communication network. This client-server relationship arises through computer programs running on the respective computers, creating a client-server relationship. The server may be a cloud server, also known as a cloud computing server or cloud host. This server is a hosting product within the cloud computing service ecosystem that addresses the management difficulties and limited scalability of traditional physical hosting and VPS services.

[0082] The server provided in this embodiment includes: a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, the method for predicting the thickness of the salt layer in the salt structure restoration is implemented.

[0083] Unless otherwise specifically stated, terms such as process, calculate, compute, determine, display, and the like may refer to the actions and / or processes of one or more processing or computing systems, or similar devices, that manipulate and convert data represented as physical (e.g., electronic) quantities within registers or memories of a processing system into other data similarly represented as physical quantities within the memories, registers, or other such information storage, transmission, or display devices of the processing system. Information and signals may be represented using any of a variety of different techniques and methods. For example, data, instructions, commands, information, signals, bits, symbols, and chips referred to throughout the above description may be represented by voltages, currents, electromagnetic waves, magnetic fields or particles, light fields or particles, or any combination thereof.

[0084] Those skilled in the art will also appreciate that the various illustrative logic blocks, modules, circuits, and algorithmic steps described in conjunction with the embodiments of the present invention may be implemented as electronic hardware, computer software, or a combination thereof. In order to clearly illustrate the interchangeability between hardware and software, the various illustrative components, blocks, modules, circuits, and steps described above are generally described with respect to their functions. Whether such functions are implemented as hardware or software depends on the specific application and the design constraints imposed on the entire system. A skilled person may implement the described functions in a flexible manner for each specific application, but such implementation decisions should not be interpreted as departing from the scope of protection of the present invention.

[0085] The steps of the methods or algorithms described in conjunction with the embodiments herein may be directly embodied as hardware, software modules executed by a processor, or a combination thereof. The software module may be located in a RAM memory, flash memory, ROM memory, EPROM memory, EEPROM memory, register, hard disk, removable disk, CD-ROM, or any other form of storage medium well known in the art. An exemplary storage medium is connected to the processor so that the processor can read information from the storage medium and write information to the storage medium. Of course, the storage medium may also be an integral part of the processor. The processor and storage medium may be located in an ASIC. The ASIC may be located in a user terminal. Of course, the processor and storage medium may also be present in a user terminal as discrete components.

[0086] For software implementation, the techniques described herein can be implemented using modules (e.g., procedures, functions, etc.) that perform the functions described herein. These software codes can be stored in a memory unit and executed by a processor. The memory unit can be implemented within the processor or external to the processor. In the latter case, it is communicatively coupled to the processor via various means, which are well known in the art.

[0087] Furthermore, in order to illustrate the effect that can be achieved by the prediction method provided by the present invention, actual examples are used for illustration, as follows:

[0088] Example 1

[0089] This embodiment provides a process for predicting the thickness of salt layers in salt structure restoration, which is as follows:

[0090] A seismic profile of a salt basin was selected, and through geological investigation and analysis and interpretation of geological data, the following basic information was obtained: the profile is located in the extension zone of the salt basin, which means that the thickness of the salt layer has been decreasing from ancient times to the present. Secondly, the initial thickness of the salt layer was set to 1.5km, but due to the existence of salt welding, the current thickness is about 50m. At the same time, a large amount of carbonate rocks developed in the early post-salt period. The evolution time of the entire salt basin is 115myr, and there was a significant basement uplift caused by thermal subsidence at about 100myr. Therefore, three special events were defined, namely salt rock welding, basement uplift, and carbonate rock development. Figure 4 and Figure 5 The calculation process of salt rock thickness change is as follows:

[0091] (1) The average rate of change S is 1.45 km / 115 myr = 0.017 km / myr.

[0092] (2) Through high-precision seismic identification, the salt weld location and thickness are determined. After determining the salt weld location, the characteristic stratum corresponding to the upper part of the salt weld is determined, that is, the structurally stable stratum. It is found that the sedimentary strata above it from 85myr to the present show good stability. Therefore, it can be considered that the salt weld was basically formed at 85myr, that is, t3 = 85myr. According to the current residual thickness of the salt weld identified by seismic identification is 50m, > 20m, at this time S 焊接后 =0.05×0.013km / myr=0.00065km / myr.

[0093] (3) The location of the profile in the study area is at the edge of the tension zone, and the left side is close to the transition zone. Therefore, the area affected by thermal subsidence is mainly in the right area outside the profile. When the basement on the right side is uplifted due to the influence of thermal subsidence, the amount of salt rock migrating downhill increases, and the thickness of the salt layer changes rapidly. The seismic profile shows that the current inclination angle of the basin is 3°. Since the change rates at 1°, 2° and 3° are different, they should be calculated in sections. The main period of thermal subsidence is 100myr-85myr, so t2=15myr. Assuming that the change rate of thermal subsidence for each degree is the same, according to the subsidence rates corresponding to different angles, the rate at 0° to 1° is S 热沉降1 =2×0.013km / myr=0.026km / myr, the speed at 1° to 2° is S 热沉降2=3×0.013km / myr=0.039km / myr, the speed at 2° to 3° is S 热沉降3 =5×0.013km / myr=0.065km / myr.

[0094] (4) According to the core and logging data, thick carbonate rocks developed in the early post-salt stage, and no clastic rocks developed. The development time of carbonate rocks is from 115myr to 100myr, so t1 = 15myr. At the same time, the post-salt carbonate rocks in this area are thick and dense. Therefore, the carbonate rocks in the early post-salt stage can quickly drive the deformation and flow of salt rocks, which is the main driving factor for the formation of salt diapir and dominates the large changes in the thickness of the salt layer. Therefore, substituting the data into S 碳酸 =0.1×(h 碳酸 / 400m)×H / t 初期碳酸 , then S 碳酸 =0.1×(2km / 0.4km)×1.5km / 15myr=0.05km / myr.

[0095] (5) Calculate the total thickness of each event: S 初期 ×t1+S 热沉降 ×t2+S 焊接后 ×t3=S 碳酸 ×t1+(S 热沉降1 ×t2 / 3+S 热沉降2 ×t2 / 3+S 热沉降3 ×t2 / 3)+S 焊接后 ×t3=0.05km / myr×15myr+(0.026km / myr×5myr+0.039km / myr×5myr+0.065km / myr×5myr)+0.00065km / myr×85myr=1.457km≈H.

[0096] It can be seen that the total thickness variation is basically consistent. Finally, based on the calculated salt layer thickness variation, the specific profile restoration is brought in and the local details can be adjusted by removing the fault distance and layer flattening.

[0097] It can be seen that the prediction method provided by the present invention, by considering the impact of special events in salt structures, targets the thickness changes caused by the non-conservation and viscoplastic deformation of salt rocks during the restoration of salt structures, constrains the thickness of the salt layer through geological data such as geodynamic background and lithology changes, realizes the evaluation of the salt layer thickness at different periods, and realizes the efficient prediction of the salt layer thickness changes during the restoration of salt structures.

[0098] It should be noted that the present invention uses the above-described embodiments to illustrate the detailed structural features of the present invention. However, the present invention is not limited to these detailed structural features, and this does not mean that the present invention must rely on these detailed structural features to be implemented. Those skilled in the art should understand that any improvements to the present invention, equivalent replacements for selected components, additions of auxiliary components, and selection of specific methods, etc., fall within the scope of protection and disclosure of the present invention.

[0099] The preferred embodiments of the present invention are described in detail above. However, the present invention is not limited to the specific details in the above embodiments. Within the technical concept of the present invention, various simple modifications can be made to the technical solution of the present invention, and these simple modifications all fall within the scope of protection of the present invention.

[0100] It should also be noted that the various specific technical features described in the above specific embodiments can be combined in any appropriate manner without contradiction. In order to avoid unnecessary repetition, the present invention will not further describe various possible combinations.

[0101] In addition, the various embodiments of the present invention may be arbitrarily combined, and as long as they do not violate the concept of the present invention, they should also be regarded as the contents disclosed by the present invention.

Claims

1. A method for predicting salt layer thickness in salt structure restoration, characterized in that: The prediction method comprises: Based on the initial salt layer thickness and the current salt layer thickness of the salt structure in the study area, the average salt layer change rate S is calculated; Based on the basin evolution background, special events in the salt basin evolution process are defined, and the variation coefficient of salt layer thickness under different special events is obtained; Based on the average change rate of the salt layer and the change coefficient as a constraint, the salt layer thickness corresponding to the special event period is predicted.

2. The prediction method according to claim 1, wherein: The calculation formula of the average change rate S of the salt layer is as follows: S=|Nn| / T Where N is the initial salt layer thickness, m; n is the current salt layer thickness, m; and T is the evolution time to date, myr.

3. The prediction method according to claim 1 or 2, wherein: The definition of special events in the evolution of salt basins based on the basin evolution background includes: establishing a complete basin dynamics framework and sedimentary stratigraphic model based on the geophysical data of the study area, and determining the key influencing events in the development of salt basins in different periods; Preferably, the geophysical data includes: well logging data and seismic data.

4. The prediction method according to any one of claims 1 to 3, wherein: The variation coefficient includes: an initial coefficient a, a heat sink coefficient b, or a welding coefficient c, or a combination of at least two of them.

5. The prediction method according to claim 4, wherein: The initial coefficient a includes: initial clastic rock layer coefficient a 碎屑 ; Preferably, the initial clastic rock layer coefficient a 碎屑 =β×S, where, if the thickness of the clastic rock is less than the initial salt layer thickness N, β=k×clastic rock thickness, k=0.00025, β<0.5; if the thickness of the clastic rock is greater than or equal to the initial salt layer thickness N, β is equal to 0.

5.

6. The prediction method according to claim 4 or 5, characterized in that: The heat sink coefficient b includes: based on the base lifting angle during the heat sinking period, when the base lifting angle is 2°, b is 3, based on the base lifting angle of 2°, b increases by 2 for every 1° increase in the base lifting angle; when the base lifting angle is greater than 4°, b is 7; when the base lifting angle is less than 2°, b decreases by 1 for every 1° decrease in the base lifting angle; Preferably, the welding coefficient c includes: based on the welding residual thickness, if the welding residual thickness is ≤20m, then c is 0.01; if the welding residual thickness is >20m, then c=0.

05.

7. The prediction method according to any one of claims 1 to 6, wherein: The relationship used to predict the salt layer thickness corresponding to the special event period includes: S 初期 ×t1+S 热沉降 ×t2+S 焊接后 ×t3+S * ×t n =H In the formula, when the rock layer is a carbonate rock layer, S 初期 =0.1×(h 初期碳酸 / 400m)×H / t 初期碳酸 , H is the total salt layer thickness change, m; h 初期碳酸 is the thickness of the carbonate formation determined by actual core sampling or seismic profile, m; t 初期碳酸 It is determined by actual core sampling or seismic profile, which is obtained by subtracting the age of the bottom from the age of the top of the carbonate rock. If the rock layer is a clastic rock layer, then S 初期 = initial coefficient a × average rate of change S; S 热沉降 = heat sink coefficient b×average rate of change S; S 焊接后 = welding coefficient c×average rate of change S; S * = other coefficients of change × average rate of change S; t1 + t2 + t3 + ... t n =T, t1 is the sedimentation time corresponding to the post-salt carbonate stratum or the sedimentation time corresponding to the thickness of the clastic sedimentary rock reaching the critical point of density inversion, t2 is the duration of basin thermal subsidence, t3 is the duration from the formation of salt welding to the present, t n is the era corresponding to other coefficients of variation.

8. A device for predicting salt layer thickness in salt structure restoration, characterized in that: The prediction device comprises: The change rate acquisition module is used to calculate the average change rate S of the salt layer based on the initial salt layer thickness and the current salt layer thickness of the salt structure in the study area; The variation coefficient acquisition module is used to define special events in the evolution of salt basins based on the basin evolution background and obtain the variation coefficient of salt layer thickness under different special events; The prediction module is used to predict the salt layer thickness corresponding to the special event period based on the average change rate of the salt layer and the change coefficient as a constraint.

9. An electronic device, characterized in that: The electronic device comprises: at least one processor; and a memory communicatively coupled to the at least one processor; Wherein, the memory stores a computer program that can be executed by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the method for predicting salt layer thickness in salt structure restoration according to any one of claims 1-7.

10. A computer storage medium, characterized in that The computer storage medium stores computer executable instructions, which, when executed by a processor, implement the method for predicting salt layer thickness in salt structure restoration according to any one of claims 1 to 7.

Images