Well control chromatography static correction method and system for loess highland and medium

Abstract

The invention provides a well control chromatography static correction method and system for a loess highland and a medium, and belongs to the field of oil and gas geophysical exploration. The method comprises the following steps: step 100, inputting seismic data; step 200, performing tomographic static correction processing on the seismic data to obtain a tomographic static correction value; step 300, acquiring a well-to-seismic error value of the whole area; step 400, applying the well-to-seismic error value of the whole area to the chromatography static correction value to obtain a new chromatography static correction value; step 500, performing inversion according to the new tomography static correction value to obtain a new near-surface velocity model; and step 600, adopting the new near-surface velocity model to calculate and obtain a new static correction value. According to the method, the seismic condition, the current collection situation and the chromatography inversion characteristics of the loess highland are comprehensively considered, the stratified data of the well are utilized to calibrate, constrain and correct the static correction value error caused after chromatography static correction, finally, the well seismic error is reduced, and reasonable seismic structure imaging is ensured.

Description

Technical Field

[0001] This invention belongs to the field of oil and gas geophysical exploration, specifically relating to a well-controlled tomography static correction method, system, and medium in the Loess Plateau region. Background Technology

[0002] In the processing of 3D seismic data in the Loess Plateau region, we often face the problem of static correction caused by the lateral changes in near-surface structure and the complex and variable surface excitation and reception conditions. This can lead to inaccurate static correction, causing well-seismic errors and reducing the reliability of the data.

[0003] Currently, various tomographic static correction methods (such as refractive tomographic static correction and grid tomographic static correction) are generally used to address static correction problems in the Loess Plateau region, but each method has its own limitations. In general areas, due to the constraints of low-gauge data (usually micro-well data), micro-well-constrained tomographic static correction or travel-time tomographic static correction methods can effectively correct near-surface influences and achieve good results. However, in the Loess Plateau region, the thick loess layer causes distortion in near-surface low-gauge data. The constrained inversion based on this low-gauge data cannot accurately reflect near-surface velocity changes, leading to inaccurate static correction values ​​in the tomographic static correction. Ultimately, this results in excessive well-seismic errors, making it difficult to achieve the desired correction effect. Summary of the Invention

[0004] The purpose of this invention is to solve the problems existing in the prior art and provide a well-controlled tomography static correction method, system and medium for the Loess Plateau region. It comprehensively considers the seismic conditions, current acquisition status and tomographic inversion characteristics of the Loess Plateau region, uses the layered data of the well to calibrate constraints, corrects the static correction error brought about by the tomographic static correction, and ultimately reduces well-seismic error to ensure reasonable seismic tectonic imaging.

[0005] This invention is achieved through the following technical solution:

[0006] The first aspect of the present invention provides a well-controlled tomography static correction method for the Loess Plateau region, comprising the following steps:

[0007] Step 100: Input earthquake data;

[0008] Step 200: Perform tomographic static correction on the seismic data to obtain the tomographic static correction amount;

[0009] Step 300: Obtain the well seismic error values ​​for the entire area;

[0010] Step 400: Apply the well seismic error values ​​of the entire area to the tomographic static correction value to obtain a new tomographic static correction value;

[0011] Step 500: Obtain a new near-surface velocity model by inverting the new tomographic static correction value;

[0012] Step 600: The new static correction is calculated using the new near-surface velocity model.

[0013] A further improvement of the present invention is that:

[0014] Step 300: Obtain the well seismic error value for the entire area. Specific operations include:

[0015] Step 310: Apply the tomographic static correction to the seismic data to generate stacked data, and load the generated stacked data into the seismic data interpretation system to pick the target horizon.

[0016] Step 320: Obtain the well vibration error value;

[0017] Step 330: Calculate the well-seismic error value for the entire area.

[0018] A further improvement of the present invention is that:

[0019] Step 320 involves obtaining the well vibration error value, specifically including the following operations:

[0020] Subtract the time of the target layer obtained after the well logging curve is calibrated by the synthetic record from the time of the target layer picked in step 310 to obtain the error value between the well's layer data and the picked layer, i.e., the well vibration error value.

[0021] A further improvement of the present invention is that:

[0022] Step 330 calculates the well-seismic error value for the entire area. Specific operations include:

[0023] Steps 310 and 320 are performed on all wells within the work area to obtain the corresponding error values ​​for all wells. The error values ​​of all wells are then linearly interpolated to generate the overall well seismic error value for the entire area.

[0024] A further improvement of the present invention is that:

[0025] In step 500, a new near-surface velocity model is obtained by inversion based on the new tomographic static correction. Specific operations include:

[0026] The new tomographic static corrections are inverted to generate a shallow near-surface velocity model.

