Construction information extraction method and device based on threshold constraint

Abstract

The invention provides a construction information extraction method and device based on threshold constraint. The method comprises the following steps: acquiring 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 construction information according to the first weighted difference operator and the second weighted difference operator; calculating a threshold value according to the initial construction information; and based on a threshold value, performing normalization processing on the initial construction information to obtain final construction information. The device comprises a data acquisition module, a difference operator calculation module, an initial construction information calculation module, a threshold calculation module and a final construction information output module. According to the method, the influence of noise or small-scale structures with small contribution to tomography inversion and the like is eliminated, only the structure information contributing to tomography velocity modeling is left, and the structure information extraction oriented to depth domain tomography velocity modeling is realized.

Description

Technical Field

[0001] The present application belongs to the field of geophysical exploration technology, and specifically relates to a method and device for extracting structural information based on threshold constraints. Background Art

[0002] At present, conventional deep-domain tomographic velocity modeling methods usually do not consider structural information, the inversion velocity model lacks geological significance, and the inversion accuracy is low. To address this issue, how to add geological structural information to constrain the tomographic velocity modeling process is the key to solving this problem. In this technical process, how to extract structural information efficiently and with high quality is the first problem that needs to be solved. The current method of extracting structure is to obtain structural information by means of slope or local energy superposition. This method is usually cumbersome to calculate, and the calculation results are not stable enough. The extracted structural information cannot be quantitatively controlled by certain parameters, and additional filtering and other technologies are required for modification and modification. In actual application, considering issues such as processing cycle and complexity, it is less used. Summary of the invention

[0003] Based on the above technical problems, the present application proposes a method and device for extracting structural information based on threshold constraints.

[0004] In a first aspect, the present application proposes a method for extracting structural information based on threshold constraints, comprising:

[0005] Acquire 3D 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 construction information according to the first weighted difference operator and the second weighted difference operator;

[0008] Calculate a threshold value according to the initial construction information;

[0009] Based on the threshold, the initial construction information is normalized to obtain final construction information.

[0010] 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 connecting 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 by W(i,j,k), then the first weighted difference operator in the first direction is the result of the weighted difference of the 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 the result of the weighted difference of the N points adjacent to the point W(i,j,k) in the y-axis direction, where N is an integer multiple of 2.

[0011] The N points adjacent to the 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 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 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] 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 three-dimensional seismic data, n2 is the number of data on the y-axis of the three-dimensional seismic data, and n3 is the number of data on the z-axis of the three-dimensional seismic data.

[0027] The initial structure information is normalized by a threshold to obtain the final structure information, and the calculation formula is 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, used for taking a point of the three-dimensional seismic data, and calculating 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 construction information calculation module, used to obtain initial construction information according to the first weighted difference operator and the second weighted difference operator;

[0034] A threshold calculation module, used to calculate a threshold according to the initial construction information;

[0035] The final structure information output module is used to normalize the initial structure information through a threshold to obtain the final structure information.

[0036] 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 threshold constraint-based construction information extraction method is executed.

[0037] 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 threshold constraint-based construction information extraction method.

[0038] Beneficial effects:

[0039] The present application proposes a method and device for extracting structural information based on threshold constraints. The initial structural information is calculated using a weighted difference operator, and then a threshold is calculated based on the initial structural information. The final structural information is obtained through threshold constraints. The method of the present application eliminates the influence of noise or small-scale structures that contribute little to tomographic inversion, leaving only structural information that contributes to tomographic velocity modeling, thereby realizing structural information extraction for deep-domain tomographic velocity modeling. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] Figure 1 This is a flowchart of a method for extracting structural information based on threshold constraints 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 This is a schematic diagram of the structure information extracted when the threshold value is 200 according to the embodiment of the present application;

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

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

[0045] Figure 6 This is a principle block diagram of the device for extracting structural information based on threshold constraints of this embodiment. DETAILED DESCRIPTION

[0046] The present disclosure is further described below in conjunction with the embodiments shown in the accompanying drawings.

[0047] At present, conventional depth domain tomography velocity modeling methods usually do not consider structural information, the inversion velocity model lacks geological significance, and the inversion accuracy is low. The current method of extracting structural information is to obtain structural information by slope or local energy superposition. This method is usually cumbersome to calculate, and the calculation results are not stable enough. The extracted structural information cannot be quantitatively controlled by a certain parameter, and additional filtering and other technologies are required for modification and modification. In actual application, considering issues such as processing cycle and complexity, it is rarely used.

[0048] The present application proposes a structural information extraction method based on threshold constraints, which adopts a two-orientation N-point weighted difference operator to calculate the initial structural information. Based on the initial structural information, a threshold is set for screening to eliminate the influence of noise or small-scale structures that have a small 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.

[0049] Embodiment 1

[0050] This embodiment proposes a construction information extraction method based on threshold constraints, such as Figure 1 As shown, including:

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

