A method and device for extracting construction information based on threshold constraints

Abstract

The application provides a method and device for extracting structural information based on threshold constraint. The method comprises the following steps: obtaining three-dimensional seismic data; taking a point of the three-dimensional seismic data, and calculating a first weighted difference operator of the point in a first direction and a second weighted difference operator of the point in a second direction; obtaining initial structural information according to the first weighted difference operator and the second weighted difference operator; calculating a threshold according to the initial structural information; and performing normalization processing on the initial structural information based on the threshold to obtain final structural information. The device comprises a data acquisition module, a difference operator calculation module, an initial structural information calculation module, a threshold calculation module and a final structural information output module. The method provided by the application eliminates the influence of noise or small-scale structures with less contribution to tomographic inversion, and only retains structural information that contributes to tomographic velocity modeling, so that the structural information extraction for depth domain tomographic velocity modeling is realized.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of geophysical exploration, and particularly relates to a structural information extraction method and device based on threshold constraint. BACKGROUND

[0002] At present, a conventional depth domain tomographic velocity modeling method generally does not consider structural information, and an inversion velocity model lacks geological significance and has low inversion precision. In view of the problem, how to add geological structural information for constraint in the process of tomographic velocity modeling is the key to solving the problem. In the technical process, how to efficiently and high-quality extract structural information is the first problem to be solved. At present, a method for extracting structures is to obtain structural information in a slope or local energy stacking manner. The method is generally cumbersome to calculate, and the calculation result is not stable enough. The extracted structural information cannot be quantitatively controlled by a certain parameter, and needs to be additionally modified and modified by using filtering and other technologies. In the actual application process, considering the processing period and complexity and other problems, the method is less applied. SUMMARY

[0003] Based on the above technical problems, the application provides a structural information extraction method and device based on threshold constraint.

[0004] In a first aspect, the application provides a structural information extraction method based on threshold constraint, comprising:

[0005] acquiring three-dimensional seismic data;

[0006] taking a point of the three-dimensional seismic data, calculating a first weighted difference operator of the point in a first direction and a second weighted difference operator in a second direction;

[0007] obtaining initial structural information according to the first weighted difference operator and the second weighted difference operator;

[0008] calculating a threshold value according to the initial structural information;

[0009] based on the threshold value, performing normalization processing on the initial structural information to obtain final structural information.

[0010] The x-axis direction of the three-dimensional seismic data is data in the seismic main line direction, the y-axis direction is data in the seismic contact line direction, and the z-axis direction is data in the seismic depth direction. A point of the three-dimensional seismic data is represented as W(i,j,k). The first weighted difference operator in the first direction is a weighted difference result of N points adjacent to the point W(i,j,k) in the x-axis direction. The second weighted difference operator in the second direction is a weighted difference result of N points adjacent to the point W(i,j,k) in the y-axis direction, wherein N is an integer multiple of 2.

[0011] The N points adjacent to point W(i,j,k) in the x-axis direction are: W(i-1,j,k)…

[0012] as well as

[0013] The N points adjacent to point W(i,j,k) in the y-axis direction are: W(i,j-1,k)…

[0014] as well as

[0015] The first weighted difference operator is calculated as follows:

[0016]

[0017] Among them, C x (i, j, k) is the first weighted difference operator, q1 is the first set weight, q2 is the first set weight Set the weight, q3 is the Set the weight, W(i-1,j,k) is the first point adjacent to point W(i,j,k) in the negative direction of the x-axis, is the first point adjacent to point W(i,j,k) in the negative direction of the x-axis Points, is the first point adjacent to point W(i,j,k) in the negative direction of the x-axis points, W(i+1,j,k) is the first point adjacent to point W(i,j,k) in the positive direction of the x-axis, is the first point adjacent to point W(i,j,k) in the positive direction of the x-axis Points, is the first point adjacent to point W(i,j,k) in the positive direction of the x-axis points.

