Method and device for correcting a point of closure of a volcanic edifice trap

Abstract

The present application provides a kind of volcano mechanism trap boundary point correction method and correction device, the correction method includes: obtaining the initial boundary point of volcano mechanism;Based on well logging data and three-dimensional seismic data, establish three-dimensional geological model;Based on the three-dimensional geological model, seismic forward processing is carried out, and the seismic attribute parameter based on the three-dimensional geological model is obtained;Based on each initial boundary point and the seismic attribute parameter based on the three-dimensional geological model, according to the preset rule, each initial boundary point is corrected, and the initial boundary point after correction is used as volcano mechanism boundary point.The volcano mechanism trap boundary point correction method and correction device have the advantages of reducing the identification error of volcano mechanism trap boundary, improving the identification accuracy, thereby more accurately identifying the trap boundary of volcano mechanism, effectively improving the identification accuracy and reliability of volcanic rock oil and gas accumulation boundary, reducing the exploration risk, and laying a foundation for subsequent deployment and more accurate calculation resource quantity.

Description

Technical Field

[0001] This invention relates to the field of volcanic structure boundary recognition technology, specifically to a method for correcting the closed boundary points of a volcanic structure, a device for correcting the closed boundary of a volcanic structure, a computer device, and a machine-readable storage medium. Background Technology

[0002] Due to the unique formation mechanism of volcanic structures, they often exhibit unconformity with sedimentary strata. Conventional seismic interpretation traces changes in the reflection phase axis of sedimentary strata surrounding volcanic rocks as the boundary of the volcanic structure to delineate its extent. However, because volcanic and sedimentary rocks often exhibit strong reflection interfaces, the unconformity points are located within the strong reflection zone, resulting in significant errors in conventional seismic interpretation.

[0003] In the mid-to-late 1990s, the exploration of igneous rock oil and gas reservoirs made rapid progress. During this period, numerous identification methods emerged, including horizontal slice comparison analysis, three-instantaneous seismic attribute analysis, coherence volume and curvature, and absorption attenuation detection, all of which were widely used. However, there were not many methods for characterizing the extent of volcanic rock traps because there were no clear geophysical parameters that could directly identify volcanic rocks. Most of the characterization was done qualitatively in attribute volumes and inversion volumes, and quantitative characterization methods have not yet been studied. This has led to difficulties in identifying the extent of volcanic rock traps, low identification accuracy, and increased exploration risks. Summary of the Invention

[0004] The purpose of this invention is to provide a method and device for correcting the boundary points of volcanic traps, so as to at least solve the problems of difficulty in identifying the range of volcanic rock traps, low identification accuracy, and increased exploration risks.

[0005] To achieve the above objectives, a first aspect of the present invention provides a method for correcting the closed boundary points of a volcanic mechanism, comprising:

[0006] Obtain the initial boundary points of the volcanic structure;

[0007] A three-dimensional geological model was established based on well logging data and three-dimensional seismic data.

[0008] Seismic forward modeling is performed based on the three-dimensional geological model to obtain seismic attribute parameters based on the three-dimensional geological model.

[0009] Based on each initial boundary point and the seismic attribute parameters based on the three-dimensional geological model, each initial boundary point is corrected according to a preset rule, and the corrected initial boundary point is used as the boundary point of the volcanic mechanism.

[0010] Optionally, the logging data includes logging data from volcanic structures and logging data from unconformities.

[0011] Optionally, the three-dimensional geological model includes: an x-axis, a y-axis, and a z-axis;

[0012] The process of obtaining seismic attribute parameters based on the three-dimensional geological model includes:

[0013] Based on each initial boundary point, a section is cut perpendicular to the y-axis of the three-dimensional geological model to obtain multiple two-dimensional slice layers;

[0014] The seismic attribute parameters of each two-dimensional slice layer are obtained as the seismic attribute parameters of the three-dimensional geological model.

[0015] Optionally, the seismic attribute parameters include at least: layer velocity, dominant seismic frequency, angle between the unconformity and the horizontal plane, and angle between the unconformity and the surface of the volcanic structure;