[0027] The shallow near-surface velocity model generated by the inversion is subjected to constrained tomographic inversion to generate a new near-surface velocity model.

[0028] A further improvement of the present invention is that:

[0029] The well-controlled tomography static correction method for the Loess Plateau region also includes:

[0030] Step 700: Apply the new static correction value calculated in step 600 to the seismic data and repeat the operations of steps 200-600 until the error value of the new tomographic static correction value obtained in step 600 is less than the set value.

[0031] A second aspect of the present invention provides a well-controlled tomography static correction system for the Loess Plateau region, comprising:

[0032] The input module is used to input seismic data;

[0033] The tomographic static correction module is used to perform tomographic static correction processing on seismic data to obtain the tomographic static correction amount;

[0034] The well vibration error acquisition module is used to acquire well vibration error values ​​for the entire area.

[0035] The tomographic static correction update module is used to apply the well-seismic error values ​​of the entire area to the tomographic static correction values ​​to obtain new tomographic static correction values.

[0036] The inversion module is used to invert the new near-surface velocity model based on the new tomographic static correction.

[0037] The static correction calculation module uses a new near-surface velocity model to calculate new static corrections.

[0038] A further improvement of the present invention is that:

[0039] The well vibration error acquisition module includes:

[0040] The target horizon picking submodule is used to apply tomographic static corrections to seismic data, generate stacked data, and load the generated stacked data into the seismic data interpretation system for target horizon picking.

[0041] The first acquisition submodule is used to acquire well vibration error values;

[0042] The second acquisition submodule is used to calculate the well seismic error value for the entire area.

[0043] A further improvement of the present invention is that:

[0044] The inversion module includes:

[0045] The first inversion submodule is used to invert the new tomographic static corrections to generate a shallow near-surface velocity model.

[0046] The second inversion submodule is used to perform constrained tomographic inversion on the shallow near-surface velocity model generated by the inversion, and generate a new near-surface velocity model.

[0047] A third aspect of the present invention provides a computer-readable storage medium storing at least one computer-executable program, which, when executed by the computer, causes the computer to perform the steps in the well-controlled tomography static correction method for the Loess Plateau region.

[0048] Compared with the prior art, the beneficial effects of the present invention are:

[0049] This invention comprehensively considers the seismic conditions, current acquisition status, and tomographic inversion characteristics of the Loess Plateau region. It utilizes layered well data to calibrate constraints, corrects the static correction error caused by tomographic static correction, and ultimately reduces well-seismic error, ensuring reasonable seismic tectonic imaging. Attached Figure Description

[0050] Figure 1 This is a flowchart of a well-controlled tomography static correction method for the Loess Plateau region according to the present invention;

[0051] Figure 2 This is a diagram showing the effect of static elevation correction.

[0052] Figure 3 This is a graph showing the effect of static calibration in chromatography;

[0053] Figure 4 This is a diagram showing the static correction effect of well-controlled tomography according to the present invention. Detailed Implementation

[0054] The present invention will now be described in further detail with reference to the accompanying drawings:

[0055] In the processing of 3D seismic data in the Loess Plateau region, the lateral variations of near-surface structures and the complex and variable surface excitation and reception conditions often lead to static correction issues. This can result in inaccurate static corrections, causing well-to-seismic errors and reducing the reliability of the data. This invention addresses these problems by comprehensively considering the seismic conditions, current acquisition status, and tomographic inversion characteristics of the Loess Plateau region. It utilizes layered well data to calibrate the tomographically corrected seismic profile, correcting the static correction errors introduced by the tomographic static correction, ultimately reducing well-to-seismic errors and ensuring accurate seismic tectonic imaging.

[0056]

Example 1

[0057] This invention provides a well-controlled tomography static correction method for the Loess Plateau region, such as... Figure 1 As shown, the specific steps include:

[0058] Step 100: Input earthquake data;

[0059] Step 200: Perform tomographic static correction on the seismic data to obtain the tomographic static correction amount;

[0060] Step 300: Obtain the well seismic error for the entire area;

[0061] Step 400: Apply the well seismic error of the entire area to the tomographic static correction value to obtain a new tomographic static correction value;

[0062] Step 500: Obtain a new near-surface velocity model by inverting the new tomographic static correction value;

[0063] Step 600: The new static correction is calculated using the new near-surface velocity model.

[0064] The static calibration process in step 200 uses existing techniques and will not be described in detail here.

[0065] This invention comprehensively considers the seismic conditions, current acquisition status, and tomographic inversion characteristics of the Loess Plateau region. It uses layered well data to calibrate the seismic profile after static tomographic correction, corrects the static correction error caused by static tomographic correction, and ultimately reduces well-seismic error, ensuring reasonable seismic tectonic imaging.

