Fracture-controlled beamlet tomography regularization method, apparatus, computer and storage medium

Abstract

The application provides a fracture-controlled beam tomography regularization method, device, computer and storage medium, which comprises the following steps: acquiring a seismic profile; extracting a coherent attribute according to the seismic profile; resampling the coherent attribute to obtain a resampled coherent attribute corresponding to a tomography grid; obtaining a ray density through ray tracing; calculating a fracture-controlled beam tomography operator according to the resampled coherent attribute and the ray density; introducing the fracture-controlled beam tomography operator into a tomography inversion target function for regularization constraint to obtain the tomography inversion target function after regularization constraint. The fracture-controlled beam tomography operator controls the direction of a ray beam to avoid an abnormal body, avoids unmodelable scattering, and ensures the accuracy of a ray tracing path. Meanwhile, the fracture-controlled beam tomography operator monitors an abnormal ray beam path through data driving, further eliminates a multi-path ray beam caused by a velocity anomaly, and improves the stability and precision of tomography.

Description

Technical Field

[0001] This invention relates to the field of geophysical exploration technology, and in particular to a fracture control bundle tomography regularization method, apparatus, computer, and storage medium. Background Technology

[0002] Currently, depth-domain velocity modeling technology is widely used and mature in both commercial software and theoretical research. Among these, migration velocity modeling, also known as tomographic inversion, remains the most widely applied and technologically mature. Tomographic inversion typically relies on ray-based algorithms, dividing the velocity model into numerous regular grids and then using ray tracing to invert the velocity within each grid. This technique can invert large-scale velocity models, primarily producing low-wavenumber components of the velocity. For characterizing large-scale structures, its inversion accuracy is acceptable. However, with increasing exploration depth and shrinking exploration targets, conventional ray-based tomography is no longer sufficient. Therefore, Gaussian beam tomography, which is closer to wave-based algorithms, and wave-based full-waveform inversion techniques have emerged. However, full-waveform inversion currently has low computational efficiency and limited effectiveness in practical data applications, with a long theoretical road ahead. Therefore, its application in real-world data is limited. Thus, Gaussian beam tomography, a technique somewhere between the two, performs well in practical data applications and is being widely adopted. Technically speaking, Gaussian beam tomography is currently the most efficient and accurate velocity modeling technique. However, this technique lacks specificity and has a lot of room for optimization. In particular, for complex structures such as fractures, the presence of anomalies often leads to many abnormal or erroneous beams in Gaussian beam tracing, affecting velocity inversion and limiting modeling accuracy.

[0003] Furthermore, from an exploration perspective, investigating the structural characteristics of underground faults has become a crucial aspect of seismic exploration. Achieving high-precision velocity characterization on both sides of a fault is of great significance for high-precision imaging. Currently, depth-domain velocity modeling techniques for faults mainly involve introducing fault information for constraint, falling into two main categories: one is to introduce manually interpreted fault information for velocity modeling constraints. This method provides strong constraints on the velocity model but carries the risk of errors in fault interpretation. If misinterpreted, the modeling results will also develop in the wrong direction, and faults are difficult to interpret completely correctly manually; therefore, this technique has not been widely adopted. The other category involves iterative solutions in migration velocity modeling to guide the modeling in the correct direction. In each iteration, a purely data-driven algorithm extracts fault features for fault-controlled constraint inversion. This technique is relatively more widely used, but due to the presence of erroneous rays or ray tracing results with large errors during ray tracing, even with data-driven fault-controlled constraints, significant errors are unavoidable. In conclusion, current fault modeling results do not fundamentally solve the problem; they only improve the results through constraints. Therefore, there is still considerable room for improvement in the modeling results. Summary of the Invention

[0004] Therefore, it is necessary to provide a method, apparatus, computer, and storage medium for fracture control bundle tomography regularization to address the aforementioned technical problems.

[0005] A fracture control bundle tomography regularization method includes:

[0006] Obtain seismic profiles;

[0007] Extract coherent properties from the seismic profile;

[0008] The resampled coherent properties are obtained by resampling the coherent properties to obtain the resampled coherent properties corresponding to the tomographic grid;

