Automatic layer tracking method, device and readable storage medium

Abstract

The present invention provides a method, device and readable storage medium for automatic layer tracking, which belongs to the field of petroleum geophysical exploration technology. The method includes: obtaining seismic data and the depth of the same layer in different input wells; determining the dip angle of the stratigraphic profile of each seismic trace based on the seismic data; determining the initial layer depth of each seismic trace based on the depth of the same layer in different input wells; constructing an objective function based on the dip angle of the stratigraphic profile of each seismic trace and the initial layer depth of each seismic trace, with the depth of the same layer in different input wells as a constraint condition; iteratively calculating the objective function, and iteratively updating the layer depth of each seismic trace until the objective function converges, and outputting the updated layer depth of each seismic trace. The operating steps of the present invention are simple, and can take into account the continuity of the seismic layer within the entire work area and the consistency with the well calibration layer information, and the tracking results are stable, fit the real layer trend, and have high accuracy.

Description

Technical Field

[0001] The present invention relates to the technical field of petroleum geophysical exploration, and in particular to a layer automatic tracking method, a layer automatic tracking device and a readable storage medium. Background Art

[0002] Identifying and tracking horizons in seismic data is a crucial step in seismic exploration. Structural interpretation is fundamental to seismic interpretation and also provides a crucial data basis for seismic data processing. Take, for example, depth-domain modeling of seismic data. Depth-domain modeling often employs a layer-by-layer approach, modeling from shallow to deep layers. Horizon interpretation is the first and most important step in depth-domain modeling. High-precision horizon interpretation significantly improves the accuracy and efficiency of velocity modeling.

[0003] Currently, two common horizon interpretation methods are manual picking and automatic tracking. Manual picking involves manually analyzing the continuity of events along a seismic profile and then manually picking horizons along the crests (or troughs, or zero-amplitude lines) of the seismic wavelets. Its advantages are that the interpretation results meet user expectations and are relatively intuitive. However, to achieve high-precision picking results, the user must perform relatively precise picking, which results in a significant workload for manual interpretation and horizon adjustment. Furthermore, the interpretation results rely heavily on the interpreter's geological experience. As seismic exploration areas expand, the volume of seismic data also increases. Relying solely on manual picking is not only labor-intensive but also inefficient in producing high-precision interpretation results. Automatic tracking interpretation involves tracking horizons based on the similarity of adjacent seismic traces after a seed point is given. Specifically, a seed point is manually assigned. At the sampling location of the seed point, the horizon waveform is extracted within a specified time window. Starting from the seed point, a search is conducted outward toward adjacent seismic traces for a waveform segment similar to the seed point waveform. The midpoint of this waveform segment is the depth of the horizon to be tracked in the adjacent traces. This point can be used as a new seed point to search for the waveform of the next adjacent channel. This method is convenient and easy to use, but the tracking results are easily affected by the seed point and the signal-to-noise ratio of the data. Tracking results using different seed points at the same location can vary. Cross-correlation and false faults can occur far from the initial artificial seed point, resulting in an incomplete layer plane. The tracking range is uncontrolled, and manual correction is impossible during the tracking process.

[0004] Although a grid-based horizon tracking method has been proposed, it involves gridding seismic data and automatically tracking within the grid. The gridded tracking results are then quality-controlled and modified, and reverse grid-corrected tracking is performed before continuing forward. However, this tracking method is not only complex but also fails to ensure the continuity of seismic horizons across the entire work area and their consistency with well calibration horizon information. Furthermore, the tracking results are often unstable and often lack accuracy that matches the true horizon direction. Summary of the Invention

[0005] The purpose of the embodiments of the present invention is to provide a method, device and readable storage medium for automatic layer tracking, so as to at least solve the problem that the above-mentioned existing layer tracking method is not only complicated in operation steps, but also cannot take into account the continuity of the seismic layer in the entire work area and the consistency with the well calibration layer information, and the tracking results are often not stable enough, and the overall accuracy problem of not being consistent with the actual layer trend often occurs.