[0016] The unconformity surface is the contact surface between the unconformity stratum and the overlying stratum.

[0017] Optionally, obtaining the seismic attribute parameters for each two-dimensional slice layer includes:

[0018] Based on 3D seismic data, the dominant seismic frequency of each 2D slice layer is extracted. The extraction time window includes one apparent period of the seismic reflection axis but does not include the reflection axes of other 2D slice layers.

[0019] Based on well logging data, the layer velocity of each two-dimensional slice layer is read;

[0020] Based on well logging data and seismic data, the angle between the unconformity surface and the horizontal plane, as well as the angle between the unconformity surface and the surface of the volcanic structure, of each two-dimensional slice layer were calculated.

[0021] Optionally, the step of correcting each initial boundary point according to a preset rule includes:

[0022]

[0023] Among them, X i The corrected initial boundary point, i.e., the boundary point of the volcanic mechanism; x i α is the initial boundary point; f is the dominant seismic frequency; v is the layer velocity; α i The angle between the unconformity surface and the horizontal plane; β i The angle between the unconformity surface and the surface of the volcanic structure.

[0024] A second aspect of the present invention provides a volcanic mechanism trap boundary point correction device, comprising:

[0025] The initial point acquisition module is used to acquire the initial boundary points of the volcano mechanism;

[0026] The model building module is used to build three-dimensional geological models based on well logging data and seismic data;

[0027] The parameter acquisition module is used to perform seismic forward modeling based on the three-dimensional geological model and obtain seismic attribute parameters based on the three-dimensional geological model.

[0028] The correction module is used to correct each initial boundary point according to a preset rule based on each initial boundary point and the seismic attribute parameters based on the three-dimensional geological model, and to use the corrected initial boundary point as the boundary point of the volcanic mechanism.

[0029] Optionally, the three-dimensional geological model includes: an x-axis, a y-axis, and a z-axis;

[0030] The process of obtaining seismic attribute parameters based on the three-dimensional geological model includes:

[0031] Based on each initial boundary point, a section is cut perpendicular to the y-axis of the three-dimensional geological model to obtain multiple two-dimensional slice layers;

[0032] The seismic attribute parameters of each two-dimensional slice layer are obtained as the seismic attribute parameters of the three-dimensional geological model.

[0033] Optionally, the seismic attribute parameters include at least: layer velocity, dominant seismic frequency, angle between the unconformity and the horizontal plane, and angle between the unconformity and the surface of the volcanic structure;

[0034] The unconformity surface is the contact surface between the unconformity stratum and the overlying stratum.

[0035] Optionally, the parameter acquisition module includes:

[0036] The seismic dominant frequency acquisition module is used to extract the seismic dominant frequency of each two-dimensional slice layer based on three-dimensional seismic data. The extraction time window includes one apparent period of the seismic reflection axis and does not include the reflection axes of other two-dimensional slice layers.

[0037] The layer velocity acquisition module is used to read the layer velocity of each two-dimensional slice layer based on well logging data;

[0038] The angle acquisition module is used to calculate the angle between the unconformity and the horizontal plane, as well as the angle between the unconformity and the surface of the volcanic structure, for each two-dimensional slice layer based on well logging data and seismic data.

[0039] Optionally, the step of correcting each initial boundary point according to a preset rule includes:

[0040]

[0041] Among them, X i The corrected initial boundary point, i.e., the boundary point of the volcanic mechanism; x i α is the initial boundary point; f is the dominant seismic frequency; v is the layer velocity; αi The angle between the unconformity surface and the horizontal plane; β i The angle between the unconformity surface and the surface of the volcanic structure.

[0042] A third aspect of the present invention provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the above-described method for correcting the closed boundary points of a volcano mechanism.

[0043] On the other hand, the present invention provides a machine-readable storage medium storing instructions for causing a machine to execute the above-described volcano mechanism trap boundary point correction method.