[0009] The ray density is obtained through ray tracing;

[0010] The crack control bundle tomography operator is calculated based on the resampling coherence properties and the ray density.

[0011] The crack control bundle tomography operator is introduced into the tomography inversion objective function for regularization constraint to obtain the tomography inversion objective function after regularization constraint.

[0012] In one embodiment, the step of extracting coherence attributes from the seismic profile includes:

[0013] Based on the seismic profile, the imaging results of each wheel in the migration velocity modeling are obtained;

[0014] The coherence properties are extracted based on the previous imaging results in the offset velocity modeling.

[0015] In one embodiment, in the step of resampling the coherent properties to obtain the resampled coherent properties corresponding to the tomographic grid, the resampling calculation formula is:

[0016]

[0017] in This indicates that the position after resampling to the tomographic grid size is a point. Relevance properties This indicates that the location point was calculated directly using the results of the previous imaging round. Relevance properties This represents the resampling operator.

[0018] In one embodiment, the step of calculating the crack control beam tomography operator based on the resampling coherence property and the ray density includes:

[0019] Based on the resampling coherence properties and the ray density The fracture control bundle tomography operator is calculated using the following formula:

[0020]

[0021] in, The grid points are The location of the crack control bundle tomography operator; Indicates that it is located at grid point The calculation time window is The resampling coherence property, Indicates that it is located at grid point Time window The ray density value, , and These are the half-time window lengths in three directions.

[0022] In one embodiment, the step of introducing the crack control bundle tomography operator into the tomography inversion objective function for regularization constraints is implemented by using a preconditional regularization operator, calculated as follows:

[0023]

[0024] in, For the tomographic inversion linearization operator, For the crack control bundle tomography operator, The solution is given by the preconditions.

[0025] In one embodiment, after the step of obtaining the tomographic inversion objective function with regularization constraints, the method further includes:

[0026] Solving the tomographic inversion objective function after regularization constraints yields the updated velocity around the fracture.

[0027] In one embodiment, the step of solving the tomographic inversion objective function after regularization constraints to obtain the updated velocity around the fracture includes:

[0028] The tomographic inversion objective function after regularization constraints is solved using the preconditional conjugate gradient method to obtain the velocity update around the fracture.

[0029] A fracture-controlled bundle tomography regularization device, comprising:

[0030] The seismic profile acquisition module is used to acquire seismic profiles.

[0031] The coherence attribute extraction module is used to extract coherence attributes based on the seismic profile;

[0032] The resampling module is used to resample the coherent properties to obtain the resampled coherent properties corresponding to the tomographic grid;

[0033] A ray density acquisition module is used to obtain ray density through ray tracing.

[0034] The control beam tomography operator calculation module is used to calculate the crack control beam tomography operator based on the resampling coherence properties and the ray density;

[0035] The regularization module is used to introduce the fracture control bundle tomography operator into the tomography inversion objective function for regularization constraints, so as to obtain the tomography inversion objective function after regularization constraints.

[0036] A computer device includes a memory and a processor, the memory storing a computer program, characterized in that the processor executes the computer program to perform the following steps:

[0037] Obtain seismic profiles;

[0038] Extract coherent properties from the seismic profile;

[0039] The resampled coherent properties are obtained by resampling the coherent properties to obtain the resampled coherent properties corresponding to the tomographic grid;

[0040] The ray density is obtained through ray tracing;

[0041] The crack control bundle tomography operator is calculated based on the resampling coherence properties and the ray density.

[0042] The crack control bundle tomography operator is introduced into the tomography inversion objective function for regularization constraint to obtain the tomography inversion objective function after regularization constraint.

[0043] A computer-readable storage medium having a computer program stored thereon, the computer program performing the following steps when executed by a processor:

[0044] Obtain seismic profiles;

[0045] Extract coherent properties from the seismic profile;

[0046] The resampled coherent properties are obtained by resampling the coherent properties to obtain the resampled coherent properties corresponding to the tomographic grid;

[0047] The ray density is obtained through ray tracing;

[0048] The crack control bundle tomography operator is calculated based on the resampling coherence properties and the ray density.