[0006] In order to achieve the above object, the present invention provides a first aspect of a layer automatic tracking method, the method comprising:

[0007] Obtain seismic data and the depth of the same layer in different input wells;

[0008] Determining the dip angle of the stratigraphic section of each seismic trace based on the seismic data;

[0009] Determine the initial depth of each seismic trace based on the depth of the same horizon in different input wells;

[0010] Based on the dip angle of the stratigraphic section of each seismic trace and the initial layer depth of each seismic trace, the objective function is constructed with the depth of the same layer in different input wells as the constraint condition;

[0011] The objective function is iteratively calculated, and the layer depth of each seismic trace is iteratively updated until the objective function converges, and the updated layer depth of each seismic trace is output.

[0012] Optionally, determining the formation dip based on the seismic data includes:

[0013] Based on the seismic data, obtaining seismic profile data;

[0014] The seismic profile data is scanned for seismic dip angles to determine the stratum profile dip angle of each seismic trace in the seismic profile.

[0015] Optionally, the method further includes:

[0016] The seismic profile data is subjected to noise removal and amplitude equalization processing.

[0017] Optionally, the input well includes a first input well and a second input well;

[0018] Based on the depth of the same layer of the input well, determine the initial layer depth of each seismic trace, including:

[0019]

[0020] in, For the i The initial horizon depth of the seismic trace; is the depth of the horizon at the well location of the first input well; is the depth of the horizon at the well location of the second input well; N is the number of seismic traces between the first input well and the second input well.

[0021] Optionally, the objective function is expressed as:

[0022]

[0023]

[0024]

[0025] in, is the objective function; is the formation dip information matching function; N is the number of seismic traces between the first input well and the second input well; is the weight factor of the formation dip information matching function; For the i The initial horizon depth of the seismic trace; For the i -1 initial horizon depth of seismic trace; For the i The seismic trace and i -The lateral distance between seismic traces; The first i The dip angle of the stratigraphic profile of the seismic trace; is the damping function; is the weight factor of the damping function; For the i +1 initial horizon depth of seismic trace.

[0026] Optionally, the constraints of the objective function are:

[0027]

[0028] in, is the depth of the horizon at the well location of the first input well; is the depth of the layer at the well location of the second input well.

[0029] Optionally, iteratively calculating the objective function and iteratively updating the layer depth of each seismic trace until the objective function converges, and outputting the updated layer depth of each seismic trace, including:

[0030] The objective function is solved by using the least square method, and the layer depth of each seismic trace is updated in each iteration;

[0031] When the objective function converges, the layer depth of each seismic trace updated in this iteration is output.

[0032] Optionally, update the horizon depth of each seismic trace using the following formula:

[0033]

[0034] in, For the m After +1 iteration i The depth of the horizon of the seismic trace; For the m After the iteration i The depth of the horizon of the seismic trace; For the m The depth correction amount after iterations; is the depth of the horizon at the well location of the first input well; is the depth of the layer at the well location of the second input well.

[0035] A second aspect of the present invention provides a layer automatic tracking device, the device comprising:

[0036] Data acquisition module, used to obtain seismic data and the depth of the same layer in different input wells;

[0037] a stratigraphic section dip determination module, configured to determine the stratigraphic section dip of each seismic trace based on the seismic data;

[0038] An initial horizon depth determination module is used to determine the initial horizon depth of each seismic trace based on the depths of the same horizon in different input wells;

[0039] An objective function construction module is used to construct an objective function based on the stratigraphic profile dip of each seismic trace and the initial layer depth of each seismic trace, with the depth of the same layer in different input wells as a constraint condition;

[0040] The layer depth output module is used to iteratively calculate the objective function and iteratively update the layer depth of each seismic trace until the objective function converges, and output the updated layer depth of each seismic trace.

[0041] On the other hand, the present invention further provides a readable storage medium having stored thereon instructions for causing a machine to execute the above-mentioned automatic layer tracking method.

[0042] This technical solution has simple operating steps and can take into account the continuity of seismic layers throughout the entire work area and the consistency with the well calibration layer information. The tracking results are stable and fit the actual layer trend, achieving efficient and high-precision automatic layer tracking.

