Seismic data smoothing method and device based on slope information
Through a method based on slope information, the slope information in the seismic data is obtained and the Gaussian weighted smoothing process is performed, which solves the problem of excessive smoothing damage to the data information in the prior art, and achieves efficient data smoothing and signal detail protection.
Patent Information
- Application Number
- CN202311126137.7
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2023-09-01
- Publication Date
- 2025-03-04
- Estimated Expiration
- 2043-09-01
AI Technical Summary
When seismic data smoothing is performed, excessive smoothing will destroy the original data information and the structure pattern is not considered, resulting in blurred signal details.
A seismic data smoothing method based on slope information is proposed. By obtaining the slope information of each point in the three-dimensional seismic data, the inverse tangent function value is calculated to obtain the inclination angle, the smooth input point is constructed, and the Gaussian coefficient is used for Gaussian weighted smoothing processing.
It achieves a better smoothing effect, while protecting the original seismic data information, not destroying the tectonic form, and improving the signal-to-noise ratio.
Smart Images

Figure CN119556346A_ABST
Abstract
Description
Technical Field
[0001] The present application belongs to the field of geophysical exploration technology, and specifically relates to a seismic data smoothing method and device based on slope information. Background Art
[0002] For the actual seismic data collected, there are often noise, abnormal signals and other influences, which lead to the decline of data quality and make it difficult to process and obtain better results. At this time, data smoothing technology is usually used to suppress or eliminate the influence of noise, abnormal signals and other factors to make the data more regular. However, excessive smoothing will destroy a lot of information in the original data and blur the details in the signal. At present, conventional smoothing methods have not considered the structural morphology, but directly use multi-point smoothing operators to perform filtering along three dimensions. Although it can achieve a good smoothing effect, it will destroy the original data information. Summary of the invention
[0003] Based on the above technical problems, the present application proposes a seismic data smoothing method and device based on slope information.
[0004] In a first aspect, the present application proposes a seismic data smoothing method based on slope information, comprising:
[0005] Step S1: Obtain slope information for each point in the three-dimensional seismic data;
[0006] Step S2: Calculate the inverse tangent function value of the slope information to obtain the inclination angle of the corresponding point;
[0007] Step S3: constructing smoothing input points required for smoothing according to the inclination angle;
[0008] Step S4: using Gaussian coefficients to perform Gaussian weighted smoothing on the smoothed input points to obtain smoothed results of the smoothed input points. The smoothed results of the smoothed input points of each point constitute the smoothed results of the seismic data.
[0009] The x-axis direction of the three-dimensional seismic data is data in the direction of the seismic main survey line, the y-axis direction is data in the direction of the seismic contact survey line, and the z-axis direction is data in the direction of the seismic depth;
[0010] The step of obtaining slope information for each point in the three-dimensional seismic data includes:
[0011] For each point in the three-dimensional seismic data, obtain slope information on the plane formed by the x-axis and the z-axis;
[0012] For each point in the 3D seismic data, the slope information on the plane formed by the y-axis and the z-axis is obtained.
[0013] The calculating the inverse tangent function value of the slope information to obtain the inclination angle of the corresponding point includes:
[0014] Using the slope information on the plane formed by the x-axis and the z-axis, calculate the inclination angle of the corresponding point on the plane formed by the x-axis and the z-axis;
[0015] Using the slope information on the plane formed by the y-axis and the z-axis, the inclination angle of the corresponding point on the plane formed by the y-axis and the z-axis is calculated. The calculation formula is as follows:
[0016]
[0017] Among them, S xz (i, j, k) is the slope information on the plane formed by the x-axis and the z-axis, S yz (i, j, k) is the slope information on the plane formed by the y-axis and the z-axis, θ xz (i,j,k) is the inclination angle of point (i,j,k) on the plane formed by the x-axis and the z-axis, θ yz (i,j,k) is the inclination angle of point (i,j,k) on the plane formed by the y-axis and the z-axis, and (i,j,k) is the coordinates of a point in the underground stratum.
[0018] The step of constructing a smoothing input point required for smoothing according to the inclination angle comprises:
[0019] Obtaining a smooth input point in the x-axis direction according to the inclination angle of the corresponding point on the plane formed by the x-axis and the z-axis;
[0020] According to the inclination angle of the corresponding point on the plane formed by the y-axis and the z-axis, a smooth input point in the y-axis direction is obtained.
[0021] The smooth input point in the x-axis direction is expressed as follows:
[0022]
[0023] Where N is the number of smoothing input points, N-1 is an integer multiple of 4, (i, j, k) is the coordinate of a point in the underground stratum, as the first smoothing input point, (i-1, j, k-ceil{1×tan[θ xz (i,j,k)]}) is the second smooth input point in the x-axis direction, is the first Smooth input point, (i+1,j,k+ceil{1×tan[θ xz (i,j,k)]}) is the first smooth input points, In the x-axis direction smooth input points, tan[θxz (i, j, k)] is the tangent function value of the inclination angle on the plane formed by the x-axis and the z-axis, and ceil{} represents rounding the content in {}.
[0024] The smooth input point in the y-axis direction is calculated as follows:
[0025]
[0026] Where N is the number of smoothing input points, N-1 is an integer multiple of 4, (i, j, k) is the coordinate of a point in the underground stratum, which is used as the first smoothing input point. is the second smooth input point in the y-axis direction, is the first Smooth input point, (i,j+1,k+ceil{1×tan[θ yz (i,j,k)]}) is the first smooth input points, is the first smooth input points, tan[θ yz (i, j, k)] is the tangent function value of the inclination angle on the plane formed by the y-axis and the z-axis, and ceil{} means rounding the content in {}.
[0027] The Gaussian coefficient is:
[0028]
[0029] Among them, G(i,j,k) is the Gaussian coefficient of point (i,j,k), σ is the standard deviation of the smoothed input point, l is the window width,
[0030] In a second aspect, the present application proposes a seismic data smoothing device based on slope information, comprising:
[0031] An information acquisition module, used for acquiring slope information for each point in the three-dimensional seismic data;
[0032] The inclination angle calculation module is used to calculate the inverse tangent function value of the slope information to obtain the inclination angle of the corresponding point;
[0033] A smoothing point construction module, used for constructing smoothing input points required for smoothing according to the inclination angle;
[0034] The Gaussian weighting module is used to perform Gaussian weighted smoothing processing on the smoothing input points using Gaussian coefficients to obtain smoothing results of the smoothing input points. The smoothing results of the smoothing input points of each point constitute the smoothing results of the seismic data.
[0035] In a third aspect, the present application proposes an electronic device, comprising: a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the seismic data smoothing method based on slope information is executed.
[0036] In a fourth aspect, the present application proposes a computer-readable storage medium storing executable instructions, which, when executed, enable a processor to execute the seismic data smoothing method based on slope information.
[0037] Beneficial effects:
[0038] The present application proposes a seismic data smoothing method and device based on slope information, which realizes a seismic data smoothing method based on slope information, which can achieve a better smoothing effect without destroying the original seismic data information. BRIEF DESCRIPTION OF THE DRAWINGS
[0039] Figure 1 A flow chart of a seismic data smoothing method based on slope information according to an embodiment of the present application;
[0040] Figure 2 A schematic diagram of the formation slope in the x direction of an embodiment of the present application;
[0041] Figure 3 A schematic diagram of the formation slope in the y direction of an embodiment of the present application;
[0042] Figure 4 A schematic diagram of data to be smoothed obtained by constructing an embodiment of the present application;
[0043] Figure 5 This is a schematic diagram of a data body after smoothing based on a slope according to an embodiment of the present application;
[0044] Figure 6 This is a principle block diagram of a seismic data smoothing device based on slope information according to an embodiment of the present application. DETAILED DESCRIPTION
[0045] The present disclosure is further described below in conjunction with the embodiments shown in the accompanying drawings.
[0046] In the prior art, data smoothing technology is usually used to suppress or eliminate the influence of noise, abnormal signals and other factors to make the data more regular. However, excessive smoothing will destroy a lot of information in the original data and blur the details in the signal. The prior art has not considered the structural morphology, but directly uses a multi-point smoothing operator to perform filtering along three dimensions. Although it can achieve a good smoothing effect, it will destroy the original data information.
[0047] The present application proposes a seismic data smoothing method based on slope information. First, the slope information of each point in the three-dimensional seismic data is read in, and then the corresponding smoothing input points are calculated and extracted along the slope information. The Gaussian coefficient is used to perform multi-point smoothing on the smoothing input points under slope constraints, thereby realizing a data smoothing method based on slope information. This method can achieve a better smoothing effect without destroying the original seismic data information.
[0048] Embodiment 1
[0049] This embodiment proposes a seismic data smoothing method based on slope information, such as Figure 1 As shown, including:
[0050] Step S1: Obtain slope information for each point in the three-dimensional seismic data;
[0051] In this embodiment, the x-axis direction of the three-dimensional seismic data is the data in the direction of the seismic main survey line, the y-axis direction is the data in the direction of the seismic contact survey line, and the z-axis direction is the data in the direction of the seismic depth;
[0052] The step of obtaining slope information for each point in the three-dimensional seismic data includes:
[0053] For each point in the three-dimensional seismic data, obtain slope information on the plane formed by the x-axis and the z-axis;
[0054] For each point in the 3D seismic data, the slope information on the plane formed by the y-axis and the z-axis is obtained.
[0055] Step S2: Calculate the inverse tangent function value of the slope information to obtain the inclination angle of the corresponding point, including:
[0056] Using the slope information on the plane formed by the x-axis and the z-axis, calculate the inclination angle of the corresponding point on the plane formed by the x-axis and the z-axis;
[0057] Using the slope information on the plane formed by the y-axis and the z-axis, the inclination angle of the corresponding point on the plane formed by the y-axis and the z-axis is calculated. The calculation formula is as follows:
[0058]
[0059] Among them, S xz (i, j, k) is the slope information on the plane formed by the x-axis and the z-axis, S yz (i, j, k) is the slope information on the plane formed by the y-axis and the z-axis, θ xz (i,j,k) is the inclination angle of point (i,j,k) on the plane formed by the x-axis and the z-axis, θ yz(i,j,k) is the inclination angle of point (i,j,k) on the plane formed by the y-axis and the z-axis, and (i,j,k) is the coordinates of a point in the underground stratum.
[0060] Step S3: constructing smoothing input points required for smoothing according to the inclination angle, including:
[0061] Obtaining a smooth input point in the x-axis direction according to the inclination angle of the corresponding point on the plane formed by the x-axis and the z-axis;
[0062] According to the inclination angle of the corresponding point on the plane formed by the y-axis and the z-axis, a smooth input point in the y-axis direction is obtained.
[0063] The smooth input point in the x-axis direction is expressed as follows:
[0064]
[0065] Where N is the number of smoothing input points, N-1 is an integer multiple of 4, (i, j, k) is the coordinate of a point in the underground stratum, as the first smoothing input point, (i-1, j, k-ceil{1×tan[θ xz (i,j,k)]}) is the second smooth input point in the x-axis direction, In the x-axis direction Smooth input point, (i+1,j,k+ceil{1×tan[θ xz (i,j,k)]}) is the first smooth input points, In the x-axis direction smooth input points, tan[θ xz (i, j, k)] is the tangent function value of the inclination angle on the plane formed by the x-axis and the z-axis, and ceil{} means rounding the content in {}.
[0066] The smooth input point in the y-axis direction is calculated as follows:
[0067]
[0068] Where N is the number of smoothing input points, N-1 is an integer multiple of 4, (i, j, k) is the coordinate of a point in the underground stratum, which is used as the first smoothing input point. is the second smooth input point in the y-axis direction, is the first Smooth input point, (i,j+1,k+ceil{1×tan[θ yz (i,j,k)]}) is the first smooth input points, is the first smooth input points, tan[θ yz (i, j, k)] is the tangent function value of the inclination angle on the plane formed by the y-axis and the z-axis, and ceil{} means rounding the content in {}.
[0069] Step S4: using Gaussian coefficients to perform Gaussian weighted smoothing on the smoothed input points to obtain smoothed results of the smoothed input points. The smoothed results of the smoothed input points of each point constitute the smoothed results of the seismic data.
[0070] Assuming that the seismic data of point (i, j, k) is W(i, j, k), the Gaussian weighted multi-point smoothing at point (i, j, k) can be expressed as:
[0071]
[0072] in, Smoothed result of sliding input points.
[0073] The Gaussian coefficient is:
[0074]
[0075] Among them, G(i,j,k) is the Gaussian coefficient of point (i,j,k), σ is the standard deviation of the smoothed input point, l is the window width, Traverse all points in the three-dimensional seismic data to obtain the smoothing results of the smoothing input points of each point, and finally form the smoothing results of the seismic data.
[0076] The present embodiment proposes a seismic data smoothing method based on slope information. First, the slope information of each point is obtained; the inverse tangent function value of the slope information is calculated to obtain the inclination angle of the corresponding point; according to the inclination angle, the smoothing input point required for smoothing is constructed; the smoothing input point is Gaussian weighted smoothing is performed on the smoothing input point using a Gaussian coefficient to obtain the smoothing result of the smoothing input point, and the smoothing result of the smoothing input point of each point constitutes the seismic data smoothing result. The method of the present embodiment can achieve a better smoothing effect without destroying the original seismic data information.
[0077] Embodiment 2
[0078] This embodiment is a specific example of the first embodiment. The real three-dimensional seismic data is input and nine-point smoothing is used as an example to describe in detail the seismic data smoothing method based on slope information. Figure 1 As shown, including:
[0079] Step S1: Obtain slope information for each point in the three-dimensional seismic data;
[0080] In this embodiment, the x-axis direction of the three-dimensional seismic data is the data in the direction of the seismic main survey line, the y-axis direction is the data in the direction of the seismic contact survey line, and the z-axis direction is the data in the direction of the seismic depth; Figure 2 As shown in FIG. 1 , the x-direction formation slope input in this embodiment is used to calculate the specific position of the smoothing input point required for x-direction smoothing. The x-direction formation slope is the slope information on the plane formed by the x-axis and the z-axis. Figure 3 As shown, it is the formation slope in the y direction input in this embodiment, which is used to calculate the specific position of the smooth input point required in the y direction. The formation slope in the y direction is the slope information on the plane formed by the y axis and the z axis. Obtaining slope information for each point in the three-dimensional seismic data includes:
[0081] For each point in the three-dimensional seismic data, obtain slope information on the plane formed by the x-axis and the z-axis;
[0082] For each point in the 3D seismic data, the slope information on the plane formed by the y-axis and the z-axis is obtained.
[0083] Step S2: Calculate the inverse tangent function value of the slope information to obtain the inclination angle of the corresponding point, including:
[0084] Using the slope information on the plane formed by the x-axis and the z-axis, calculate the inclination angle of the corresponding point on the plane formed by the x-axis and the z-axis;
[0085] Using the slope information on the plane formed by the y-axis and the z-axis, the inclination angle of the corresponding point on the plane formed by the y-axis and the z-axis is calculated. The calculation formula is as follows:
[0086]
[0087] Among them, S xz (i, j, k) is the slope information on the plane formed by the x-axis and the z-axis, S yz (i, j, k) is the slope information on the plane formed by the y-axis and the z-axis, θ xz (i,j,k) is the inclination angle of point (i,j,k) on the plane formed by the x-axis and the z-axis, θ yz (i,j,k) is the inclination angle of point (i,j,k) on the plane formed by the y-axis and the z-axis, and (i,j,k) is the coordinates of a point in the underground stratum.
[0088] Step S3: constructing smoothing input points required for smoothing according to the inclination angle, including:
[0089] Obtaining a smooth input point in the x-axis direction according to the inclination angle of the corresponding point on the plane formed by the x-axis and the z-axis;
[0090] According to the inclination angle of the corresponding point on the plane formed by the y-axis and the z-axis, a smooth input point in the y-axis direction is obtained.
[0091] This embodiment takes nine points as an example, and the smooth input points in the x-axis direction are expressed as follows:
[0092]
[0093] Where (i, j, k) is the coordinate of a point in the underground stratum, which is used as the first smoothing input point.
[0094] (i-1,j,k-ceil{1×tan[θ xz (i,j,k)]}) is the second smooth input point in the x-axis direction,
[0095] (i-2,j,k-ceil{2×tan[θ xz (i,j,k)]}) is the third smooth input point in the x-axis direction,
[0096] (i+1,j,k+ceil{1×tan[θ xz (i,j,k)]}) is the fourth smooth input point in the x-axis direction,
[0097] (i+2,j,k+ceil{2×tan[θ xz (i,j,k)]}) is the fifth smooth input point in the x-axis direction, tan[θ xz (i, j, k)] is the tangent function value of the inclination angle on the plane formed by the x-axis and the z-axis, and ceil{} means rounding the content in {} to obtain the smallest integer greater than or equal to the specified expression.
[0098] The smooth input point in the y-axis direction is calculated as follows:
[0099]
[0100] Where (i, j, k) is the coordinate of a point in the underground stratum, which is used as the first smoothing input point.
[0101] (i,j-2,k-ceil{2×tan[θ yz (i,j,k)]}) is the second smooth input point in the y-axis direction,
[0102] (i-2,j,k-ceil{2×tan[θ xz (i,j,k)]}) is the third smooth input point in the y-axis direction,
[0103] (i,j+1,k+ceil{1×tan[θ yz(i,j,k)]}) is the fourth smooth input point in the y-axis direction,
[0104] (i,j+2,k+ceil{2×tan[θ yz (i,j,k)]}) is the fifth smooth input point in the y-axis direction, tan[θ yz (i, j, k)] is the tangent function value of the inclination angle on the plane formed by the y-axis and the z-axis, and ceil{} means rounding the content in {} to obtain the smallest integer greater than or equal to the specified expression.
[0105] like Figure 4 As shown, the data to be smoothed is constructed in this embodiment. It can be seen that there are a lot of noises in the data to be smoothed, and the signal-to-noise ratio is low.
[0106] Step S4: using Gaussian coefficients to perform Gaussian weighted smoothing on the smoothed input points to obtain smoothed results of the smoothed input points. The smoothed results of the smoothed input points of each point constitute the smoothed results of the seismic data.
[0107] Assuming that the seismic data of point (i, j, k) is W(i, j, k), the Gaussian weighted multi-point smoothing at point (i, j, k) can be expressed as:
[0108]
[0109] in, Smoothed result of sliding input points.
[0110] The Gaussian coefficient is:
[0111]
[0112] Among them, G(i,j,k) is the Gaussian coefficient of point (i,j,k), σ is the standard deviation of the smoothed input point, l is the window width, Traverse all points in the three-dimensional seismic data to obtain the smoothing results of the smoothing input points of each point, and finally form the smoothing results of the seismic data.
[0113] like Figure 5 As shown, it is the use of Figure 1 and Figure 2 From the slope-based smoothed data body calculated based on the slope information, we can see that the data is smoother, the noise is effectively removed, the signal-to-noise ratio is improved, and the structural morphology is effectively maintained.
[0114] The present embodiment proposes a seismic data smoothing method based on slope information, which calculates the inverse tangent function value of the slope information to obtain the inclination angle of the corresponding point; constructs nine smoothing input points required for smoothing according to the inclination angle; uses Gaussian coefficients to perform Gaussian weighted smoothing on the smoothing input points to obtain smoothing results of the smoothing input points, and the smoothing results of the smoothing input points of each point constitute the seismic data smoothing results. The method of the present embodiment can achieve a good smoothing effect without destroying the original seismic data information.
[0115] Embodiment 3
[0116] This embodiment proposes a seismic data smoothing device based on slope information, such as Figure 6 As shown, it includes: an information acquisition module, a dip calculation module, a smooth point construction module and a Gaussian weighting module; the information acquisition module is connected to the dip calculation module, the dip calculation module is connected to the smooth point construction module, and the smooth point construction module is connected to the Gaussian weighting module;
[0117] An information acquisition module, used for acquiring slope information for each point in the three-dimensional seismic data;
[0118] The inclination angle calculation module is used to calculate the inverse tangent function value of the slope information to obtain the inclination angle of the corresponding point;
[0119] A smoothing point construction module, used for constructing smoothing input points required for smoothing according to the inclination angle;
[0120] The Gaussian weighting module is used to perform Gaussian weighted smoothing processing on the smoothing input points using Gaussian coefficients to obtain smoothing results of the smoothing input points. The smoothing results of the smoothing input points of each point constitute the smoothing results of the seismic data.
[0121] The present embodiment proposes a seismic data smoothing device based on slope information, which uses an information acquisition module to acquire the slope information of each point in the three-dimensional seismic data; uses a dip calculation module to calculate the arctangent function value of the slope information to obtain the dip of the corresponding point; uses a smooth point construction module to construct the smooth input point required for smoothing according to the dip; uses a Gaussian weighting module to perform Gaussian weighted smoothing on the smooth input point using a Gaussian coefficient to obtain the smoothing result of the smooth input point, and the smoothing results of the smooth input points of each point constitute the seismic data smoothing result. The present embodiment can achieve a better smoothing effect without destroying the original seismic data information.
[0122] Embodiment 4:
[0123] This embodiment provides an electronic device, including: a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the seismic data smoothing method based on slope information is executed.
[0124] The electronic device may be a mobile phone, a computer or a tablet computer, etc., including a memory and a processor, wherein a computer program is stored in the memory, and when the computer program is executed by the processor, the seismic data smoothing method based on slope information as described in the embodiment is implemented. It is understood that the electronic device may also include an input / output (I / O) interface and a communication component.
[0125] The processor is used to execute all or part of the steps in the seismic data smoothing method based on slope information as described in the above embodiment. The memory is used to store various types of data, which may include instructions of any application or method in the electronic device, as well as data related to the application.
[0126] The processor can be an application specific integrated circuit (ASIC), a digital signal processor (DSP), a programmable logic device (PLD), a field programmable gate array (FPGA), a controller, a microcontroller, a microprocessor or other electronic components, and is used to execute the seismic data smoothing method based on slope information described in the above embodiments.
[0127] Embodiment 5
[0128] This embodiment provides a computer-readable storage medium storing executable instructions, which, when executed, enable a processor to perform the seismic data smoothing method based on slope information.
[0129] If implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium.
[0130] Based on this understanding, the technical solution of the present application, or the part that contributes to the prior art, or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes a number of instructions for enabling a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the seismic data smoothing method based on slope information described in each embodiment of the present application.
[0131] The aforementioned storage media include: flash memory, hard disk, multimedia card, card-type memory (for example, SD (Secure Digital Memory Card) or DX (Memory Data Register, MDR abbreviation, memory data register) memory, etc.), random access memory (RAM), static random access memory (SRAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), programmable read-only memory (PROM), magnetic memory, disk, CD, server, APP (Application, abbreviation of application software) application store and other media that can store program verification codes, on which computer programs are stored. When the computer program is executed by the processor, the various steps of the above-mentioned seismic data smoothing method based on slope information can be implemented.
[0132] The various embodiments in the present disclosure are described in a progressive manner, and the same or similar parts between the various embodiments can be referenced to each other, and each embodiment focuses on the differences from other embodiments.
[0133] The protection scope of the present disclosure is not limited to the above-mentioned embodiments. Obviously, those skilled in the art can make various changes and modifications to the present disclosure without departing from the scope and spirit of the present disclosure. If these changes and modifications fall within the scope of the claims of the present disclosure and their equivalents, the intention of the present disclosure also includes these changes and modifications.
Claims
1. A seismic data smoothing method based on slope information, characterized in that: include: Step S1: Obtain slope information for each point in the three-dimensional seismic data; Step S2: Calculate the inverse tangent function value of the slope information to obtain the inclination angle of the corresponding point; Step S3: constructing smoothing input points required for smoothing according to the inclination angle; Step S4: using Gaussian coefficients to perform Gaussian weighted smoothing on the smoothed input points to obtain smoothed results of the smoothed input points. The smoothed results of the smoothed input points of each point constitute the smoothed results of the seismic data.
2. The seismic data smoothing method based on slope information according to claim 1, characterized in that: The x-axis direction of the three-dimensional seismic data is data in the direction of the seismic main survey line, the y-axis direction is data in the direction of the seismic contact survey line, and the z-axis direction is data in the direction of the seismic depth; The step of obtaining slope information for each point in the three-dimensional seismic data includes: For each point in the three-dimensional seismic data, obtain slope information on the plane formed by the x-axis and the z-axis; For each point in the 3D seismic data, the slope information on the plane formed by the y-axis and the z-axis is obtained.
3. The seismic data smoothing method based on slope information according to claim 2, characterized in that: The calculating the inverse tangent function value of the slope information to obtain the inclination angle of the corresponding point includes: Using the slope information on the plane formed by the x-axis and the z-axis, calculate the inclination angle of the corresponding point on the plane formed by the x-axis and the z-axis; Using the slope information on the plane formed by the y-axis and the z-axis, the inclination angle of the corresponding point on the plane formed by the y-axis and the z-axis is calculated. The calculation formula is as follows: Among them, S xz (i, j, k) is the slope information on the plane formed by the x-axis and the z-axis, S yz (i, j, k) is the slope information on the plane formed by the y-axis and the z-axis, θ xz (i,j,k) is the inclination angle of point (i,j,k) on the plane formed by the x-axis and the z-axis, θ yz (i,j,k) is the inclination angle of point (i,j,k) on the plane formed by the y-axis and the z-axis, and (i,j,k) is the coordinates of a point in the underground stratum.
4. The seismic data smoothing method based on slope information according to claim 3, characterized in that: The step of constructing a smoothing input point required for smoothing according to the inclination angle comprises: Obtaining a smooth input point in the x-axis direction according to the inclination angle of the corresponding point on the plane formed by the x-axis and the z-axis; According to the inclination angle of the corresponding point on the plane formed by the y-axis and the z-axis, a smooth input point in the y-axis direction is obtained.
5. The seismic data smoothing method based on slope information according to claim 4, characterized in that: The smooth input point in the x-axis direction is expressed as follows: Where N is the number of smoothing input points, N-1 is an integer multiple of 4, (i, j, k) is the coordinate of a point in the underground stratum, as the first smoothing input point, (i-1, j, k-ceil{1×tan[θ xz (i,j,k)]}) is the second smooth input point in the x-axis direction, is the first Smooth input point, (i+1,j,k+ceil{1×tan[θ xz (i,j,k)]}) is the first smooth input points, is the first smooth input points, tan[θ xz (i, j, k)] is the tangent function value of the inclination angle on the plane formed by the x-axis and the z-axis, and ceil{} means rounding the content in {}.
6. The seismic data smoothing method based on slope information according to claim 4, characterized in that: The smooth input point in the y-axis direction is calculated as follows: Where N is the number of smoothing input points, N-1 is an integer multiple of 4, (i, j, k) is the coordinate of a point in the underground stratum, which is used as the first smoothing input point. is the second smooth input point in the y-axis direction, is the first Smooth input point, (i,j+1,k+ceil{1×tan[θ yz (i,j,k)]}) is the first smooth input points, is the first smooth input points, tan[θ yz (i, j, k)] is the tangent function value of the inclination angle on the plane formed by the y-axis and the z-axis, and ceil{} means rounding the content in {}.
7. The seismic data smoothing method based on slope information according to claim 4, characterized in that: The Gaussian coefficient is: Among them, G(i,j,k) is the Gaussian coefficient of point (i,j,k), σ is the standard deviation of the smoothed input point, l is the window width, 8. A seismic data smoothing device based on slope information, characterized in that: include: An information acquisition module, used for acquiring slope information for each point in the three-dimensional seismic data; The inclination angle calculation module is used to calculate the inverse tangent function value of the slope information to obtain the inclination angle of the corresponding point; A smoothing point construction module, used for constructing smoothing input points required for smoothing according to the inclination angle; The Gaussian weighting module is used to perform Gaussian weighted smoothing processing on the smoothing input points using Gaussian coefficients to obtain smoothing results of the smoothing input points. The smoothing results of the smoothing input points of each point constitute the smoothing results of the seismic data.
9. An electronic device, characterized in that: include: A memory and a processor, wherein a computer program is stored in the memory, and when the computer program is executed by the processor, the seismic data smoothing method based on slope information as described in any one of claims 1 to 7 is executed.
10. A computer-readable storage medium, characterized in that: It stores executable instructions, which, when executed, enable the processor to execute the seismic data smoothing method based on slope information as described in any one of claims 1-7.
Citation Information
Patent Citations
Smoothing optimization method of velocity model
CN103076627A
Speed smoothing method and system of construction constraint
CN108693556A
Vertical line deviation calculation method and system based on interference imaging altimeter
CN113900069A
Automatic Tracking of Faults by Slope Decomposition
US20150117151A1