[0049] The crack control bundle tomography operator is introduced into the tomography inversion objective function for regularization constraint to obtain the tomography inversion objective function after regularization constraint.

[0050] The aforementioned crack-controlled beam tomography regularization method, apparatus, computer, and storage medium control the direction of the ray beam to avoid anomalies through the crack-controlled beam tomography operator, thereby avoiding unsimulated scattering and ensuring the accuracy of the ray tracing path. At the same time, by monitoring abnormal ray beam paths through data-driven monitoring, multi-path (erroneous path) ray beams caused by velocity anomalies are further eliminated to improve tomographic stability and accuracy. Attached Figure Description

[0051] Figure 1 This is a flowchart illustrating the crack control bundle tomography regularization method in one embodiment;

[0052] Figure 2 A structural block diagram of a fracture control bundle tomography regularization device in one embodiment;

[0053] Figure 3 This is an internal structural diagram of a computer device in one embodiment;

[0054] Figure 4A This is a schematic diagram of coherence properties calculated using seismic profiles in one embodiment.

[0055] Figure 4B This is a schematic diagram of the resampled coherence properties obtained after resampling in one embodiment;

[0056] Figure 4C This is a schematic diagram of the ray density obtained by ray tracing in one embodiment;

[0057] Figure 4D This is a schematic diagram of the ray density obtained by ray beam tracing after constraining with the crack control bundle tomography regularization operator in one embodiment.

[0058] Figure 4E This is a schematic diagram of a velocity model obtained by conventional tomographic inversion in one embodiment;

[0059] Figure 4F This is a schematic diagram of the velocity model obtained by tomographic inversion after adding the crack control bundle tomographic regularization operator in one embodiment;

[0060] Figure 4G This is a schematic diagram of the imaging results corresponding to a conventional tomographic inversion velocity model in one embodiment;

[0061] Figure 4H This is a schematic diagram of the imaging results corresponding to the velocity model of tomographic inversion after adding the crack control bundle tomographic regularization operator in one embodiment. Detailed Implementation

[0062] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the scope of this application.

[0063] Example 1

[0064] In this embodiment, as Figure 1 As shown, a fracture control bundle tomography regularization method is provided, which includes:

[0065] Step 110: Obtain the seismic profile.

[0066] Step 120: Extract coherent attributes based on the seismic profile.

[0067] In this embodiment, coherent attributes are extracted using seismic profiles. .

[0068] In one embodiment, the step of extracting coherence attributes from the seismic profile includes: obtaining the imaging results of each round in the migration velocity modeling based on the seismic profile; and extracting the coherence attributes based on the imaging results of the previous round in the migration velocity modeling. .

[0069] Step 130: Resample the coherent properties to obtain the resampled coherent properties corresponding to the tomographic grid.

[0070] In one embodiment, in the step of resampling the coherent properties to obtain the resampled coherent properties corresponding to the tomographic grid, the resampling calculation formula is:

[0071] (1)

[0072] in This indicates that the position after resampling to the tomographic grid size is a point. Relevance properties This indicates that the location point was calculated directly using the results of the previous imaging round. The coherence property, which can be calculated using the coherence property formula, is... This represents the resampling operator.

[0073] In this embodiment, the average resampling operator is used for resampling calculation. For example, to... The three dimensions are resampled at intervals of 1 to form a tomographic grid. Therefore, the formula for resampling can be written as:

[0074] (2)

[0075] in, , , , This is the floor operator.

[0076] Step 140: Obtain the ray density through ray tracing.

[0077] In this embodiment, after obtaining the coherent properties after resampling, the ray density needs to be obtained through the previous round of ray tracing. .

[0078] Step 150: Calculate the crack control bundle tomography operator based on the resampling coherence properties and the ray density.

[0079] In one embodiment, the step of calculating the crack control beam tomography operator based on the resampling coherence property and the ray density includes:

[0080] Based on the resampling coherence properties and the ray density The fracture control bundle tomography operator is calculated using the following formula:

[0081] (3)

[0082] in, The grid points are The location of the crack control bundle tomography operator; Indicates that it is located at grid point The calculation time window is The resampling coherence property, Indicates that it is located at grid point Time window The ray density value, , and These are the half-time window lengths in three directions.