[0043] Other features and advantages of the embodiments of the present invention will be described in detail in the subsequent detailed description. BRIEF DESCRIPTION OF THE DRAWINGS

[0044] The accompanying drawings are used to provide a further understanding of the embodiments of the present invention and constitute a part of the specification. Together with the following detailed description, they are used to explain the embodiments of the present invention, but do not constitute a limitation of the embodiments of the present invention. In the accompanying drawings:

[0045] Figure 1 It is a flow chart of the automatic layer tracking method provided by the present invention;

[0046] Figure 2 is a schematic diagram of a seismic profile provided by the present invention;

[0047] Figure 3 is a schematic diagram of the dip angle of the stratum section provided by the present invention;

[0048] Figure 4 It is a schematic diagram of the initial layer depth provided by the present invention;

[0049] Figure 5 is a schematic diagram of the updated horizon depth provided by the present invention;

[0050] Figure 6 It is a schematic diagram of the convergence of the objective function provided by the present invention;

[0051] Figure 7 It is a structural schematic diagram of the layer automatic tracking device provided by the present invention.

[0052] Description of Reference Numerals

[0053] 10-data acquisition module; 20-stratum profile dip determination module;

[0054] 30-initial layer depth determination module; 40-objective function construction module;

[0055] 50-bit depth output module. DETAILED DESCRIPTION

[0056] The following describes the specific embodiments of the present invention in detail with reference to the accompanying drawings. It should be understood that the specific embodiments described herein are only used to illustrate and explain the present invention and are not intended to limit the present invention.

[0057] Figure 1 It is a flow chart of the automatic layer tracking method provided by the present invention; Figure 2 is a schematic diagram of a seismic profile provided by the present invention; Figure 3 is a schematic diagram of the dip angle of the stratum section provided by the present invention; Figure 4 It is a schematic diagram of the initial layer depth provided by the present invention; Figure 5 is a schematic diagram of the updated horizon depth provided by the present invention; Figure 6 It is a schematic diagram of the convergence of the objective function provided by the present invention; Figure 7 It is a structural schematic diagram of the layer automatic tracking device provided by the present invention.

[0058] like Figure 1 As shown, an embodiment of the present invention provides a layer automatic tracking method, the method comprising:

[0059] Step 1: Obtain seismic data and the depth of the same layer in different input wells;

[0060] In this example, seismic waves are excited to the surface and recorded, and the seismic data collected in the field are processed according to a conventional seismic data processing procedure to obtain seismic data.

[0061] Step 2: determining the dip angle of the stratigraphic profile of each seismic trace based on the seismic data;

[0062] Specifically, in this embodiment, determining the formation dip based on the seismic data includes:

[0063] Step 201: obtaining seismic profile data based on the seismic data;

[0064] Step 202: Scan the seismic profile data for seismic dip angles to determine the dip angle of each seismic trace in the seismic profile and record it as .

[0065] In another embodiment, after obtaining the seismic profile data in step 201, the method further includes:

[0066] The seismic profile data is subjected to noise removal and amplitude equalization processing.

[0067] By adopting the above method, the interference of noise signals can be eliminated, the accuracy of data can be guaranteed, and the stratigraphic profile dip of each seismic trace obtained after seismic dip scanning can be more accurate, thereby improving the accuracy of layer tracking.

[0068] Step 3: Determine the initial depth of each seismic trace based on the depths of the same layer in different input wells;

[0069] In this embodiment, the depth information of the same layer at different input well locations is obtained based on the calibrated layer information of the input well. Taking two wells (the first input well and the second input well) as an example, assuming that the depth of a certain layer at the first input well is , the depths at the first input well are ,based on and The initial value of the seismic horizon depth is constructed using linear interpolation. Specifically, the following calculation formula is used to calculate the initial horizon depth of each seismic trace:

[0070]

[0071] in, For the i The initial horizon depth of the seismic trace; is the depth of the horizon at the well location of the first input well; is the depth of the horizon at the well location of the second input well; N is the number of seismic traces between the first input well and the second input well; i is the index of the seismic trace.