[0018] The second weighted difference operator is calculated as follows:

[0019]

[0020] Among them, C y (i, j, k) is the second weighted difference operator, q1 is the first set weight, q2 is the second Set the weight, q3 is the Set the weight, W(i,j-1,k) is the first point adjacent to point W(i,j,k) in the negative direction of the y-axis, is the first point adjacent to point W(i,j,k) in the negative direction of the y-axis Points, is the first point adjacent to point W(i,j,k) in the negative direction of the y-axis points, W(i,j+1,k) is the first point adjacent to point W(i,j,k) in the positive direction of the y-axis, is the first point adjacent to point W(i,j,k) in the positive direction of the y-axis Points, is the first point adjacent to point W(i,j,k) in the positive direction of the y-axis points.

[0021] The initial construction information is obtained according to the first weighted difference operator and the second weighted difference operator, and the calculation formula is as follows:

[0022]

[0023] Among them, S(i,j,k) is the initial construction information, C x (i, j, k) is the first weighted difference operator, C y (i, j, k) is the second weighted difference operator.

[0024] The threshold is calculated as follows:

[0025]

[0026] Where thre is the threshold, S(i, j, k) is the initial structural information, ε is the set control coefficient, n1 is the number of data on the x-axis of the 3D seismic data, n2 is the number of data on the y-axis of the 3D seismic data, and n3 is the number of data on the z-axis of the 3D seismic data.

[0027] The initial structure information is normalized by a threshold to obtain the final structure information, which is calculated as follows:

[0028]

[0029] Among them, S(i,j,k) is the initial construction information, thre is the threshold, The final construction information.

[0030] In a second aspect, the present application proposes a structure information extraction device based on threshold constraints, comprising:

[0031] A data acquisition module, used for acquiring three-dimensional seismic data;

[0032] A differential operator calculation module, configured to take a point of the three-dimensional seismic data and calculate a first weighted differential operator of the point in a first direction and a second weighted differential operator in a second direction;

[0033] An initial structure information calculation module is configured to obtain initial structure information according to the first weighted difference operator and the second weighted difference operator.

[0034] A threshold calculation module is configured to calculate a threshold according to the initial structure information.

[0035] A final structure information output module is configured to normalize the initial structure information by the threshold to obtain final structure information.

[0036] In a third aspect, the present application 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 perform the method for extracting structure information based on threshold constraint.

[0037] 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 perform the method for extracting structure information based on threshold constraint.

[0038] Beneficial effects:

[0039] The present application provides a method and device for extracting structure information based on threshold constraint, which obtains initial structure information by using a weighted difference operator, calculates a threshold based on the initial structure information, and obtains final structure information by threshold constraint. The method of the present application eliminates the influence of noise or small-scale structure which has little contribution to tomographic inversion, and only retains structure information which contributes to tomographic velocity modeling, thereby realizing structure information extraction for deep domain tomographic velocity modeling. BRIEF DESCRIPTION OF DRAWINGS

[0040] Figure 1 A flowchart of the method for extracting structure information based on threshold constraint according to an embodiment of the present application;

[0041] Figure 2 A schematic diagram of three-dimensional seismic data according to an embodiment of the present application;

[0042] Figure 3 A schematic diagram of structure information extracted when the threshold is 200 according to an embodiment of the present application;

[0043] Figure 4 A schematic diagram of structure information extracted when the threshold is 400 according to an embodiment of the present application;

[0044] Figure 5 A schematic diagram of structure information extracted when the threshold is 600 according to an embodiment of the present application;

[0045] Figure 6 A principle block diagram of the device for extracting structure information based on threshold constraint according to an embodiment of the present application. DETAILED DESCRIPTION

[0046] The present disclosure will be further described with reference to the embodiments shown in the drawings.