[0044] This technical solution corrects the boundary points of volcanic structures obtained from seismic data to obtain updated boundary points, which can reduce the identification error of volcanic structure trap boundaries, improve identification accuracy, and thus more accurately identify the trap boundaries of volcanic structures. This effectively improves the identification accuracy and reliability of volcanic rock oil and gas accumulation boundaries, reduces exploration risks, and lays the foundation for subsequent deployment and more accurate calculation of resource quantities.

[0045] Other features and advantages of the embodiments of the present invention will be described in detail in the following detailed description section. Attached Figure Description

[0046] The accompanying drawings are provided to further illustrate embodiments of the present invention and form part of the specification. They are used together with the following detailed description to explain the embodiments of the present invention, but do not constitute a limitation thereof. In the drawings:

[0047] Figure 1 This is a flowchart of the volcanic mechanism trap boundary point correction method provided by the present invention;

[0048] Figure 2 This is a schematic diagram of the two-dimensional slice layer provided by the present invention;

[0049] Figure 3 This is a simplified schematic diagram of a two-dimensional slice layer provided by the present invention;

[0050] Figure 4 This is a schematic diagram of the closed boundary of the volcanic mechanism provided by the present invention;

[0051] Figure 5 This is a schematic diagram of the volcanic mechanism trap boundary point correction device provided by the present invention;

[0052] Figure 6 This is a schematic diagram of the parameter acquisition module of the volcanic mechanism trap boundary point correction device provided by the present invention.

[0053] Explanation of reference numerals in the attached figures

[0054] 11-Initial point acquisition module; 12-Model building module; 13-Parameter acquisition module;

[0055] 14-Correction Module; 131-Seismic Dominant Frequency Acquisition Module; 132-Layer Velocity Acquisition Module;

[0056] 133 - Angle Acquisition Module. Detailed Implementation

[0057] The specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings. It should be understood that the specific embodiments described herein are for illustration and explanation only and are not intended to limit the present invention.

[0058] Figure 1 This is a flowchart of the volcanic mechanism trap boundary point correction method provided by the present invention. Figure 1 As shown, this invention provides a method for correcting the boundary points of a volcanic mechanism trap, the method comprising:

[0059] Step 101: Obtain the initial boundary points of the volcano mechanism;

[0060] Step 102: Based on well logging data and 3D seismic data, establish a 3D geological model;

[0061] Step 103: Perform seismic forward modeling based on the three-dimensional geological model to obtain seismic attribute parameters based on the three-dimensional geological model;

[0062] Step 104: Based on each initial boundary point and the seismic attribute parameters based on the three-dimensional geological model, each initial boundary point is corrected according to a preset rule, and the corrected initial boundary point is used as the boundary point of the volcanic mechanism.

[0063] Specifically, obtaining the initial boundary points of the volcanic structure includes: initial boundary points of the volcanic structure obtained from seismic data. These initial boundary points are synthesized by seismic interpretation software. Since they are obtained solely from seismic data and are affected by factors such as human experience, there is a significant error. Therefore, correction is required to obtain more accurate boundary points.

[0064] Three-dimensional geological models can be established using a geogeological modeling system based on well logging data and 3D seismic data. This system is a geological modeling software designed for describing the structure and reservoir morphology of oil and gas reservoirs during the development stage of oil and gas fields. It utilizes geostatistics, planar and spatial network technologies, and integrates technologies such as well logging, well logging, comprehensive geological analysis, and rapid and efficient 3D visualization. Relying on an efficient and integrated interpretation process, a unified application environment, and an optimized data management framework, it ultimately achieves a quantitative spatial description of oil and gas reservoirs. The system provides functions such as data management and work area management for application, fault modeling, structural modeling, stratigraphic facies modeling, attribute modeling, oil and gas reservoir modeling, and digital automatic modeling.

[0065] Seismic attribute parameters can be directly read from well logging data and 3D seismic data or obtained through data calculation. Furthermore, the 3D geological model contains seismic attribute parameters. Each initial boundary point may have the same or different seismic attribute parameters, which are determined based on the location of the initial boundary point.