[0052] In this embodiment, the three-dimensional seismic data is three-dimensional seismic data obtained under a three-dimensional coordinate system consisting of an x-axis as the horizontal axis, a y-axis as the vertical axis, and a z-axis as the vertical axis, wherein the x-axis direction is data in the direction of the main seismic survey line, the y-axis direction is data in the direction of the seismic connecting survey line, and the z-axis direction is data in the direction of the seismic depth.

[0053] Step S2: 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 in a second direction;

[0054] In this embodiment, a point of the three-dimensional seismic data is represented by W(i,j,k), then the first weighted difference operator in the first direction is the result of the weighted difference 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 the result of the weighted difference of N points adjacent to the point W(i,j,k) in the y-axis direction, where 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 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 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: Calculate a threshold value according to 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 three-dimensional seismic data, n2 is the number of data on the y-axis of the three-dimensional seismic data, and n3 is the number of data on the z-axis of the three-dimensional seismic data. It can be seen from the above formula that the data of the initial structural information S(i,j,k) in the x-direction, y-direction and z-direction are added 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 screened by the threshold to eliminate the influence of noise or small-scale structures that contribute little to tomographic inversion, leaving only the structural information that contributes to tomographic velocity modeling, thereby realizing the extraction of structural information for deep domain tomographic velocity modeling.

[0073] Step S5: Based on the threshold, the initial construction information is normalized to obtain the final construction information, and the calculation formula is as follows:

[0074]

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

[0076] This embodiment proposes a structural information extraction method based on threshold constraints, taking a point of three-dimensional seismic data, calculating the first weighted difference operator of the point in the first direction and the second weighted difference operator in the 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 normalizing the initial structural information based on the threshold to obtain the final structural information. This embodiment uses the threshold to filter the initial structural information, eliminate the influence of noise or small-scale structures that contribute little to tomographic inversion, and only leave structural information that contributes to tomographic velocity modeling, thereby realizing structural information extraction for tomographic velocity modeling in the depth domain.

[0077] Embodiment 2

[0078] This embodiment is a specific example of the first embodiment. The initial construction information is calculated by using a two-direction six-point weighted difference operator. Then, a threshold is calculated based on the initial construction information. The initial construction information is screened by setting the threshold. The detailed steps are as follows:

[0079] Step S100: Acquire three-dimensional seismic data;

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

[0081] Step S101: 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 in a second direction;

[0082] In this embodiment, a point of the three-dimensional seismic data is represented as W(i,j,k), then the first weighted difference operator in the first direction is the result of the weighted difference of the six points adjacent to the point W(i,j,k) in the x-axis direction; the second weighted difference operator in the second direction is the result of the weighted difference of the 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 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 using six points adjacent to point W(i,j,k) in the x-axis direction. The calculation formula is as follows:

[0087]

[0088] Among them, C x (i,j,k) is the first weighted difference operator, q1 is the first set weight, and the first set weight in this embodiment is 0.6, q2 is the second set weight, and the second set weight in this embodiment is 0.2, q3 is the third set weight, and the third set weight in this embodiment is 0.3, W(i-1,j,k) is the first point adjacent to point W(i,j,k) in the negative direction of the x-axis, W(i-2,j,k) is the second point adjacent to point W(i,j,k) in the negative direction of the x-axis, W(i-3,j,k) is the third point adjacent to point W(i,j,k) in the negative direction of the x-axis, W(i+1,j,k) is the first point adjacent to point W(i,j,k) in the positive direction of the x-axis, W(i+2,j,k) is the second point adjacent to 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 point W(i,j,k) in the positive direction of the x-axis.

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

[0090]

[0091] Among them, C y(i,j,k) is the second weighted difference operator, q1 is the first set weight, which is 0.6 in this embodiment, q2 is the second set weight, which is 0.2 in this embodiment, q3 is the third set weight, which is 0.2 in this embodiment, W(i,j-1,k) is the first point adjacent to point W(i,j,k) in the negative direction of the y-axis, W(i,j-2,k) is the first point adjacent to point W(i,j,k) in the negative direction of the y-axis, k), W(i,j-3,k) is the third point adjacent to point W(i,j,k) in the negative direction of the y-axis, W(i,j+1,k) is the first point adjacent to point W(i,j,k) in the positive direction of the y-axis, W(i,j+2,k) is the second point adjacent to 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 point W(i,j,k) in the positive direction of the y-axis.

[0092] 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:

[0093]

[0094] 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.

[0095] Step S103: Calculate a threshold value according to the initial construction information;

[0096] The threshold is calculated as follows:

[0097]

[0098] 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 three-dimensional seismic data, n2 is the number of data on the y-axis of the three-dimensional seismic data, and n3 is the number of data on the z-axis of the three-dimensional seismic data. It can be seen from the above formula that the data of the initial structural information S(i,j,k) in the x-direction, y-direction and z-direction are added 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 screened by the threshold to eliminate the influence of noise or small-scale structures that contribute little to tomographic inversion, leaving only the structural information that contributes to tomographic velocity modeling, thereby realizing the extraction of structural information for deep domain tomographic velocity modeling.

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

[0100]

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