[0047] The conventional depth domain tomographic velocity modeling method does not generally consider the structural information, the inverted velocity model lacks geological significance, and the inversion accuracy is low. The current method for extracting structural information is to obtain structural information in the form of slope or local energy stack. This method is generally more complicated to calculate, and the calculation result is not stable enough. The extracted structural information cannot be quantitatively controlled by a certain parameter, and additional filtering and other techniques need to be used for modification and modification. In the actual application process, considering the processing period and complexity, the application is less.

[0048] The present application provides a threshold constraint based structural information extraction method. The initial structural information is calculated by using a double-azimuth N-point weighted difference operator. Based on the initial structural information, the threshold is set for screening, the influence of noise or small-scale structures with small contribution to tomographic inversion is removed, and only the structural information contributing to tomographic velocity modeling is left to realize the structural information extraction for depth domain tomographic velocity modeling.

[0049] Embodiment one,

[0050] The present embodiment provides a threshold constraint based structural information extraction method, as shown in Figure 1 , which comprises:

[0051] Step S1: obtaining three-dimensional seismic data;

[0052] In the present embodiment, the three-dimensional seismic data is obtained in a three-dimensional coordinate system with x-axis as the horizontal axis, y-axis as the vertical axis, and z-axis as the vertical axis. The x-axis direction is the data of the seismic main line direction, the y-axis direction is the data of the seismic contact line direction, and the z-axis direction is the data of the seismic depth direction.

[0053] Step S2: taking a point of the three-dimensional seismic data, calculating a first weighted difference operator of the point in a first direction and a second weighted difference operator of the point in a second direction;

[0054] In the present embodiment, the point of the three-dimensional seismic data is represented as W(i,j,k), the first weighted difference operator in the first direction is the weighted difference result of N points adjacent to the point W(i,j,k) in the x-axis direction, and the second weighted difference operator in the second direction is the weighted difference result of N points adjacent to the point W(i,j,k) in the y-axis direction, wherein N is an integer multiple of 2.

[0055] The N points adjacent to the point W(i,j,k) in the x-axis direction are: W(i-1,j,k)

[0056] as well as

[0057] The N points adjacent to point W(i,j,k) in the y-axis direction are: W(i,j-1,k)…

[0058] as well as

[0059] Step S3: Obtaining initial construction information according to the first weighted difference operator and the second weighted difference operator;

[0060] The first weighted difference operator is calculated as follows:

[0061]

[0062] Among them, C x (i, j, k) is the first weighted difference operator, q1 is the first set weight, q2 is the first set weight Set the weight, q3 is the Set the weight, W(i-1,j,k) is the first point adjacent to point W(i,j,k) in the negative direction of the x-axis, is the first point adjacent to point W(i,j,k) in the negative direction of the x-axis Points, is the first point adjacent to point W(i,j,k) in the negative direction of the x-axis points, W(i+1,j,k) is the first point adjacent to point W(i,j,k) in the positive direction of the x-axis, is the first point adjacent to point W(i,j,k) in the positive direction of the x-axis Points, is the first point adjacent to point W(i,j,k) in the positive direction of the x-axis points.

[0063] The second weighted difference operator is calculated as follows:

[0064]

[0065] Among them, C y (i, j, k) is the second weighted difference operator, q1 is the first set weight, q2 is the second Set the weight, q3 is the Set the weight, W(i,j-1,k) is the first point adjacent to point W(i,j,k) in the negative direction of the y-axis, is the first point adjacent to point W(i,j,k) in the negative direction of the y-axis Points, is the first point adjacent to point W(i,j,k) in the negative direction of the y-axis points, W(i,j+1,k) is the first point adjacent to point W(i,j,k) in the positive direction of the y-axis, is the first point adjacent to point W(i,j,k) in the positive direction of the y-axis Points, is the first point adjacent to point W(i,j,k) in the positive direction of the y-axis points.

[0066] The initial construction information is obtained according to the first weighted difference operator and the second weighted difference operator, and the calculation formula is as follows:

[0067]