[0066] In another implementation, before establishing a three-dimensional geological model, the acquired logging data undergoes consistency processing. Due to the large time span of logging data acquisition, numerous logging instrument signals, inconsistent calibration standards, and inconsistent operating methods, in order to eliminate the different systematic errors between logging data measured at different times and with different instruments, consistency processing is usually performed on the logging data before using the data for modeling. Specifically, this includes using the mean-variance method to analyze the mean and variance of the logging curves, and relatively unifying the logging data measured at different times and with different instruments.

[0067] Furthermore, the three-dimensional geological model is established based on well logging data and three-dimensional seismic data through well logging interpretation and seismic forward modeling techniques. The well logging data includes well logging data of volcanic structures and well logging data of unconformities.

[0068] The unconformity refers to strata with sedimentary discontinuities or missing strata.

[0069] Specifically, in the process of oil and gas exploration, well logging data and 3D seismic data are usually used. Well logging data is usually obtained by setting up at least one acquisition well in different strata at the acquisition site to obtain well logging data of each stratum and obtain well logging curves. 3D seismic data is usually obtained by laying out a certain number of excitation points and receiver points in an area manner according to the pre-set acquisition requirements and collecting seismic data on the strata.

[0070] After obtaining the logging data, it is processed using logging interpretation, also known as well logging comprehensive interpretation. The core of well logging interpretation is to determine the relationship between well logging information and geological information. The correct method is used to process the well logging information into geological information, which is often used for oil and gas exploration. Various physical parameters in the well logging data, such as resistivity, spontaneous potential, sonic velocity, and rock bulk density, are converted into geological information such as lithology, clay content, water saturation, and permeability.

[0071] Seismic forward modeling, also known as forward simulation, is a method that uses numerical calculations to study the propagation of seismic waves in the subsurface medium when the structure and parameters of the subsurface medium are known, thereby obtaining theoretical seismic records.

[0072] A high-precision model is established using 3D seismic data and well logging data to provide a 3D visualization platform for further reservoir research. In this embodiment, a sequence structure model is established based on each formation stratum, and an attribute model is established using the total high-resolution well logging curves.

[0073] Furthermore, the three-dimensional geological model includes: an x-axis, a y-axis, and a z-axis;

[0074] The process of obtaining seismic attribute parameters based on the three-dimensional geological model includes:

[0075] Based on each initial boundary point, a section is cut perpendicular to the y-axis of the three-dimensional geological model to obtain multiple two-dimensional slice layers;

[0076] The seismic attribute parameters of each two-dimensional slice layer are obtained as the seismic attribute parameters of the three-dimensional geological model.

[0077] Specifically, the x-axis, y-axis, and z-axis in the three-dimensional geological model are set in the directions conventionally used by those skilled in the art. After obtaining the three-dimensional geological model, in order to further understand the spatial distribution of the three-dimensional geological model in this embodiment, a section is cut perpendicular to the y-axis of the three-dimensional geological model based on each initial boundary point to obtain multiple geological profiles of the three-dimensional geological model as two-dimensional slices. The obtained multiple geological profiles can clearly show the spatial relationship between the various strata. By conducting close and detailed analysis on each two-dimensional slice, the position of the initial boundary point contained in the corresponding two-dimensional slice can be corrected.

[0078] Since the dip angles between the volcanic structure and the unconformity layer, and between the overlying strata and the unconformity layer, may vary slightly at different locations, the three-dimensional geological model can be sectioned based on the initial boundary points to obtain two-dimensional slices. This allows us to obtain the seismic attribute parameters contained in the two-dimensional slices corresponding to each initial boundary point. The seismic attribute parameters of the two-dimensional slices corresponding to each initial boundary point can then be used as the seismic attribute parameters of the three-dimensional geological model for subsequent calculations, thereby correcting the initial boundary points in the corresponding two-dimensional slices.