[0072] Step 4: Based on the dip angle of the stratigraphic profile of each seismic trace and the initial layer depth of each seismic trace, the objective function is constructed with the depth of the same layer in different input wells as a constraint condition;

[0073] In this embodiment, the depth of the same layer in different input wells is the well information matching information. Then, an objective function based on the well information matching information and least squares optimization can be constructed. The expression of the objective function is:

[0074]

[0075]

[0076]

[0077] in, is the objective function; is the formation dip information matching function; N is the number of seismic traces between the first input well and the second input well; is the weight factor of the formation dip information matching function; For the i The initial horizon depth of the seismic trace; For the i -1 initial horizon depth of seismic trace; For the i The seismic trace and i -The lateral distance between seismic traces; The first i The dip angle of the stratigraphic profile of the seismic trace; is the damping function; is the weight factor of the damping function; For the i +1 initial horizon depth of seismic trace.

[0078] More specifically, It can be expressed as , which means that according to i Horizon depth of seismic traces Hedi i - Horizon depth of 1 seismic trace The calculated i The calculated dip angle of the channel is consistent with the i The dip angles of the stratigraphic sections of the seismic traces are compared iteratively; similarly, It can be expressed as , which means that according to i +1 seismic trace depth Hedi i Horizon depth of seismic traces The calculated i The calculated dip angle of the +1 channel is the same as that of the first channel in the seismic profile. i The dip angles of the stratigraphic sections of the seismic traces are iteratively compared.

[0079] The constraints of the objective function are:

[0080]

[0081] in, is the depth of the horizon at the well location of the first input well; is the depth of the layer at the well location of the second input well.

[0082] In this embodiment, the set constraint conditions can ensure that the depth of the same layer at the well position of the first input well and the depth at the well position of the second input well remain unchanged after the update.

[0083] Step 5: Iteratively calculate the objective function and iteratively update the layer depth of each seismic trace until the objective function converges, and then output the updated layer depth of each seismic trace.

[0084] Specifically, the method includes: solving the objective function by using the least square method, and updating the layer depth of each seismic trace in each iteration;

[0085] When the objective function converges, the layer depth of each seismic trace updated during this iteration is output, and the updated layer depth of each seismic trace is obtained. Then, based on the layer depth of each seismic trace, the layer in the entire stratum can be obtained.

[0086] In this embodiment, the least squares iterative optimization of the constructed initial layer depth is performed using the objective function F. The first-order derivative of the objective function F with respect to the layer depth is calculated to construct the equation:

[0087]

[0088] in, For the i The initial horizon depth of the seismic trace; i is 0, 1, ..., N -1;

[0089] The equation constructed by the first-order derivative is a nonlinear equation and can be solved by the Newton iteration method. First, use the initial value of the layer depth obtained above as Then, in each iteration m, the layer depth correction amount Accumulated to the current layer depth Get the updated layer depth:

[0090]

[0091] in, For the m After +1 iteration i The depth of the horizon of the seismic trace; For the m After the iteration i The depth of the horizon of the seismic trace; For the m The depth correction amount after iterations; is the depth of the horizon at the well location of the first input well; is the depth of the layer at the well location of the second input well.

[0092] Specifically, This can be obtained by solving the following five diagonal equations:

[0093]

[0094] in, is the objective function F The second derivative with respect to the depth of the horizon, j is the index of the seismic trace. And, j The value of is: , therefore, through able to obtain The value of .