[0068] Among them, S(i,j,k) is the initial construction information, C x (i, j, k) is the first weighted difference operator, C y (i, j, k) is the second weighted difference operator. From the above formula, it can be seen that the sum of the cube of the first weighted difference operator and the cube of the second weighted difference operator is the initial construction information.

[0069] Step S4: Calculating a threshold value based on the initial construction information;

[0070] The threshold is calculated as follows:

[0071]

[0072] Among them, thre is the threshold, S(i,j,k) is the initial structural information, ε is the set control coefficient, n1 is the number of data on the x-axis of the 3D seismic data, n2 is the number of data on the y-axis of the 3D seismic data, and n3 is the number of data on the z-axis of the 3D seismic data. From the above formula, it can be seen that the data of the initial structural information S(i,j,k) in the x-direction, y-direction, and z-direction are summed and averaged, and then the set control system is used to multiply the average value to obtain the final threshold. The initial structural information is filtered by this threshold to eliminate the influence of noise or small-scale structures that have little contribution to tomographic inversion, leaving only the structural information that contributes to tomographic velocity modeling, thereby realizing structural information extraction for deep-domain tomographic velocity modeling.

[0073] Step S5: Based on the threshold, the initial structure information is normalized to obtain the final structure information, which is calculated as follows:

[0074]

[0075] Among them, S(i,j,k) is the initial construction information, thre is the threshold, to obtain final structural information.

[0076] The embodiment proposes a threshold constraint-based structural information extraction method, which takes a point of three-dimensional seismic data, calculates a first weighted difference operator of the point in a first direction and a second weighted difference operator of the point in a second direction, obtains initial structural information according to the first weighted difference operator and the second weighted difference operator, calculates a threshold according to the initial structural information, and performs normalization processing on the initial structural information based on the threshold to obtain final structural information. The embodiment filters the initial structural information through the threshold, eliminates the influence of noise or small-scale structures with less contribution to tomographic inversion, and only leaves structural information that contributes to tomographic velocity modeling, thereby realizing structural information extraction for depth domain tomographic velocity modeling.

[0077] Embodiment two,

[0078] The embodiment is a specific example of the embodiment one, adopts a dual-azimuth six-point weighted difference operator to calculate initial structural information, then calculates a threshold based on the initial structural information, and filters the initial structural information by setting the threshold. The detailed steps are as follows:

[0079] Step S100: Obtain three-dimensional seismic data.

[0080] In the embodiment, the three-dimensional seismic data is obtained in a three-dimensional coordinate composed of an x-axis as a horizontal axis, a y-axis as a vertical axis, and a z-axis as a vertical axis, wherein the x-axis is data in a seismic main line direction, the y-axis is data in a seismic contact line direction, and the z-axis is data in a seismic depth direction.

[0081] Step S101: Take a point of the three-dimensional seismic data, calculate a first weighted difference operator of the point in a first direction and a second weighted difference operator of the point in a second direction.

[0082] In the embodiment, the point of the three-dimensional seismic data is represented as W(i,j,k), the first weighted difference operator in the first direction is the weighted difference result of six points adjacent to the point W(i,j,k) in the x-axis direction, and the second weighted difference operator in the second direction is the weighted difference result of six points adjacent to the point W(i,j,k) in the y-axis direction.

[0083] The six points adjacent to the point W(i,j,k) in the x-axis direction are W(i-1,j,k), W(i-2,j,k), W(i-3,j,k), W(i+1,j,k), W(i+2,j,k), and W(i+3,j,k).

[0084] The six points adjacent to the point W(i,j,k) in the y-axis direction are: W(i,j-1,k), W(i,j-2,k), W(i,j-3,k), W(i,j+1,k), W(i,j+2,k) and W(i,j+3,k).

[0085] Step S102: obtaining initial construction information according to the first weighted difference operator and the second weighted difference operator;

[0086] The first weighted difference operator is constructed by using the six points adjacent to the point W(i,j,k) in the x-axis direction, and the calculation formula is as follows:

[0087]

[0088] wherein, C x The first weighted difference operator is W(i,j,k), q1 is the first setting weight, the first setting weight in the embodiment is 0.6, q2 is the second setting weight, the second setting weight in the embodiment is 0.2, q3 is the third setting weight, the third setting weight in the embodiment is 0.3, W(i-1,j,k) is the first point adjacent to the point W(i,j,k) in the negative direction of the x-axis, W(i-2,j,k) is the second point adjacent to the point W(i,j,k) in the negative direction of the x-axis, W(i-3,j,k) is the third point adjacent to the point W(i,j,k) in the negative direction of the x-axis, W(i+1,j,k) is the first point adjacent to the point W(i,j,k) in the positive direction of the x-axis, W(i+2,j,k) is the second point adjacent to the point W(i,j,k) in the positive direction of the x-axis, and W(i+3,j,k) is the third point adjacent to the point W(i,j,k) in the positive direction of the x-axis.

[0089] The second weighted difference operator is constructed by using the six points adjacent to the point W(i,j,k) in the y-axis direction, and the calculation formula is as follows:

[0090]

[0091] wherein, C y(i,j,k) is the second weighted difference operator, q1 is the first setting weight, the first setting weight of the embodiment is 0.6, q2 is the second setting weight, the second setting weight of the embodiment is 0.2, q3 is the third setting weight, the third setting weight of the embodiment is 0.2, W(i,j-1,k) is the first point adjacent to the point W(i,j,k) in the negative direction of the y axis, W(i,j-2,k) is the second point adjacent to the point W(i,j,k) in the negative direction of the y axis, W(i,j-3,k) is the third point adjacent to the point W(i,j,k) in the negative direction of the y axis, W(i,j+1,k) is the first point adjacent to the point W(i,j,k) in the positive direction of the y axis, W(i,j+2,k) is the second point adjacent to the point W(i,j,k) in the positive direction of the y axis, and W(i,j+3,k) is the third point adjacent to the point W(i,j,k) in the positive direction of the y axis.

[0092] The initial structure information is obtained according to the first weighted difference operator and the second weighted difference operator, and the calculation formula is as follows:

[0093]

[0094] wherein S(i,j,k) is the initial structure information, C x (i,j,k) is the first weighted difference operator, C y (i,j,k) is the second weighted difference operator, and it can be seen from the above formula that the sum of the cube of the first weighted difference operator and the cube of the second weighted difference operator is the initial structure information.

[0095] Step S103: calculating the threshold value according to the initial structure information;

[0096] The threshold value is calculated according to the following formula:

[0097]

[0098] wherein thre is the threshold value, S(i,j,k) is the initial structure information, ε is the set control coefficient, n1 is the data number of the x axis of the three-dimensional seismic data, n2 is the data number of the y axis of the three-dimensional seismic data, and n3 is the data number of the z axis of the three-dimensional seismic data, it can be seen from the above formula that the data in the x direction, the y direction and the z direction of the initial structure information S(i,j,k) is summed and averaged, then the set control system is multiplied by the average value to obtain the final threshold value, the initial structure information is screened through the threshold value, the influence of noise or small-scale structure with less contribution to tomographic inversion is removed, only the structure information with contribution to tomographic velocity modeling is left, and the structure information extraction for the depth domain tomographic velocity modeling is realized.

[0099] Step S104: based on the threshold, normalizing the initial structure information to obtain the final structure information, the calculation formula is as follows:

[0100]

[0101] Wherein, S(i,j,k) is the initial structure information, thre is the threshold, is the final structure information.

[0102] To prove the correctness and effectiveness of the method proposed in this example, and show that the method has higher precision, the following is explained by actual data test.