[0079] Furthermore, the seismic attribute parameters include at least: layer velocity, dominant seismic frequency, angle between the unconformity and the horizontal plane, and angle between the unconformity and the surface of the volcanic structure;

[0080] The unconformity surface is the contact surface between the unconformity stratum and the overlying stratum.

[0081] Further, obtaining the seismic attribute parameters of the three-dimensional geological model includes:

[0082] Based on 3D seismic data, the dominant seismic frequency of each 2D slice layer is extracted. The extraction time window includes one apparent period of the seismic reflection axis but does not include the reflection axes of other 2D slice layers.

[0083] Based on well logging data, the layer velocity of each two-dimensional slice layer is read;

[0084] Based on well logging data and seismic data, the angle between the unconformity surface and the horizontal plane, as well as the angle between the unconformity surface and the surface of the volcanic structure, of each two-dimensional slice layer were calculated.

[0085] Specifically, earthquake attribute parameters can be directly output through earthquake interpretation software;

[0086] Based on the well logging data and seismic data, when calculating the angle between the unconformity surface of the two-dimensional slice layer corresponding to each initial boundary point and the horizontal plane, as well as the angle between the unconformity surface and the surface of the volcanic structure, a large grid spacing parameter can be initially used to ensure that low-frequency information is present when calculating the angle between the unconformity surface and the horizontal plane, as well as the angle between the unconformity surface and the surface of the volcanic structure, making the results more accurate. After the preliminary calculation, a small grid parameter is then used to calculate to obtain accurate angle information.

[0087] Furthermore, the step of correcting each initial boundary point according to a preset rule includes:

[0088]

[0089] Among them, X i The corrected initial boundary point, i.e., the boundary point of the volcanic mechanism; x i α is the initial boundary point; f is the dominant seismic frequency; v is the layer velocity; α iThe angle between the unconformity surface and the horizontal plane; β i The angle between the unconformity surface and the surface of the volcanic structure.

[0090] Specifically, Figure 2 This is a schematic diagram of the two-dimensional slice layer provided by the present invention. Figure 3 This is a simplified schematic diagram of a two-dimensional slice layer provided by the invention; such as Figure 2 As stated above, after obtaining the three-dimensional geological model, the y-axis section of the three-dimensional geological model is used to obtain... Figure 2 After the two-dimensional slice layer shown, Figure 2 Simplified schematic diagram of two-dimensional slice layers Figure 3 , Figure 3 In this context, the length of OA is the difference between the updated boundary point and the initial boundary point's x-coordinate, and the length of AC is the difference between the updated boundary point and the initial boundary point's z-coordinate.

[0091] The specific derivation steps are as follows:

[0092] Where v is in m / s and f is in Hz = 1 / s;

[0093] ∠OBD=∠AOB=α i

[0094] ∠EBC=∠BDE=∠BOD+∠OBD=α i +β i

[0095]

[0096]

[0097]

[0098]

[0099]

[0100] Similarly, the length of AC can be obtained.

[0101] AC = OC × sin(∠AOC)

[0102]

[0103] The boundary point update formula can be obtained as follows:

[0104]

[0105] Among them, X i The corrected initial boundary point, i.e., the boundary point of the volcanic mechanism; x iLet f be the initial boundary point, f be the dominant seismic frequency, v be the layer velocity, and α be the initial boundary point. i β is the angle between the unconformity surface and the horizontal plane. i The angle between the unconformity and the volcanic structure;

[0106] In another implementation, Figure 4 This is a schematic diagram of the closed boundary of the volcanic mechanism provided by the present invention; as shown. Figure 4 As shown, based on the above calculation formula, all initial boundary points can be updated, and by sequentially connecting the updated initial boundary points, the closed boundary of the volcanic mechanism can be obtained. The area of ​​the closed boundary of the volcanic mechanism can be obtained through the closed boundary of the volcanic mechanism, which can be calculated using the following formula:

[0107] S 更新后 =S 初始 -ΔS

[0108]

