A seismic data smoothing method and device based on slope information
By obtaining the slope information of 3D seismic data, calculating the dip and performing Gaussian weighted smoothing, the problem of conventional smoothing methods destroying data information is solved, and effective data smoothing and information retention are achieved.
Patent Information
- Application Number
- CN202311126137.7
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2023-09-01
- Publication Date
- 2025-10-17
- Estimated Expiration
- 2043-09-01
AI Technical Summary
When processing seismic data, conventional smoothing methods can effectively remove noise, but they will destroy the original data information, especially the structural details.
By obtaining the slope information in the 3D seismic data, calculating the inverse tangent function value to get the dip angle, constructing the smooth input points, and using the Gaussian coefficient to perform Gaussian weighted smoothing, a seismic data smoothing method based on slope information is formed.
An effective data smoothing effect is achieved while maintaining the information integrity of the original seismic data, especially the details of the structural morphology.
Smart Images

Figure CN119556346B_ABST
Abstract
Description
TECHNICAL FIELD
[0001] The application belongs to the technical field of geophysical exploration, and particularly relates to a slope information-based seismic data smoothing method and device. BACKGROUND
[0002] For actual collected seismic data, noise, abnormal signals and the like often exist, which leads to a decline in data quality and makes it difficult to obtain good results through processing. At this time, data smoothing technology is usually used to suppress or eliminate the influence of noise, abnormal signals and the like, so that the data is more regular. However, excessive smoothing can destroy much information in the original data and blur details in the signals. The conventional smoothing method does not consider the structural morphology, but directly uses a multi-point smoothing operator to perform filtering processing along three dimensions. Although a good smoothing effect can be achieved, the original data information is destroyed. SUMMARY
[0003] In view of the above technical problems, the application provides a slope information-based seismic data smoothing method and device.
[0004] In a first aspect, the application provides a slope information-based seismic data smoothing method, which comprises the following steps:
[0005] Step S1: acquiring slope information for each point in three-dimensional seismic data;
[0006] Step S2: calculating an inverse tangent function value of the slope information to obtain an inclination angle of the corresponding point;
[0007] Step S3: constructing a smoothing input point required for smoothing according to the inclination angle;
[0008] Step S4: performing Gaussian weighted smoothing processing on the smoothing input point by using a Gaussian coefficient to obtain a smoothing result of the smoothing input point, and the smoothing result of the smoothing input point of each point constitutes a seismic data smoothing result.
[0009] The x-axis direction of the three-dimensional seismic data is data in a main seismic line direction, the y-axis direction is data in a seismic contact line direction, and the z-axis direction is data in a seismic depth direction;
[0010] The acquiring of the slope information for each point in the three-dimensional seismic data comprises the following steps:
[0011] For each point in the three-dimensional seismic data, slope information on a plane composed of the x-axis and the z-axis is acquired;
[0012] For each point in the three-dimensional seismic data, slope information on a plane composed of the y-axis and the z-axis is acquired.
[0013] The arctangent function value of the slope information is calculated to obtain the dip angle of the corresponding point, including:
[0014] The slope information on the plane composed of the x-axis and the z-axis is used to calculate the dip angle of the corresponding point on the plane composed of the x-axis and the z-axis;
[0015] The slope information on the plane composed of the y-axis and the z-axis is used to calculate the dip angle of the corresponding point on the plane composed of the y-axis and the z-axis, and the calculation formula is as follows:
[0016]
[0017] Wherein, S xz (i,j,k) is the slope information on the plane composed of the x-axis and the z-axis, S yz (i,j,k) is the slope information on the plane composed of the y-axis and the z-axis, θ xz (i,j,k) is the dip angle of the point (i,j,k) on the plane composed of the x-axis and the z-axis, θ yz (i,j,k) is the dip angle of the point (i,j,k) on the plane composed of the y-axis and the z-axis, and (i,j,k) is the coordinates of a point in the underground stratum.
[0018] The dip angle is used to construct the smoothing input point required for smoothing, including:
[0019] The dip angle of the corresponding point on the plane composed of the x-axis and the z-axis is used to obtain the smoothing input point in the x-axis direction;
[0020] The dip angle of the corresponding point on the plane composed of the y-axis and the z-axis is used to obtain the smoothing input point in the y-axis direction.
[0021] The expression of the smoothing input point in the x-axis direction is as follows:
[0022]
[0023] Wherein, N is the number of smoothing input points, N-1 is an integer multiple of 4, (i,j,k) is the coordinates 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 smoothing input point in the x-axis direction, is the first smoothing input point in the x-axis direction, is the second smoothing input point in the x-axis direction, xz is the first smoothing input point in the x-axis direction, is the second smoothing input point in the x-axis direction, is the first smoothing input point in the x-axis direction, is the second smoothing input point in the x-axis direction, and tan[θxz (i,j,k)] is the tangent function value of the dip angle in the plane composed of the x-axis and the z-axis, and ceil{} represents rounding the content in {} to an integer.
[0024] The smoothing input point in the y-axis direction is calculated as follows:
[0025]
[0026] wherein 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 subsurface formation, as the first smoothing input point, is the second smoothing input point in the y-axis direction, is the third smoothing input point in the y-axis direction, is the fourth smoothing input point in the y-axis direction, yz (i,j,k)] is the tangent function value of the dip angle in the plane composed of the x-axis and the z-axis, and ceil{} represents rounding the content in {} to an integer. is the fourth smoothing input point in the y-axis direction, is the fifth smoothing input point in the y-axis direction, is the fifth smoothing input point in the y-axis direction, yz (i,j,k)] is the tangent function value of the dip angle in the plane composed of the x-axis and the z-axis, and ceil{} represents rounding the content in {} to an integer.
[0027] The Gaussian coefficient is:
[0028]
[0029] wherein G(i,j,k) is the Gaussian coefficient of the point (i,j,k), σ is the standard deviation of the smoothing input points, and l represents the window width,
[0030] In a second aspect, the application provides a device for smoothing seismic data based on slope information, comprising:
[0031] An information acquisition module is configured to acquire slope information for each point in three-dimensional seismic data.
[0032] An inclination calculation module is configured to calculate the arctangent function value of the slope information to obtain the dip angle of the corresponding point.
[0033] A smoothing point construction module is configured to construct smoothing input points required for smoothing according to the dip angle.
[0034] A Gaussian weighting module is configured to perform Gaussian weighting smoothing processing on the smoothing input points using a Gaussian coefficient to obtain the smoothing result of the smoothing input points, and the smoothing result of each point constitutes the smoothing result of the seismic data.
[0035] In a third aspect, the present application provides an electronic device, comprising a memory and a processor, wherein the memory stores a computer program, and the computer program is executed by the processor to implement the slope information based seismic data smoothing method.
[0036] In a fourth aspect, the present application provides a computer readable storage medium, which stores executable instructions, and the instructions, when executed, cause a processor to implement the slope information based seismic data smoothing method.
[0037] Advantages:
[0038] The present application provides a slope information based seismic data smoothing method and device, which realizes the slope information based seismic data smoothing method, can achieve good smoothing effect, and will not damage the original seismic data information. BRIEF DESCRIPTION OF DRAWINGS
[0039] Figure 1 A slope information based seismic data smoothing method flow chart is provided for the present application;
[0040] Figure 2 An x direction stratum slope schematic diagram is provided for the present application;
[0041] Figure 3 A y direction stratum slope schematic diagram is provided for the present application;
[0042] Figure 4 A data schematic diagram to be smoothed obtained by construction is provided for the present application;
[0043] Figure 5 A slope based smoothed data body schematic diagram is provided for the present application;
[0044] Figure 6 A slope information based seismic data smoothing device principle block diagram is provided for the present application. DETAILED DESCRIPTION
[0045] The present disclosure will be further described below in conjunction with the embodiments shown in the 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, so that the data is more regular. However, excessive smoothing will damage much information in the original data and blur the details in the signal. The prior art does not consider the construction form, but directly uses a multi-point smoothing operator to perform filtering processing along three dimensions, which can achieve good smoothing effect, but will damage the original data information.
[0047] The application provides a slope information-based seismic data smoothing method, which comprises the following steps: firstly, reading slope information of each point in three-dimensional seismic data; then, calculating and extracting corresponding smoothing input points along the slope information; and finally, performing multi-point smoothing on the smoothing input points under slope constraint by using Gaussian coefficients, so as to realize the slope information-based data smoothing method.
[0048] Embodiment one,
[0049] The embodiment provides a slope information-based seismic data smoothing method, which comprises the following steps as shown in the figure: Figure 1
[0050] Step S1: acquiring slope information of each point in three-dimensional seismic data;
[0051] In the embodiment, the x-axis direction of the three-dimensional seismic data is data in a main seismic line direction, the y-axis direction is data in a seismic contact line direction, and the z-axis direction is data in a seismic depth direction.
[0052] The acquiring of the slope information of each point in the three-dimensional seismic data comprises the following steps:
[0053] For each point in the three-dimensional seismic data, slope information on a plane composed of the x-axis and the z-axis is acquired.
[0054] For each point in the three-dimensional seismic data, slope information on a plane composed of the y-axis and the z-axis is acquired.
[0055] Step S2: calculating an inverse tangent function value of the slope information to obtain a dip angle of the corresponding point, which comprises the following steps:
[0056] For each point in the three-dimensional seismic data, slope information on a plane composed of the x-axis and the z-axis is acquired.
[0057] For each point in the three-dimensional seismic data, slope information on a plane composed of the y-axis and the z-axis is acquired.
[0058]
[0059] Wherein, S xz (i,j,k) is the slope information on the plane composed of the x-axis and the z-axis, S yz (i,j,k) is the slope information on the plane composed of the y-axis and the z-axis, θ xz (i,j,k) is the dip angle of the point (i,j,k) on the plane composed of the x-axis and the z-axis, θ yz (i,j,k) is the dip angle of point (i,j,k) on the plane composed of y-axis and z-axis, (i,j,k) is the coordinate of a point in the subsurface formation.
[0060] Step S3: constructing the smoothing input points required for smoothing according to the dip angle, including:
[0061] obtaining the smoothing input points in the x-axis direction according to the dip angle of the corresponding point on the plane composed of x-axis and z-axis;
[0062] obtaining the smoothing input points in the y-axis direction according to the dip angle of the corresponding point on the plane composed of y-axis and z-axis.
[0063] The smoothing input points in the x-axis direction are expressed as follows:
[0064]
[0065] wherein 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 subsurface formation, (i-1,j,k-ceil{1xtan[θ xz (i,j,k)]}) is the second smoothing input point in the x-axis direction, is the first smoothing input point in the x-axis direction, is the second smoothing input point in the x-axis direction, xz is the first smoothing input point in the x-axis direction, is the second smoothing input point in the x-axis direction, is the first smoothing input point in the x-axis direction, is the second smoothing input point in the x-axis direction, xz tan[θ yz (i,j,k)] is the tangent value of the dip angle on the plane composed of x-axis and z-axis, and ceil{} means rounding the content in {} to an integer.
[0066] The smoothing input points in the y-axis direction are calculated as follows:
[0067]
[0068] wherein 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 subsurface formation, (i-1,j,k-ceil{1xtan[θ (i,j,k)]}) is the second smoothing input point in the y-axis direction, is the first smoothing input point in the y-axis direction, is the second smoothing input point in the y-axis direction, yz is the first smoothing input point in the y-axis direction, is the second smoothing input point in the y-axis direction, for the first smooth input point in the y-axis direction tan [θ(i,j,k)] is the tangent function value of the inclination angle in the plane composed of the y-axis and the z-axis, and ceil{} represents rounding the content in {} to an integer. yz
[0069] Step S4: Gaussian weighting smoothing is performed on the smooth input points by using the Gaussian coefficient, to obtain the smoothing result of the smooth input points, and the smoothing result of the smooth input points of each point constitutes the seismic data smoothing result.
[0070] Suppose the seismic data of point (i,j,k) is W(i,j,k), then the Gaussian weighting multi-point smoothing at point (i,j,k) can be expressed as:
[0071]
[0072] wherein, the smoothing result of the smooth input points.
[0073] The Gaussian coefficient is:
[0074]
[0075] wherein, G(i,j,k) is the Gaussian coefficient of point (i,j,k), σ is the standard deviation of the smooth input points, and l represents the window width, All points in the three-dimensional seismic data are traversed to obtain the smoothing result of the smooth input points of each point, and finally the seismic data smoothing result is constituted.
[0076] The method for smoothing seismic data based on slope information provided in this embodiment first acquires the slope information of each point, calculates the arctangent function value of the slope information to obtain the inclination angle of the corresponding point, constructs the smooth input points required for smoothing according to the inclination angle, performs Gaussian weighting smoothing on the smooth input points by using the Gaussian coefficient to obtain the smoothing result of the smooth input points, and the smoothing result of the smooth input points of each point constitutes the seismic data smoothing result. The method of this embodiment can achieve good smoothing effect and will not damage the original seismic data information.
[0077] Embodiment two,
[0078] This embodiment is a specific example of embodiment one, real three-dimensional seismic data is input, and nine-point smoothing is taken as an example to describe the method for smoothing seismic data based on slope information in detail, as shown in FIG. 1, which includes the following steps: Figure 1
[0079] Step S1: Acquire the slope information of 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, 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 The figure shows the y-direction formation slope input in this embodiment, which is used to calculate the specific location of the smooth input point required in the y-direction. The y-direction formation slope is the slope information on the plane formed by the y-axis and the z-axis. The slope information is obtained for each point in the 3D seismic data, including:
[0081] For each point in the 3D 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 as follows:
[0086]
[0087] Among them, S xz (i, j, k) is the slope information on the plane composed of 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 coordinate 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 of the corresponding point on the plane composed of the y-axis and the z-axis, the smooth input point in the y-axis direction is obtained.
[0091] The embodiment takes nine points as an example, and the smooth input point in the x-axis direction is expressed as follows:
[0092]
[0093] wherein (i, j, k) is the coordinate of a point in the underground stratum, as the first smooth input point,
[0094] (i-1, j, k-ceil{1x tan[θ xz (i, j, k)]}) is the second smooth input point in the x-axis direction,
[0095] (i-2, j, k-ceil{2x tan[θ xz (i, j, k)]}) is the third smooth input point in the x-axis direction,
[0096] (i+1, j, k+ceil{1x tan[θ xz (i, j, k)]}) is the fourth smooth input point in the x-axis direction,
[0097] (i+2, j, k+ceil{2x tan[θ xz (i, j, k)]}) is the fifth smooth input point in the x-axis direction, tan[θ xz (i, j, k)] is the tangent value of the inclination on the plane composed of the x-axis and the z-axis, and ceil{} means rounding the content in {} to the minimum 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] wherein (i, j, k) is the coordinate of a point in the underground stratum, as the first smooth input point,
[0101] (i, j-2, k-ceil{2x tan[θ yz (i, j, k)]}) is the second smooth input point in the y-axis direction,
[0102] (i-2, j, k-ceil{2x tan[θ xz (i, j, k)]}) is the third smooth input point in the y-axis direction,
[0103] (i, j+1, k+ceil{1x tan[θ yz(i,j,k)]}) is the fourth smoothing input point in the y-axis direction,
[0104] (i,j+2,k+ceil{2×tan[θ yz (i,j,k)]}) is the fifth smoothing input point in the y-axis direction, tan[θ yz (i,j,k)] is the tangent value of the inclination angle in the plane composed of the y-axis and the z-axis, and ceil{} represents rounding the content in {} to the minimum integer greater than or equal to the specified expression.
[0105] As Figure 4 shown, it is the data to be smoothed constructed in the embodiment, and it can be seen that there are many noises in the data to be smoothed, and the signal-to-noise ratio is low.
[0106] Step S4: Gaussian weighting smoothing processing is performed on the smoothing input points by using the Gaussian coefficient, to obtain the smoothing result of the smoothing input points, and the smoothing result of each point constitutes the seismic data smoothing result.
[0107] Suppose that the seismic data of point (i,j,k) is W(i,j,k), then the Gaussian weighting multi-point smoothing at point (i,j,k) can be expressed as:
[0108]
[0109] wherein, the smoothing result of the sliding input point.
[0110] The Gaussian coefficient is:
[0111]
[0112] wherein, G(i,j,k) is the Gaussian coefficient of point (i,j,k), σ is the standard deviation of the smoothing input point, and l represents the window width, All points in the three-dimensional seismic data are traversed to obtain the smoothing result of the smoothing input point of each point, and finally the seismic data smoothing result is constituted.
[0113] As Figure 5 shown, it is the data to be smoothed constructed in the embodiment, and it can be seen that there are many noises in the data to be smoothed, and the signal-to-noise ratio is low. Figure 1 and Figure 2 the slope information, and it can be seen that the data is smoother, the noise is effectively removed, the signal-to-noise ratio is improved, and the structure is effectively maintained.
[0114] The embodiment provides a slope information-based seismic data smoothing method, an inverse tangent function value of slope information is calculated to obtain an inclination angle of a corresponding point; nine smoothing input points required for smoothing are constructed according to the inclination angle; Gaussian weighting smoothing processing is performed on the smoothing input points by using Gaussian coefficients to obtain a smoothing result of the smoothing input points, and the smoothing result of the smoothing input points of each point constitutes a seismic data smoothing result. The method can achieve a good smoothing effect and does not damage original seismic data information.
[0115] Embodiment three,
[0116] The embodiment provides a slope information-based seismic data smoothing device, as shown in the figure, which comprises an information acquisition module, an inclination angle calculation module, a smoothing point construction module and a Gaussian weighting module; the information acquisition module is connected with the inclination angle calculation module, the inclination angle calculation module is connected with the smoothing point construction module, and the smoothing point construction module is connected with the Gaussian weighting module. Figure 6
[0117] The information acquisition module is used for acquiring slope information of each point in three-dimensional seismic data.
[0118] The inclination angle calculation module is used for calculating an inverse tangent function value of the slope information to obtain an inclination angle of a corresponding point.
[0119] The smoothing point construction module is used for constructing smoothing input points required for smoothing according to the inclination angle.
[0120] The Gaussian weighting module is used for performing Gaussian weighting smoothing processing on the smoothing input points by using Gaussian coefficients to obtain a smoothing result of the smoothing input points, and the smoothing result of the smoothing input points of each point constitutes a seismic data smoothing result.
[0121] The embodiment provides a slope information-based seismic data smoothing device, which acquires slope information of each point in three-dimensional seismic data by using an information acquisition module; calculates an inverse tangent function value of the slope information to obtain an inclination angle of a corresponding point by using an inclination angle calculation module; constructs smoothing input points required for smoothing according to the inclination angle by using a smoothing point construction module; and performs Gaussian weighting smoothing processing on the smoothing input points by using a Gaussian weighting module to obtain a smoothing result of the smoothing input points, and the smoothing result of the smoothing input points of each point constitutes a seismic data smoothing result. The embodiment can achieve a good smoothing effect and does not damage original seismic data information.
[0122] Embodiment four,
[0123] The embodiment provides an electronic device, which comprises a memory and a processor, and the memory stores a computer program which is executed by the processor to execute the slope information-based seismic data smoothing method.
[0124] The electronic device can be a mobile phone, a computer, a tablet computer, or the like, and includes a memory and a processor. The memory stores a computer program which, when executed by the processor, implements the slope information based seismic data smoothing method as described in the embodiments. It can be understood that the electronic device can further include an input / output (I / O) interface and a communication component.
[0125] The processor is configured to execute all or part of the steps of the slope information based seismic data smoothing method as described in the embodiments. The memory is configured to store various types of data, which can include, for example, instructions of any application program or method in the electronic device, and application program related data.
[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 elements, and is configured to execute the slope information based seismic data smoothing method as described in the embodiments.
[0127] Embodiment five,
[0128] The embodiment provides a computer readable storage medium storing executable instructions, which, when executed, cause a processor to execute the slope information based seismic data smoothing method.
[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 such understanding, the technical solutions of the present application, in essence, or the part that contributes to the prior art, or part of the technical solutions 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 causing 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 slope information based seismic data smoothing method described in the embodiments.
[0131] The aforementioned storage medium includes a flash memory, a hard disk, a multimedia card, a card-type memory (e.g., an SD (Secure Digital Memory Card) or a DX (Memory Data Register, MDR) memory, etc.), a random access memory (RAM), a static random access memory (SRAM), a read-only memory (ROM), an electrically erasable programmable read-only memory (EEPROM), a programmable read-only memory (PROM), a magnetic memory, a magnetic disk, an optical disk, a server, an APP (Application) application market, etc., which can store a program check code, and stores a computer program thereon, which, when executed by a processor, can implement each step of the seismic data smoothing method based on slope information described above.
[0132] Each of the embodiments in the present disclosure is described in a progressive manner, and the same or similar parts between the embodiments can be referred to each other. Each of the embodiments focuses on the difference from other embodiments.
[0133] The scope of protection of the present disclosure is not limited to the above-described embodiments, and it is obvious that those skilled in the art can make various modifications and changes to the present disclosure without departing from the scope and spirit of the present disclosure. If these modifications and changes belong to the scope of the claims of the present disclosure and equivalent technologies thereof, the intention of the present disclosure also includes these modifications and changes.
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 3D seismic data; 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; The obtaining of slope information for each point in the three-dimensional seismic data includes: obtaining slope information on a plane formed by an x-axis and a z-axis for each point in the three-dimensional seismic data; obtaining slope information on a plane formed by a y-axis and a z-axis 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 seismic data smoothed results.
2. The seismic data smoothing method based on slope information according to claim 1, characterized in that: 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 as follows: in, is the slope information on the plane formed by the x-axis and the z-axis, is the slope information on the plane formed by the y-axis and the z-axis, for point The inclination angle on the plane formed by the x-axis and the z-axis, for point The inclination angle on the plane formed by the y-axis and the z-axis, are the coordinates of a point in the underground stratum.
3. The seismic data smoothing method based on slope information according to claim 2, characterized in that: The step of constructing smoothing input points required for smoothing according to the inclination angle includes: 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.
4. The seismic data smoothing method based on slope information according to claim 3, characterized in that: The smooth input point in the x-axis direction is expressed as follows: Among them, N is the number of smoothing input points, N-1 is an integer multiple of 4, is the coordinate of a point in the underground stratum, which serves as the first smoothing input point, is the second smooth input point in the x-axis direction, In the x-axis direction Smooth input points, is the first Smooth input points, In the x-axis direction Smooth input points, is the tangent function value of the inclination angle on the plane formed by the x-axis and the z-axis, Express The contents in are rounded.
5. The seismic data smoothing method based on slope information according to claim 3, characterized in that: The smooth input point in the y-axis direction is calculated as follows: Among them, N is the number of smoothing input points, N-1 is an integer multiple of 4, is the coordinate of a point in the underground stratum, which serves as the first smoothing input point, is the second smooth input point in the y-axis direction, In the y-axis direction Smooth input points, is the first Smooth input points, In the y-axis direction Smooth input points, is the tangent function value of the inclination angle on the plane formed by the y-axis and the z-axis, Express The contents in are rounded.
6. The seismic data smoothing method based on slope information according to claim 3, characterized in that: The Gaussian coefficient is: in, for Gaussian coefficient of the point, the standard deviation of the smoothed input points, l Indicates the window width, .
7. A seismic data smoothing device based on slope information, characterized in that: include: An information acquisition module, for acquiring slope information for each point in the three-dimensional seismic data; 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; The obtaining of slope information for each point in the three-dimensional seismic data includes: obtaining slope information on a plane formed by an x-axis and a z-axis for each point in the three-dimensional seismic data; obtaining slope information on a plane formed by a y-axis and a z-axis 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 seismic data smoothing results.
8. 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 according to any one of claims 1 to 6 is executed.
9. A computer-readable storage medium, characterized in that The device stores executable instructions, which, when executed, enable a processor to execute the seismic data smoothing method based on slope information as described in any one of claims 1 to 6.
Citation Information
Patent Citations
Speed smoothing method and system of construction constraint
CN108693556A
Automatic Tracking of Faults by Slope Decomposition
US20150117151A1