[0095] More specifically, in this embodiment, in order to demonstrate the technical advantages of the present invention, an effect test was performed using complex seismic data. Figure 2 The image shows a seismic profile after noise removal and amplitude equalization. The horizontal axis represents the seismic trace number, while the vertical axis represents the imaging depth. As can be seen, the data is complex, with layers with steep dips, making it quite challenging to analyze. Figure 3 Yes Figure 2 The stratigraphic dip angle information obtained by scanning the seismic profile shown is the stratigraphic dip angle of each seismic trace, where the horizontal and vertical coordinate axes have the same meaning as Figure 2 , the figure shows the tangent value of the dip angle of the seismic section. Figure 4 It is the initial value of the seismic horizon depth constructed after the horizon depths of two wells are given. It is obtained by linear interpolation based on the horizon depths of the two wells. The horizontal and vertical coordinate axes have the same meaning. Figure 2 The layer depth obtained by using the automatic tracking algorithm of the present invention is as follows: Figure 5 The black solid line shows the very high layer tracking accuracy. The black solid line almost completely coincides with the crest of the seismic event. Even in steep structures, the tracked layers still closely match the event, verifying the effectiveness and high accuracy of the algorithm. Figure 6 The figure shows the convergence of the least squares iteration objective function. The horizontal axis represents the number of iterations, and the vertical axis represents the objective function value. As can be seen, after seven iterations, the objective function value met the requirements. The run time was approximately 0.5 seconds, meeting the requirements for real-time interactive layer picking. This method achieves efficient and high-precision automatic layer tracking, overcoming the challenges of poor stability and low accuracy associated with traditional methods.

[0096] like Figure 7 As shown, an embodiment of the present invention further provides a layer automatic tracking device, the device comprising:

[0097] The data acquisition module 10 is used to obtain seismic data and the depth of the same layer in different input wells;

[0098] a stratigraphic section dip determination module 20 for determining the stratigraphic section dip of each seismic trace based on the seismic data;

[0099] An initial horizon depth determination module 30 is used to determine the initial horizon depth of each seismic trace based on the depths of the same horizon in different input wells;

[0100] An objective function construction module 40 is used to construct an objective function based on the stratigraphic profile dip of each seismic trace and the initial layer depth of each seismic trace, with the depth of the same layer in different input wells as a constraint condition;

[0101] The horizon depth output module 50 is used to iteratively calculate the objective function and iteratively update the horizon depth of each seismic trace until the objective function converges, and then output the updated horizon depth of each seismic trace.

[0102] An embodiment of the present invention further provides a readable storage medium having instructions stored thereon, the instructions being used to enable a machine to execute the above-mentioned automatic layer tracking method.

[0103] Those skilled in the art will appreciate that all or part of the steps in the methods described in the aforementioned embodiments can be performed by instructing the relevant hardware through a program. The program, stored in a storage medium, includes instructions for causing a microcontroller, chip, or processor to execute all or part of the steps in 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 mobile hard drive, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk.

[0104] Those skilled in the art can clearly understand that, for the convenience and simplicity of description, only the division of the above-mentioned functional units and modules is used as an example for illustration. In actual applications, the above-mentioned functions can be distributed and completed by different functional units and modules as needed, that is, the internal structure of the device can be divided into different functional units or modules to complete all or part of the functions described above. The functional units and modules in the embodiment can be integrated into one processing unit, or each unit can exist physically alone, or two or more units can be integrated into one unit. The above-mentioned integrated unit can be implemented in the form of hardware or in the form of a software functional unit. In addition, the specific names of the functional units and modules are only for the convenience of distinguishing each other and are not used to limit the scope of protection of this application.

[0105] The above describes in detail the optional embodiments of the present invention in conjunction with the accompanying drawings. However, the embodiments of the present invention are not limited to the specific details in the above embodiments. Within the technical concept of the embodiments of the present invention, a variety of simple modifications can be made to the technical solutions of the embodiments of the present invention, and these simple modifications all fall within the scope of protection 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 unless there is any contradiction. In order to avoid unnecessary repetition, the embodiments of the present invention will no longer describe the various possible combinations separately.

[0106] 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 embodiments of the present invention, they should also be regarded as the contents disclosed in the embodiments of the present invention.

Claims