[0109] Among them, S 更新后 To update the area of ​​the closed boundary of the volcanic structure, S 初始 Let X be the area of ​​the initial closed boundary of the volcanic structure. i For the updated boundary points of the volcanic mechanism, x i Let Δy be the initial boundary point. i This represents the channel spacing between two adjacent two-dimensional slice layers.

[0110] Figure 5 This is a schematic diagram of the volcano mechanism trap boundary point correction device provided by the present invention. Figure 5 As shown, an embodiment of the present invention provides a volcanic mechanism closed boundary point correction device, the device comprising:

[0111] Initial point acquisition module 11 is used to acquire the initial boundary points of the volcano mechanism;

[0112] Model building module 12 is used to build a three-dimensional geological model based on well logging data and seismic data;

[0113] The parameter acquisition module 13 is used to perform seismic forward modeling based on the three-dimensional geological model to obtain seismic attribute parameters based on the three-dimensional geological model.

[0114] The correction module 14 is used to correct each initial boundary point according to a preset rule based on each initial boundary point and the seismic attribute parameters based on the three-dimensional geological model, and to use the corrected initial boundary point as the boundary point of the volcanic mechanism.

[0115] Furthermore, the three-dimensional geological model includes: an x-axis, a y-axis, and a z-axis;

[0116] The process of obtaining seismic attribute parameters based on the three-dimensional geological model includes:

[0117] Based on each initial boundary point, slices are made perpendicular to the y-axis of the three-dimensional geological model to obtain multiple two-dimensional slice layers;

[0118] The seismic attribute parameters of each two-dimensional slice layer are obtained as the seismic attribute parameters of the three-dimensional geological model.

[0119] Furthermore, the seismic attribute parameters include at least: layer velocity, dominant seismic frequency, angle between the unconformity and the horizontal plane, and angle between the unconformity and the surface of the volcanic structure;

[0120] The unconformity surface is the contact surface between the unconformity stratum and the overlying stratum.

[0121] Figure 6 This is a schematic diagram of the parameter acquisition module of the volcanic mechanism trap boundary point correction device provided by the present invention, as shown below. Figure 6 As shown, the parameter acquisition module 13 further includes:

[0122] The seismic dominant frequency acquisition module 131 is used to extract the seismic dominant frequency of each two-dimensional slice layer based on three-dimensional seismic data, wherein the extraction time window includes one apparent period of the seismic reflection axis and does not include the reflection axes of the other two-dimensional slice layers.

[0123] Layer velocity acquisition module 132 is used to read the layer velocity of each two-dimensional slice layer based on well logging data;

[0124] Angle acquisition module 133 is used to calculate the angle between the unconformity and the horizontal plane, and the angle between the unconformity and the surface of the volcanic structure, based on well logging data and seismic data.

[0125] Furthermore, the step of correcting each initial boundary point according to a preset rule includes:

[0126]

[0127] Among them, X i The corrected initial boundary point, i.e., the boundary point of the volcanic mechanism; x i α is the initial boundary point; f is the dominant seismic frequency; v is the layer velocity; α i The angle between the unconformity surface and the horizontal plane; β i The angle between the unconformity surface and the surface of the volcanic structure.

[0128] A third aspect of the present invention provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the above-described method for correcting the closed boundary points of a volcano mechanism.

[0129] On the other hand, the present invention provides a machine-readable storage medium storing instructions for causing a machine to execute the above-described volcano mechanism trap boundary point correction method.

[0130] Example 1

[0131] This embodiment takes the intermediate-basic volcanic rocks of the Huoshiling Formation in the Chang'anhua area as an example to correct and update the boundary points of the volcanic structure, and obtain multiple initial boundary points of the volcanic structure obtained after seismic analysis.