[0066]

Example 2

[0067] Step 300: Obtain the well seismic error for the entire area. Specific operations include:

[0068] Step 310: Apply the tomographic static correction to the seismic data to generate stacked data, and load the generated stacked data into the seismic data interpretation system to pick the target horizon.

[0069] It should be understood that at the seismic trace head, the tomographic static correction value is recorded in a trace head word. That is, the value of the tomographic static correction value is assigned to a trace head word. When applied, this trace head word will add the value to this seismic trace. After superposition processing, multiple seismic traces at the same location are superimposed into one, thereby generating the entire superposition data.

[0070] It should be understood that the seismic data interpretation system is used for geological interpretation, including stratigraphic calibration, target layer tracking, and three-dimensional depiction of structural maps. It is a mature and readily available product, and will not be elaborated upon here.

[0071] Target horizon selection specifically involves: after determining the target horizon on the seismic profile, it can be manually selected along the main survey line using the selection tools provided by the seismic data interpretation software, or automatically selected using the tools provided by the seismic data interpretation software after defining several key seed points, and finally manually modified. This is an existing technical method and will not be elaborated further here.

[0072] Step 320: Obtain the well vibration error. Specific operations include:

[0073] Subtract the time of the target layer obtained after the well logging curve is calibrated by the synthetic record from the time of the target layer picked in step 310 to obtain the error between the well's layer data and the picked layer, i.e., well vibration error.

[0074] Step 330: Calculate the well-seismic error for the entire area. Specific operations include:

[0075] Perform steps 310 and 320 on all wells within the work area to obtain the corresponding error values ​​for all wells. Then, perform linear interpolation on the error values ​​of all wells to generate the error value for the entire area.

[0076]

Example 3

[0077] Step 400: Apply the well seismic error of the entire area to the tomographic static correction value to obtain a new tomographic static correction value.

[0078] It should be understood that the well-seismic error value for the entire area refers to the corresponding error value for both the shot point and the receiver point in the entire work area. This error value is added one by one to the original tomographic static correction value decomposed to the shot point and receiver point to obtain the updated tomographic static correction value decomposed to the shot point and receiver point.

[0079]

Example 4

[0080] In step 500, a new near-surface velocity model is obtained by inversion based on the new tomographic static correction. Specific operations include:

[0081] The new tomographic static corrections are inverted to generate a shallow near-surface velocity model.

[0082] The shallow near-surface velocity model generated by the inversion is subjected to constrained tomographic inversion to generate a new near-surface velocity model.

[0083] The inversion method in this embodiment is an existing technology and will not be described in detail here.

[0084]

Example 5

[0085] Step 600: Calculate the new static correction amount using the new near-surface velocity model.

[0086] It should be understood that the near-surface velocity model can directly calculate the static correction, and the basic calculation formula is as follows:

[0087] △T=h / V

[0088] Where △T is the static correction amount, h is the model thickness, and V is the model velocity.

[0089]

Example 6

[0090] An embodiment of the present invention provides a well-controlled tomography static correction method for the Loess Plateau region, which further includes:

[0091] Step 700: Apply the new static correction value calculated in step 600 to the seismic data and repeat the operations of steps 200-600 until the error value of the new tomographic static correction value obtained in step 600 is less than the set value.

[0092] Figure 2 This is a diagram showing the effect of static elevation correction. Figure 3 This is a graph showing the effect of static calibration in chromatography; Figure 4 This is a static correction effect diagram of well-controlled tomography according to the present invention. A comparison shows that, in terms of imaging effect, Figure 4 The method of this invention provides more accurate focusing on the phase axis of the target layer. In terms of results, it can better reflect the true formation state and is more consistent with well seismic data in actual drilling data comparison with smaller errors.

[0093]

Example 7

[0094] This invention provides a well-controlled tomography static correction system for the Loess Plateau region, comprising:

[0095] The input module is used to input seismic data;

[0096] The tomographic static correction module is used to perform tomographic static correction processing on seismic data to obtain the tomographic static correction amount;

[0097] The well vibration error acquisition module is used to acquire the well vibration error of the entire area;

[0098] The tomographic static correction update module is used to apply the well seismic error of the entire area to the tomographic static correction to obtain a new tomographic static correction.

[0099] The inversion module is used to invert the new near-surface velocity model based on the new tomographic static correction.

[0100] The static correction calculation module uses a new near-surface velocity model to calculate new static corrections.

[0101] The well vibration error acquisition module includes:

[0102] The target horizon picking submodule is used to apply tomographic static corrections to seismic data, generate stacked data, and load the generated stacked data into the seismic data interpretation system for target horizon picking.

[0103] The first acquisition submodule is used to acquire well vibration error values;