[0083] Step 160: Introduce the crack control bundle tomography operator into the tomography inversion objective function for regularization constraint to obtain the tomography inversion objective function after regularization constraint.

[0084] In this embodiment, by introducing the operator F into the tomography inversion matrix, the control bundle tomography velocity modeling technique for fractures can be realized.

[0085] In one embodiment, the step of introducing the crack control bundle tomography operator into the tomography inversion objective function for regularization constraints is implemented by using a preconditional regularization operator, calculated as follows:

[0086] (4)

[0087] in, For the tomographic inversion linearization operator, The fracture control bundle tomography operator designed in equation (3) is used for calculation. The solution is given by the preconditions.

[0088] In one embodiment, after obtaining the tomographic inversion objective function after regularization constraints, the method further includes: solving the tomographic inversion objective function after regularization constraints to obtain the velocity update amount around the fracture.

[0089] In this embodiment, a high-precision velocity model for the fracture can be obtained by solving formula (4), which can significantly improve the velocity modeling accuracy on both sides of the fracture.

[0090] In one embodiment, the step of solving the tomographic inversion objective function after regularization constraints to obtain the velocity update amount around the fracture includes: solving the tomographic inversion objective function after regularization constraints using the preconditional conjugate gradient method to obtain the velocity update amount around the fracture.

[0091] In this embodiment, the calculation formula for the preconditional conjugate gradient method is as follows:

[0092] (5)

[0093] In the above embodiments, the operator controls the direction of the X-ray beam to avoid anomalous objects, avoids scattering that cannot be simulated, and ensures the accuracy of the X-ray tracing path; at the same time, by monitoring anomalous X-ray beam paths through data-driven monitoring, multi-path (erroneous path) X-ray beams caused by velocity anomalies are further eliminated to improve the stability and accuracy of tomography.

[0094] Example 2

[0095] This application presents a technology developed to address the challenge of fracture velocity modeling. The method employed involves using coherence properties to guide and adjust the ray beam distribution, forming a control beam tomography regularization operator introduced into the objective function. This regularization operator has two main effects: it controls the ray beam direction to avoid anomalous bodies, preventing unsimulable scattering and ensuring accurate ray tracing paths; simultaneously, it uses data-driven monitoring of anomalous ray beam paths to further eliminate multi-path (erroneous path) ray beams caused by velocity anomalies, thereby improving tomographic stability and accuracy.

[0096] The fracture control bundle tomography regularization method in this application is based on the migration velocity modeling process, with each update based on the previous result. First, coherent attributes are extracted from the previous imaging results in the migration velocity modeling, and then the attributes are resampled. The specific calculation formula is as follows:

[0097] (1)

[0098] in This indicates that the position after resampling to the tomographic grid size is a point. Relevance properties This indicates that the location point was calculated directly using the results of the previous imaging round. The coherent properties can be directly calculated using commercial software or coherent property calculation formulas. This patent suggests using the average resampling operator for resampling calculations, for example, to represent the resampling operator. The three dimensions are resampled at intervals of 1 to form a tomographic grid. Then the above formula can be written as

[0099] (2)

[0100] in , , , This is the rounding operator. After obtaining the resampled coherent properties, the ray density needs to be obtained from the previous round of ray tracing. Then, the fracture control bundle tomography regularization operator is calculated using the following formula:

[0101] (3)

[0102] in, The grid points are Location-based crack control bundle regularization operator. Indicates that it is located at grid point The calculation time window is The resampling coherence property, Indicates that it is located at grid point Time window The ray density value. , and These represent the half-time window lengths in three directions. Introducing operator F into the tomographic inversion matrix enables control-bundle tomographic velocity modeling for fractures. The introduction method employs a pre-conditional regularization operator, the specific calculation formula of which is as follows:

[0103] (4)

[0104] In the formula, For the tomographic inversion linearization operator, This refers to the preconditional regularization operator designed in formula (3). The key to this invention is the construction of this operator. The solution is given by the preconditions. Therefore, by solving equation (4), a high-precision velocity model for the fracture can be obtained, which can significantly improve the accuracy of velocity modeling on both sides of the fracture. The precondition conjugate gradient method can be used to solve the equation, and the specific calculation formula is as follows:

[0105] (5)

[0106] The benefits of this application: A technology developed to address the challenge of velocity modeling around fractures. The method employed is as follows: Coherence properties guide the adjustment of the ray beam distribution, forming a control beam tomography regularization operator introduced into the objective function. This regularization operator has two main effects: Firstly, it controls the ray beam direction to avoid anomalies, preventing unsimulable scattering and ensuring accurate ray tracing paths. Secondly, it uses data-driven monitoring of abnormal ray beam paths to further eliminate multi-path (erroneous path) ray beams caused by velocity anomalies, thereby improving tomographic stability and accuracy. The specific steps are as follows:

[0107] 1) Extracting coherence attributes from seismic profiles ;

[0108] 2) Resample the coherent properties to obtain the resampled coherent properties corresponding to the tomographic mesh. ;

[0109] 3) Obtaining ray density through ray tracing ;

[0110] 4) Using the resampling coherence properties and ray density obtained in steps 2-3, construct the fracture control bundle tomography operator using formula (3). This is then introduced into the torsion inversion objective function for regularization constraints;

[0111] 5) Solve the objective function to obtain the velocity update around the fracture. We completed the high-precision velocity modeling work for fracture.

[0112] To demonstrate the correctness and effectiveness of the method and to show that it has higher accuracy, the following explanation is based on actual data testing.

[0113] like Figure 4A As shown, this is the coherence property calculated using seismic profiles in the embodiment. It can be seen that the fracture morphology is clearly depicted and can be used for fracture guidance.

[0114] like Figure 4B As shown, this is the resampled coherence property obtained after resampling. It can be seen that the signal-to-noise ratio is improved after resampling, the fracture characterization is clearer, and it is more conducive to fracture guidance.

[0115] like Figure 4C As shown, this is the ray density obtained by ray tracing in the embodiment, which is used to construct the fracture control bundle tomography regularization operator. It can be seen that without the constraint of the fracture control bundle tomography regularization operator, the ray density looks rather messy overall, and there are many local anomalies, which is not conducive to high-precision velocity inversion.

[0116] like Figure 4D As shown, this is the ray density obtained by ray beam tracing after constraint using the fracture-controlled bundle tomography regularization operator. It can be seen that after adopting the fracture-controlled bundle tomography regularization technique, compared to... Figure 3 As a result, the beam density is more uniform, and the shadow of the guide fracture can be seen in the beam density.

[0117] like Figure 4E As shown, this is the velocity model obtained using conventional tomographic inversion velocity modeling technology in the embodiment. It can be seen that the large-scale structural morphology is depicted, but the velocity difference between the two sides of the fracture is not reflected.

[0118] like Figure 4F As shown, this is the velocity model obtained by tomographic inversion after adding the fracture control bundle tomographic regularization operator in the embodiment. It can be seen that after adding the fracture control bundle tomographic regularization operator, the velocity response on both sides of the fracture can be clearly seen in the velocity model, and the accuracy of the velocity model is significantly improved.

[0119] like Figure 4G As shown, this is Figure 4E The imaging profile corresponding to the conventional tomographic inversion velocity model shows that the fracture imaging accuracy is low and there is still considerable room for improvement.

[0120] like Figure 4H As shown, this is Figure 4FAfter adding the fracture control bundle tomography regularization operator, the velocity model obtained by tomography inversion shows that the imaging accuracy of the main fracture is significantly improved, the breakpoints are more brittle, and the fracture surface is clearer.

[0121] It should be understood that, although Figure 1 The steps in the flowchart are shown sequentially as indicated by the arrows, but these steps are not necessarily executed in the order indicated by the arrows. Unless otherwise specified herein, there is no strict order in which these steps are executed, and they can be performed in other orders. Figure 1 At least some of the steps in the process may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily completed at the same time, but can be executed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but can be executed in turn or alternately with other steps or at least some of the sub-steps or stages of other steps.

[0122] Example 3

[0123] In this embodiment, as Figure 2 As shown, a fracture control bundle tomography regularization device is provided, comprising:

[0124] Seismic profile acquisition module 210 is used to acquire seismic profiles;

[0125] The coherence attribute extraction module 220 is used to extract coherence attributes based on the seismic profile;

[0126] The resampling module 230 is used to resample the coherent properties to obtain the resampled coherent properties corresponding to the tomographic grid;

[0127] Module 240 for obtaining ray density via ray tracing;

[0128] Control beam tomography operator calculation module 250 is used to calculate the crack control beam tomography operator based on the resampled coherence properties and the ray density;

[0129] The regularization module 260 is used to introduce the fracture control bundle tomography operator into the tomography inversion objective function for regularization constraint, so as to obtain the tomography inversion objective function after regularization constraint.

[0130] In one embodiment, the seismic profile acquisition module includes:

[0131] The imaging result acquisition unit is used to acquire the imaging results of each wheel in the migration velocity modeling based on the seismic profile.

[0132] The coherent attribute extraction unit is used to extract the coherent attributes based on the previous imaging results in the offset velocity modeling.

[0133] In one embodiment, the resampling module is used to resample the coherence properties to obtain resampled coherence properties corresponding to the tomographic grid, wherein the resampling calculation formula is:

[0134]

[0135] in This indicates that the position after resampling to the tomographic grid size is a point. Relevance properties This indicates that the location point was calculated directly using the results of the previous imaging round. Relevance properties This represents the resampling operator.

[0136] In one embodiment, the control beam tomography operator calculation module is used for:

[0137] Based on the resampling coherence properties and the ray density The fracture control bundle tomography operator is calculated using the following formula:

[0138]

[0139] in, The grid points are The location of the crack control bundle tomography operator; Indicates that it is located at grid point The calculation time window is The resampling coherence property, Indicates that it is located at grid point Time window The ray density value, , and These are the half-time window lengths in three directions.

[0140] In one embodiment, the regularization module is used to introduce the fracture control bundle tomography operator into the tomography inversion objective function for regularization constraints. The reference method is to use a pre-conditional regularization operator, and the calculation formula is:

[0141]

[0142] in, For the tomographic inversion linearization operator, For the crack control bundle tomography operator, The solution is given by the preconditions.

[0143] In one embodiment, the fracture control bundle tomography regularization device further includes:

[0144] The function solving module is used to solve the tomographic inversion objective function after regularization constraints to obtain the velocity update around the fracture.

[0145] In one embodiment, the function solving module is further used to solve the tomographic inversion objective function after regularization constraints using the preconditioned conjugate gradient method to obtain the velocity update amount around the fracture.

[0146] Specific limitations regarding the fracture control bundle tomography regularization device can be found in the limitations of the fracture control bundle tomography regularization method above, and will not be repeated here. Each unit in the aforementioned fracture control bundle tomography regularization device can be implemented entirely or partially through software, hardware, or a combination thereof. These units can be embedded in or independent of the processor in a computer device in hardware form, or stored in the memory of a computer device in software form, so that the processor can call and execute the operations corresponding to each unit.

[0147] Example 4

[0148] In this embodiment, a computer device is provided. Its internal structure diagram can be shown as follows: Figure 3 As shown, the computer device includes a processor, memory, network interface, display screen, and input devices connected via a system bus. The processor provides computational and control capabilities. The memory includes a non-volatile storage medium and internal memory. The non-volatile storage medium stores an operating system and computer programs, and also deploys a database for seismic data. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage medium. The network interface is used to communicate with other computer devices that have deployed application software. When executed by the processor, the computer program implements a fracture control bundle tomography regularization method. The display screen can be a liquid crystal display (LCD) or an e-ink display. The input devices can be a touch layer covering the display screen, buttons, a trackball, or a touchpad on the computer device's casing, or an external keyboard, touchpad, or mouse.

[0149] Those skilled in the art will understand that Figure 3 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0150] In one embodiment, a computer device is provided, including a memory and a processor, the memory storing a computer program, the processor executing the computer program to perform the following steps:

[0151] Obtain seismic profiles;

[0152] Extract coherent properties from the seismic profile;