[0132] A well, named Well 1, was set up in the Chang'anhua area to obtain logging data for the region. After analyzing the data from Well 1, the top interface of the volcanic structure was marked vertically based on the analysis results. After marking, the range of the volcanic structure and the unconformity strata of Well 1 were initially delineated based on the Chaganhua 3D seismic data. The 3D geological model formed by the volcanic structure and the unconformity strata was determined. This 3D geological model has x-axis, y-axis and z-axis, and the directions of the x-axis, y-axis and z-axis are set according to the conventional directions of those skilled in the art. With the y-axis as the reference axis, the 3D geological model was sliced ​​based on the initial boundary points to obtain multiple 2D slice layers. The x-axis coordinates and z-axis coordinates of each initial boundary point can be obtained through the 2D slice layers.

[0133] The dominant seismic frequency of the target layer was extracted from the seismic data. The time window range was 52ms, and the obtained dominant seismic frequency was 20Hz. There were no adjacent wells in the unconformity strata in this area. Therefore, the average layer velocity of 4000m / s was taken from the Yingcheng Formation above well Cha 1. The dip angle of the strata was calculated using OpenDtect software. The results showed that the average dip angle of the upper strata of the volcanic rocks in the north of well Cha 1 was 7.4°, which is the angle between the unconformity surface and the horizontal plane. The average dip angle of the unconformity strata was 5.3°, which is the angle between the unconformity surface and the surface of the volcanic structure.

[0134] The seismic dominant frequency of 20Hz, the formation velocity of 4000m / s, the angle between the unconformity and the horizontal plane of 7.4°, and the angle between the unconformity and the surface of the volcanic structure of 5.3° are substituted into the volcanic structure boundary point update formula for calculation. The calculated values ​​are used as the coordinate axes of the updated boundary points to correct the initial boundary points. The x and y coordinates of the initial boundary points are (588434, 4935543). Substituting them into the above formula, the x and y coordinates of the updated boundary points are (588173, 4935543). Since the profile is obtained along the y-axis, the y-axis coordinates of the two boundary points are the same. Since the volcanic structure trap boundary is usually calculated by projecting the shape of the volcanic structure in the three-dimensional geological model onto a two-dimensional plane to form the trap boundary, the z-axis coordinate is not needed in this embodiment. Therefore, the x-axis coordinate of the initial boundary point 588434 is updated using the x-axis coordinate of the updated boundary point 588173.

[0135] Those skilled in the art will understand that all or part of the steps in the methods of the above embodiments can be implemented by a program instructing related hardware. This program is stored in a storage medium and includes several instructions to cause a microcontroller, chip, or processor to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as a USB flash drive, a portable hard drive, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk.

[0136] The optional embodiments of the present invention have been described in detail above with reference to the accompanying drawings. However, the embodiments of the present invention are not limited to the specific details described above. Within the scope of the technical concept of the embodiments of the present invention, various simple modifications can be made to the technical solutions of the embodiments of the present invention, and these simple modifications all fall within the protection scope of the embodiments of the present invention. It should also be noted that the various specific technical features described in the above specific embodiments can be combined in any suitable manner without contradiction. To avoid unnecessary repetition, the embodiments of the present invention will not further describe the various possible combinations.

[0137] Furthermore, various different embodiments of the present invention can be combined in any way, as long as they do not violate the spirit of the embodiments of the present invention, they should also be regarded as the content disclosed by the embodiments of the present invention.

Claims

1. A method for correcting the boundary points of a volcanic structure trap, characterized in that, The method includes: Obtain multiple initial boundary points of the volcanic structure; A three-dimensional geological model was established based on well logging data and three-dimensional seismic data. Seismic forward modeling is performed based on the three-dimensional geological model to obtain seismic attribute parameters based on the three-dimensional geological model. The seismic attribute parameters include: layer velocity, seismic dominant frequency, angle between the unconformity and the horizontal plane, and angle between the unconformity and the surface of the volcanic structure. The unconformity is the contact surface between the unconformity strata and the overlying strata. Based on each initial boundary point and the seismic attribute parameters based on the three-dimensional geological model, each initial boundary point is corrected according to a preset rule, and the corrected initial boundary point is used as the boundary point of the volcanic structure; wherein, the correction of each initial boundary point according to the preset rule includes: in, The corrected initial boundary point, i.e., the boundary point of the volcanic mechanism; These are the initial boundary points; This is the dominant frequency of the earthquake. For layer velocity; The angle between the unconformity surface and the horizontal plane; The angle between the unconformity surface and the surface of the volcanic structure.