[0103] As Figure 2 shown, this is the three-dimensional seismic data input in the embodiment, used to verify whether the structure extracted by the patent is accurate. As Figure 3 shown, by setting the control coefficient epsilon to adjust the threshold value in this embodiment to 200, the corresponding extracted structure information can be seen. When the threshold value is low, the extracted structure information is rich in details, and all the structure forms are extracted, but the too detailed structure information is not conducive to tomographic velocity modeling. As Figure 4 shown, by setting the control coefficient epsilon to adjust the threshold value in this embodiment to 400, the corresponding extracted structure information can be seen. With the increase of the threshold value, the details of the extracted structure information are reduced. As Figure 5 shown, by setting the control coefficient epsilon to adjust the threshold value in this embodiment to 600, the corresponding extracted structure information can be seen. When the threshold value continues to increase, the details of the extracted structure information are reduced, but the main structure is retained, which also shows that the method of this embodiment is flexible, and the degree of detail of the extracted structure information can be adjusted according to needs.

[0104] The structure information extraction method based on threshold constraint proposed in this embodiment adopts a double-azimuth six-point weighted difference operator to calculate the initial structure information, and then calculates the threshold based on the initial structure information. The threshold is calculated based on the initial structure information; based on the threshold, the initial structure information is normalized to obtain the final structure information. This embodiment filters the initial structure information through the threshold, eliminates the influence of noise or small-scale structures with small contribution to tomographic inversion, and only leaves the structure information that contributes to tomographic velocity modeling, realizing the extraction of structure information for depth domain tomographic velocity modeling.

[0105] Example three,

[0106] The embodiment proposes a structure information extraction device based on threshold constraint, as Figure 6The shown include: data acquisition module, difference operator calculation module, initial structure information calculation module, threshold calculation module and final structure information output module; the data acquisition module is connected with the difference operator calculation module, the difference operator calculation module is connected with the initial structure information calculation module, the initial structure information calculation module is connected with the threshold calculation module, and the threshold calculation module is connected with the final structure information output module;

[0107] The data acquisition module is used for acquiring three-dimensional seismic data.

[0108] The difference operator calculation module is used for taking a point of the three-dimensional seismic data, calculating a first weighted difference operator of the point in a first direction and a second weighted difference operator in a second direction.

[0109] The initial structure information calculation module is used for obtaining initial structure information according to the first weighted difference operator and the second weighted difference operator.

[0110] The threshold calculation module is used for calculating a threshold according to the initial structure information.

[0111] The final structure information output module is used for normalizing the initial structure information by the threshold to obtain final structure information.

[0112] Embodiment four:

[0113] 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 threshold-constrained structure information extraction method.

[0114] The electronic device can be a mobile phone, a computer or a tablet computer, and comprises a memory and a processor, and the memory stores a computer program which is executed by the processor to implement the threshold-constrained structure information extraction method as described in the embodiments. It can be understood that the electronic device can further comprise an input / output (I / O) interface and a communication component.

[0115] The processor is configured to execute all or part of the steps of the threshold-constrained structure information extraction method as described in the above 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.

[0116] 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, for executing the threshold constraint based construction information extraction method described in the above embodiments.

[0117] Embodiment Five:

[0118] The embodiment provides a computer readable storage medium storing executable instructions, which, when executed, cause a processor to perform the threshold constraint based construction information extraction method.

[0119] If implemented in the form of a software functional unit and sold or used as an independent product, the software functional unit can be stored in a computer readable storage medium.

[0120] Based on such understanding, the technical solution of the present application, essentially or the part that contributes to the prior art, or part of the technical solution can be embodied in the form of a software product, which is stored in a storage medium, includes a plurality of instructions to make a computer device (which can be a personal computer, a server, or a network device, etc.) execute all or part of the steps of the threshold constraint based construction information extraction method described in the embodiments of the present application.