[0153] The resampled coherent properties are obtained by resampling the coherent properties to obtain the resampled coherent properties corresponding to the tomographic grid;

[0154] The ray density is obtained through ray tracing;

[0155] The crack control bundle tomography operator is calculated based on the resampling coherence properties and the ray density.

[0156] The crack control bundle tomography operator is introduced into the tomography inversion objective function for regularization constraint to obtain the tomography inversion objective function after regularization constraint.

[0157] In one embodiment, the processor, when executing a computer program, also performs the following steps:

[0158] Based on the seismic profile, the imaging results of each wheel in the migration velocity modeling are obtained;

[0159] The coherence properties are extracted based on the previous imaging results in the offset velocity modeling.

[0160] In one embodiment, when the processor executes a computer program to implement the step of resampling the coherent attributes to obtain the resampled coherent attributes corresponding to the tomographic grid, the resampling calculation formula is:

[0161]

[0162] in This indicates that the position after resampling to the tomographic grid size is a point. Relevance properties This indicates that the location point was calculated directly using the results of the previous imaging round. Relevance properties This represents the resampling operator.

[0163] In one embodiment, the processor, when executing a computer program, also performs the following steps:

[0164] Based on the resampling coherence properties and the ray density The fracture control bundle tomography operator is calculated using the following formula:

[0165]

[0166] in, The grid points are The location of the crack control bundle tomography operator; Indicates that it is located at grid point The calculation time window is The resampling coherence property, Indicates that it is located at grid point Time window The ray density value, , and These are the half-time window lengths in three directions.

[0167] In one embodiment, when the processor executes the computer program, the step of introducing the crack control bundle tomography operator into the tomography inversion objective function for regularization constraints is implemented, wherein the reference method is to use a preconditional regularization operator, and the calculation formula is:

[0168]

[0169] in, For the tomographic inversion linearization operator, For the crack control bundle tomography operator, The solution is given by the preconditions.

[0170] In one embodiment, the processor, when executing a computer program, also performs the following steps:

[0171] Solving the tomographic inversion objective function after regularization constraints yields the updated velocity around the fracture.

[0172] In one embodiment, the processor, when executing a computer program, also performs the following steps:

[0173] The tomographic inversion objective function after regularization constraints is solved using the preconditional conjugate gradient method to obtain the velocity update around the fracture.

[0174] Example 5

[0175] In this embodiment, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, it performs the following steps:

[0176] Obtain seismic profiles;

[0177] Extract coherent properties from the seismic profile;

[0178] The resampled coherent properties are obtained by resampling the coherent properties to obtain the resampled coherent properties corresponding to the tomographic grid;

[0179] The ray density is obtained through ray tracing;

[0180] The crack control bundle tomography operator is calculated based on the resampling coherence properties and the ray density.

[0181] The crack control bundle tomography operator is introduced into the tomography inversion objective function for regularization constraint to obtain the tomography inversion objective function after regularization constraint.

[0182] In one embodiment, when the computer program is executed by a processor, it also performs the following steps:

[0183] Based on the seismic profile, the imaging results of each wheel in the migration velocity modeling are obtained;

[0184] The coherence properties are extracted based on the previous imaging results in the offset velocity modeling.

[0185] In one embodiment, when the computer program is executed by a processor, it further implements the step of resampling the coherent properties to obtain the resampled coherent properties corresponding to the tomographic grid, wherein the resampling calculation formula is:

[0186]

[0187] in This indicates that the position after resampling to the tomographic grid size is a point. Relevance properties This indicates that the location point was calculated directly using the results of the previous imaging round. Relevance properties This represents the resampling operator.

[0188] In one embodiment, when the computer program is executed by a processor, it also performs the following steps:

[0189] Based on the resampling coherence properties and the ray density The fracture control bundle tomography operator is calculated using the following formula:

[0190]

[0191] in, The grid points are The location of the crack control bundle tomography operator; Indicates that it is located at grid point The calculation time window is The resampling coherence property, Indicates that it is located at grid point Time window The ray density value, , and These are the half-time window lengths in three directions.