1. A method for automatic layer tracking, characterized in that: The method comprises: Obtain seismic data and the depth of the same layer in different input wells; Determining the dip angle of the stratigraphic section of each seismic trace based on the seismic data; Determine the initial depth of each seismic trace based on the depth of the same horizon in different input wells; Based on the dip angle of the stratigraphic section of each seismic trace and the initial layer depth of each seismic trace, the objective function is constructed with the depth of the same layer in different input wells as the constraint condition; Iteratively calculating the objective function and iteratively updating the layer depth of each seismic trace until the objective function converges, and outputting the updated layer depth of each seismic trace; Wherein, the expression of the objective function is: in, is the objective function; is the formation dip information matching function; N is the number of seismic traces between the first input well and the second input well; is the weight factor of the formation dip information matching function; For the i +1 initial horizon depth of seismic trace; For the i The initial horizon depth of the seismic trace; For the i +1 seismic channel and i The lateral distance between seismic traces; The first i +1 stratigraphic profile dip angle of seismic trace; is the damping function; is the weight factor of the damping function; For the i + Initial horizon depth of 2 seismic traces.

2. The automatic layer tracking method according to claim 1, characterized in that: Determine the dip angle of the stratigraphic section of each seismic trace based on the seismic data, including: Based on the seismic data, obtaining seismic profile data; The seismic profile data is scanned for seismic dip angles to determine the stratum profile dip angle of each seismic trace in the seismic profile.

3. The automatic layer tracking method according to claim 2, characterized in that: The method further comprises: The seismic profile data is subjected to noise removal and amplitude equalization processing.

4. The automatic layer tracking method according to claim 1, characterized in that: The input wells include a first input well and a second input well; Determine the initial depth of each seismic trace based on the depth of the same horizon in different input wells, including: in, For the i +1 initial horizon depth of seismic trace; is the depth of the horizon at the well location of the first input well; is the depth of the horizon at the well location of the second input well; N is the number of seismic traces between the first input well and the second input well.

5. The automatic layer tracking method according to claim 4, characterized in that: The constraints of the objective function are: in, is the depth of the horizon at the well location of the first input well; is the depth of the horizon at the well location of the second input well; is the initial horizon depth of the first seismic trace; For the N The initial horizon depth of the seismic trace.

6. The automatic layer tracking method according to claim 1, characterized in that: Iteratively calculating the objective function and iteratively updating the layer depth of each seismic trace until the objective function converges, and outputting the updated layer depth of each seismic trace, including: The objective function is solved by using the least square method, and the layer depth of each seismic trace is updated in each iteration; When the objective function converges, the layer depth of each seismic trace updated in this iteration is output.

7. The automatic layer tracking method according to claim 1, characterized in that: The following formula is used to update the horizon depth of each seismic trace: in, For the m After +1 iteration i +1 seismic trace depth; For the m After the iteration i +1 seismic trace depth; For the m The depth correction amount after iterations; is the depth of the horizon at the well location of the first input well; is the depth of the layer at the well location of the second input well.

8. A layer automatic tracking device, characterized in that: The device comprises: Data acquisition module, used to obtain seismic data and the depth of the same layer in different input wells; a stratigraphic section dip determination module, configured to determine the stratigraphic section dip of each seismic trace based on the seismic data; An initial horizon depth determination module is used to determine the initial horizon depth of each seismic trace based on the depths of the same horizon in different input wells; An objective function construction module is used to construct an objective function based on the stratigraphic profile dip of each seismic trace and the initial layer depth of each seismic trace, with the depth of the same layer in different input wells as a constraint condition; A horizon depth output module is used to iteratively calculate the objective function and iteratively update the horizon depth of each seismic trace until the objective function converges, and output the updated horizon depth of each seismic trace; Wherein, the expression of the objective function is: in, is the objective function; is the formation dip information matching function; N is the number of seismic traces between the first input well and the second input well; is the weight factor of the formation dip information matching function; For the i +1 initial horizon depth of seismic trace; For the i The initial horizon depth of the seismic trace; For the i +1 seismic channel and i The lateral distance between seismic traces; The first i +1 stratigraphic profile dip angle of seismic trace; is the damping function; is the weight factor of the damping function; For the i + Initial horizon depth of 2 seismic traces.

9. A readable storage medium having instructions stored thereon, wherein the instructions are used to enable a machine to execute the automatic layer tracking method according to any one of claims 1 to 7.

Images