[0121] The aforementioned storage medium includes a flash memory, a hard disk, a multimedia card, a card-type memory (for example, 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 store, and various media that can store program check codes, on which a computer program is stored, and the computer program, when executed by a processor, can implement each step of the threshold constraint based construction information extraction method described above.

[0122] 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 referred to each other, and each embodiment focuses on the difference from other embodiments.

[0123] The scope 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 method for extracting structural information based on threshold constraints, characterized in that: include: Acquire 3D seismic data; Taking a point of the three-dimensional seismic data, calculating a first weighted difference operator of the point in a first direction and a second weighted difference operator in a second direction; Obtaining initial construction information according to the first weighted difference operator and the second weighted difference operator; Calculate a threshold value according to the initial construction information; Based on the threshold, the initial construction information is normalized to obtain final construction information.

2. The method for extracting structural information based on threshold constraints according to claim 1, characterized in that: 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. A point of the three-dimensional seismic data is represented as , then the first weighted difference operator in the first direction is the difference between the point The result of weighted difference of the adjacent N points; the second weighted difference operator in the second direction is the weighted difference operator in the y-axis direction with the point The result of weighted difference of N adjacent points, where N is an integer multiple of 2.

3. The method for extracting structural information based on threshold constraints according to claim 2, characterized in that: The point in the x-axis direction The N adjacent points are: … 、 、 … as well as ; The point in the y-axis direction The N adjacent points are: … 、 、 … as well as .

4. The method for extracting structural information based on threshold constraints according to claim 1, characterized in that: The first weighted difference operator is calculated as follows: in, is the first weighted difference operator, Set the weight for the first one, For the Set weights, For the Set weights, The point in the negative direction of the x-axis The first adjacent point, The point in the negative direction of the x-axis The neighboring Points, The point in the negative direction of the x-axis The neighboring Points, is the point in the positive direction of the x-axis The first adjacent point, is the point in the positive direction of the x-axis The neighboring Points, is the point in the positive direction of the x-axis The neighboring points.

5. The method for extracting structural information based on threshold constraints according to claim 1, characterized in that: The second weighted difference operator is calculated as follows: in, is the second weighted difference operator, Set the weight for the first one, For the Set weights, For the Set weights, In the negative direction of the y-axis The first adjacent point, In the negative direction of the y-axis The neighboring Points, In the negative direction of the y-axis The neighboring Points, In the positive direction of the y-axis The first adjacent point, In the positive direction of the y-axis The neighboring Points, In the positive direction of the y-axis The neighboring points.

6. The method for extracting structural information based on threshold constraints according to claim 1, characterized in that: The initial construction information is obtained according to the first weighted difference operator and the second weighted difference operator, and the calculation formula is as follows: in, is the initial construction information, is the first weighted difference operator, is the second weighted difference operator.

7. The method for extracting structural information based on threshold constraints according to claim 1, characterized in that: The threshold is calculated as follows: in, is the threshold, is the initial construction information, is the set control coefficient, is the number of data on the x-axis of the 3D seismic data, The number of data on the y-axis of 3D seismic data, is the number of data on the z-axis of the 3D seismic data.

8. The method for extracting structural information based on threshold constraints according to claim 1, characterized in that: The initial structure information is normalized by a threshold to obtain the final structure information, which is calculated as follows: in, is the initial construction information, is the threshold, The final construction information.

9. A structural information extraction device based on threshold constraints, characterized in that: include: A data acquisition module, used for acquiring three-dimensional seismic data; A differential operator calculation module, configured to take a point of the three-dimensional seismic data and calculate a first weighted differential operator of the point in a first direction and a second weighted differential operator in a second direction; an initial construction information calculation module, configured to obtain initial construction information according to the first weighted difference operator and the second weighted difference operator; A threshold calculation module, configured to calculate a threshold based on the initial configuration information; The final structure information output module is used to normalize the initial structure information through a threshold to obtain the final structure information.

10. 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 method for extracting structural information based on threshold constraints according to any one of claims 1 to 8 is executed.

11. A computer-readable storage medium, characterized in that The device stores executable instructions, which, when executed, enable a processor to execute the construction information extraction method based on threshold constraints according to any one of claims 1 to 8.

Images