[0192] In one embodiment, when the computer program is executed by the processor, it further implements the introduction of the crack control bundle tomography operator into the tomography inversion objective function for regularization constraints. The reference method is to use a preconditional regularization operator, and the calculation formula is:

[0193]

[0194] in, For the tomographic inversion linearization operator, For the crack control bundle tomography operator, The solution is given by the preconditions.

[0195] In one embodiment, when the computer program is executed by a processor, it also performs the following steps:

[0196] Solving the tomographic inversion objective function after regularization constraints yields the updated velocity around the fracture.

[0197] In one embodiment, when the computer program is executed by a processor, it also performs the following steps:

[0198] The tomographic inversion objective function after regularization constraints is solved using the preconditional conjugate gradient method to obtain the velocity update around the fracture.

[0199] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium. When executed, the computer program can include the processes of the embodiments of the above methods. Any references to memory, storage, databases, or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile memory may include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory may include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), dual data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), RAMbus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0200] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0201] The embodiments described above are merely illustrative of several implementation methods of this application, and while the descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the invention patent. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this application, and these all fall within the protection scope of this application. Therefore, the protection scope of this patent application should be determined by the appended claims.

Claims

1. A fracture control bundle tomography regularization method, characterized in that, include: Obtain seismic profiles; Extract coherent properties from the seismic profile; The resampled coherent properties are obtained by resampling the coherent properties to obtain the resampled coherent properties corresponding to the tomographic grid; The ray density is obtained through ray tracing; The crack control bundle tomography operator is calculated based on the resampling coherence properties and the ray density. The crack control bundle tomography operator is introduced into the tomography inversion objective function for regularization constraint to obtain the tomography inversion objective function after regularization constraint; The step of extracting coherent attributes from the seismic profile includes: Based on the seismic profile, the imaging results of each wheel in the migration velocity modeling are obtained; The coherence attribute is extracted based on the previous imaging result in the offset velocity modeling. In the step of resampling the coherent attributes to obtain the resampled coherent attributes corresponding to the tomographic grid, the calculation formula for resampling is: in This indicates that the position after resampling to the tomographic grid size is a point. Relevance properties This indicates that the location point was calculated directly using the results of the previous imaging round. Relevance properties Represents the resampling operator; The step of calculating the crack control beam tomography operator based on the resampling coherence properties and the ray density includes: Based on the resampling coherence properties and the ray density The fracture control bundle tomography operator is calculated using the following formula: in, The grid points are The location of the crack control bundle tomography operator; Indicates that it is located at grid point The calculation time window is The resampling coherence property, Indicates that it is located at grid point Time window The ray density value, , and These are the half-time window lengths in three directions.

2. The method according to claim 1, characterized in that, In the step of introducing the fracture control bundle tomography operator into the tomography inversion objective function for regularization constraints, the method of reference is to use a preconditional regularization operator, and the calculation formula is: in, For the tomographic inversion linearization operator, For the crack control bundle tomography operator, The solution is given by the preconditions.

3. The method according to claim 1, characterized in that, The step of obtaining the tomographic inversion objective function after obtaining the regularization constraint further includes: Solving the tomographic inversion objective function after regularization constraints yields the updated velocity around the fracture.

4. The method according to claim 3, characterized in that, The step of solving the tomographic inversion objective function after regularization constraints to obtain the updated velocity around the fracture includes: The tomographic inversion objective function after regularization constraints is solved using the preconditional conjugate gradient method to obtain the velocity update around the fracture.

5. A fracture control bundle tomography regularization device based on the fracture control bundle tomography regularization method according to claim 1, characterized in that, include: The seismic profile acquisition module is used to acquire seismic profiles. The coherence attribute extraction module is used to extract coherence attributes based on the seismic profile; The resampling module is used to resample the coherent properties to obtain the resampled coherent properties corresponding to the tomographic grid; A ray density acquisition module is used to obtain ray density through ray tracing. The control beam tomography operator calculation module is used to calculate the crack control beam tomography operator based on the resampling coherence properties and the ray density; The regularization module is used to introduce the fracture control bundle tomography operator into the tomography inversion objective function for regularization constraints, so as to obtain the tomography inversion objective function after regularization constraints.

6. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 4.

7. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 4.

Images