[0102] In order to prove the correctness and effectiveness of the method proposed in this example and to show that the method has higher accuracy, an actual data test is used to illustrate it.

[0103] like Figure 2 As shown in FIG. 1 , this is the 3D seismic data input in the embodiment, which is used to verify whether the structure extracted by the patent of the present invention is accurate. Figure 3 As shown in FIG. 1 , by adjusting the threshold value of 200 in this embodiment by setting the control coefficient ε, it can be seen that when the threshold value is low, the extracted structural information is rich in details and all structural forms are extracted, but too detailed structural information is not conducive to tomographic velocity modeling. Figure 4 As shown in FIG. 1 , by adjusting the structure information extracted when the threshold is 400 in this embodiment through the set control coefficient ε, it can be seen that as the threshold is increased, the details of the extracted structure information are reduced. Figure 5 As shown, by setting the control coefficient ε to adjust the corresponding extracted structural information when the threshold value is 600 in the embodiment, it can be seen that when the threshold value continues to increase, the details of the extracted structural information are reduced, but the main structure is retained, which also shows that the method of this embodiment is more flexible and the degree of detail of the extracted structural information can be adjusted as needed.

[0104] This embodiment proposes a structural information extraction method based on threshold constraints, which uses a two-directional six-point weighted difference operator to calculate the initial structural information, and then calculates the threshold based on the initial structural information, and calculates the threshold based on the initial structural information; based on the threshold, the initial structural information is normalized to obtain the final structural information. This embodiment uses the threshold to filter the initial structural information, eliminate the influence of noise or small-scale structures that contribute little to tomographic inversion, and only retains the structural information that contributes to tomographic velocity modeling, thereby realizing structural information extraction for deep-domain tomographic velocity modeling.

[0105] Embodiment 3

[0106] This embodiment proposes a structure information extraction device based on threshold constraints, such as Figure 6The system includes: a data acquisition module, a differential operator calculation module, an initial structure information calculation module, a threshold calculation module and a final structure information output module; the data acquisition module is connected to the differential operator calculation module, the differential operator calculation module is connected to the initial structure information calculation module, the initial structure information calculation module is connected to the threshold calculation module, and the threshold calculation module is connected to the final structure information output module;

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

[0108] A differential operator calculation module, used for taking a point of the three-dimensional seismic data, and calculating a first weighted differential operator of the point in a first direction and a second weighted differential operator in a second direction;

[0109] An initial construction information calculation module, used to obtain initial construction information according to the first weighted difference operator and the second weighted difference operator;

[0110] A threshold calculation module, used to calculate a threshold according to the initial construction information;

[0111] The final structure information output module is used to normalize the initial structure information through a threshold to obtain the final structure information.

[0112] Embodiment 4:

[0113] This embodiment provides an electronic device, 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 method for extracting construction information based on threshold constraints is executed.

[0114] 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 on the memory, and when the computer program is executed by the processor, the construction information extraction method based on threshold constraints 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.

[0115] The processor is used to execute all or part of the steps in the construction information extraction method based on threshold constraints 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.

[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 components, and is used to execute the threshold constraint-based construction information extraction method described in the above embodiment.

[0117] Embodiment five:

[0118] This embodiment provides a computer-readable storage medium storing executable instructions, which, when executed, enable a processor to execute 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, it can be stored in a computer-readable storage medium.

[0120] Based on this understanding, the technical solution of the present application can essentially or partly be embodied in the form of a software product that contributes to the prior art. 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 threshold constraint-based construction information extraction method described in each embodiment of the present application.

[0121] 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 threshold constraint-based construction information extraction method can be implemented.

[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 referenced to each other, and each embodiment focuses on the differences from other embodiments.

[0123] 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 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 connecting 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 by W(i,j,k), then the first weighted difference operator in the first direction is the result of the weighted difference of the 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 the result of the weighted difference of the N points adjacent to the point W(i,j,k) in the y-axis direction, 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 N points adjacent to point W(i,j,k) in the x-axis direction are: as well as The N points adjacent to point W(i,j,k) in the y-axis direction are: W(i,j-1,k)… 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: Among them, C x (i,j,k) is the first weighted difference operator, q1 is the first set weight, q2 is the first 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.

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: Among them, C y (i,j,k) is the second weighted difference operator, q1 is the first set weight, q2 is the 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.

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: 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.

7. The method for extracting structural information based on threshold constraints according to claim 1, characterized in that: The threshold is calculated as follows: 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 three-dimensional seismic data, n2 is the number of data on the y-axis of the three-dimensional seismic data, and n3 is the number of data on the z-axis of the three-dimensional 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, and the calculation formula is as follows: Among them, S(i,j,k) is the initial construction information, thre 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, used for taking a point of the three-dimensional seismic data, and calculating 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, used to obtain initial construction information according to the first weighted difference operator and the second weighted difference operator; A threshold calculation module, used to calculate a threshold according to the initial construction 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: It stores executable instructions, which, when executed, enable a processor to execute the construction information extraction method based on threshold constraints as described in any one of claims 1 to 8.

Images