[0104] The second acquisition submodule is used to calculate the well seismic error value for the entire area.

[0105] The inversion module includes:

[0106] The first inversion submodule is used to invert the new tomographic static corrections to generate a shallow near-surface velocity model.

[0107] The second inversion submodule is used to perform constrained tomographic inversion on the shallow near-surface velocity model generated by the inversion, and generate a new near-surface velocity model.

[0108]

Example 8

[0109] This invention provides a computer-readable storage medium storing at least one computer-executable program, which, when executed by the computer, causes the computer to perform the steps in the well-controlled tomography static correction method for the Loess Plateau region.

[0110] 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, and when executed, it 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 can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various 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.

[0111] The above technical solution is only one embodiment of the present invention. For those skilled in the art, based on the principles disclosed in the present invention, it is easy to make various types of improvements or modifications, and not limited to the technical solutions described in the specific embodiments of the present invention. Therefore, the foregoing description is only a preferred option and is not restrictive.

Claims

1. A well-controlled tomography static correction method for the Loess Plateau region, characterized in that, Includes the following steps: Step 100: Input earthquake data; Step 200: Perform tomographic static correction on the seismic data to obtain the tomographic static correction amount; Step 300: Obtain the well seismic error values ​​for the entire area; Step 400: Apply the well seismic error values ​​of the entire area to the tomographic static correction value to obtain a new tomographic static correction value; Step 500: Obtain a new near-surface velocity model by inverting the new tomographic static correction value; Step 600: The new static correction is calculated using the new near-surface velocity model.

2. The method according to claim 1, characterized in that, Step 300: Obtain the well seismic error value for the entire area. Specific operations include: Step 310: Apply the tomographic static correction to the seismic data to generate stacked data, and load the generated stacked data into the seismic data interpretation system to pick the target horizon. Step 320: Obtain the well vibration error value; Step 330: Calculate the well-seismic error value for the entire area.

3. The method according to claim 2, characterized in that, Step 320 involves obtaining the well vibration error value, specifically including the following operations: Subtract the time of the target layer obtained after the well logging curve is calibrated by the synthetic record from the time of the target layer picked in step 310 to obtain the error value between the well's layer data and the picked layer, i.e., the well vibration error value.

4. The method according to claim 3, characterized in that, Step 330 calculates the well-seismic error value for the entire area. Specific operations include: Steps 310 and 320 are performed on all wells within the work area to obtain the corresponding error values ​​for all wells. The error values ​​of all wells are then linearly interpolated to generate the overall well seismic error value for the entire area.

5. The method according to claim 1, characterized in that, In step 500, a new near-surface velocity model is obtained by inversion based on the new tomographic static correction. Specific operations include: The new tomographic static corrections are inverted to generate a shallow near-surface velocity model. The shallow near-surface velocity model generated by the inversion is subjected to constrained tomographic inversion to generate a new near-surface velocity model.

6. The method according to claim 1, characterized in that, The well-controlled tomography static correction method for the Loess Plateau region also includes: Step 700: Apply the new static correction value calculated in step 600 to the seismic data and repeat the operations of steps 200-600 until the error value of the new tomographic static correction value obtained in step 600 is less than the set value.

7. A well-controlled tomography static correction system for the Loess Plateau region, characterized in that, include: The input module is used to input seismic data; The tomographic static correction module is used to perform tomographic static correction processing on seismic data to obtain the tomographic static correction amount; The well vibration error acquisition module is used to acquire well vibration error values ​​for the entire area. The tomographic static correction update module is used to apply the well-seismic error values ​​of the entire area to the tomographic static correction values ​​to obtain new tomographic static correction values. The inversion module is used to invert the new near-surface velocity model based on the new tomographic static correction. The static correction calculation module uses a new near-surface velocity model to calculate new static corrections.

8. The system according to claim 7, characterized in that, The well vibration error acquisition module includes: The target horizon picking submodule is used to apply tomographic static corrections to seismic data, generate stacked data, and load the generated stacked data into the seismic data interpretation system for target horizon picking. The first acquisition submodule is used to acquire well vibration error values; The second acquisition submodule is used to calculate the well seismic error value for the entire area.

9. The system according to claim 7, characterized in that, The inversion module includes: The first inversion submodule is used to invert the new tomographic static corrections to generate a shallow near-surface velocity model. The second inversion submodule is used to perform constrained tomographic inversion on the shallow near-surface velocity model generated by the inversion, and generate a new near-surface velocity model.

10. A computer-readable storage medium storing at least one computer-executable program, which, when executed by the computer, causes the computer to perform the steps of the well-controlled tomography static correction method for the Loess Plateau region as described in any one of claims 1-6.

Images