2. The method according to claim 1, characterized in that, The logging data includes logging data from volcanic structures and logging data from unconformities.

3. The method according to claim 1, characterized in that, The three-dimensional geological model includes an x-axis, a y-axis, and a z-axis; The process of obtaining seismic attribute parameters based on the three-dimensional geological model includes: Based on each initial boundary point, a section is cut perpendicular to the y-axis of the three-dimensional geological model to obtain multiple two-dimensional slice layers; The seismic attribute parameters of each two-dimensional slice layer are obtained as the seismic attribute parameters of the three-dimensional geological model.

4. The method according to claim 3, characterized in that, The process of obtaining the seismic attribute parameters for each two-dimensional slice layer includes: Based on 3D seismic data, the dominant seismic frequency of each 2D slice layer is extracted. The extraction time window includes one apparent period of the seismic reflection axis but does not include the reflection axes of other 2D slice layers. Based on well logging data, the layer velocity of each two-dimensional slice layer is read; Based on well logging data and 3D seismic data, the angle between the unconformity surface and the horizontal plane of each 2D slice layer, as well as the angle between the unconformity surface and the surface of the volcanic structure, are calculated.

5. A volcanic mechanism closed boundary point correction device, characterized in that, The device includes: The initial point acquisition module is used to acquire multiple initial boundary points of the volcano mechanism; The model building module is used to build three-dimensional geological models based on well logging data and seismic data; The parameter acquisition module is used to perform seismic forward modeling based on the three-dimensional geological model and acquire seismic attribute parameters based on the three-dimensional geological model. The seismic attribute parameters include: layer velocity, seismic dominant frequency, angle between the unconformity and the horizontal plane, and angle between the unconformity and the surface of the volcanic structure. The unconformity is the contact surface between the unconformity strata and the overlying strata. The correction module is used to correct each initial boundary point according to a preset rule based on each initial boundary point and the seismic attribute parameters based on the three-dimensional geological model, and to use the corrected initial boundary point as the volcanic mechanism boundary point; wherein, the correction of each initial boundary point according to the preset rule includes: in, The corrected initial boundary point, i.e., the boundary point of the volcanic mechanism; These are the initial boundary points; This is the dominant frequency of the earthquake. For layer velocity; The angle between the unconformity surface and the horizontal plane; The angle between the unconformity surface and the surface of the volcanic structure.

6. The apparatus according to claim 5, characterized in that, The three-dimensional geological model includes an x-axis, a y-axis, and a z-axis; The process of obtaining seismic attribute parameters based on the three-dimensional geological model includes: Based on each initial boundary point, a section is cut perpendicular to the y-axis of the three-dimensional geological model to obtain multiple two-dimensional slice layers; The seismic attribute parameters of each two-dimensional slice layer are obtained as the seismic attribute parameters of the three-dimensional geological model.

7. The apparatus according to claim 6, characterized in that, The parameter acquisition module includes: The seismic dominant frequency acquisition module is used to extract the seismic dominant frequency of each two-dimensional slice layer based on three-dimensional seismic data. The extraction time window includes one apparent period of the seismic reflection axis and does not include the reflection axes of other two-dimensional slice layers. The layer velocity acquisition module is used to read the layer velocity of each two-dimensional slice layer based on well logging data; The angle acquisition module is used to calculate the angle between the unconformity and the horizontal plane, as well as the angle between the unconformity and the surface of the volcanic structure, for each two-dimensional slice layer based on well logging data and three-dimensional seismic data.

8. A computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the volcanic mechanism closed boundary point correction method as described in any one of claims 1-4.

9. A machine-readable storage medium storing instructions for causing a machine to perform the volcano mechanism trap boundary point correction method according to any one of claims 1-4.

Images