Method, device, computer device and storage medium for determining fault throw

Abstract

The application discloses a method and device for determining fault throw, a computer device and a storage medium, and belongs to the technical field of seismic data interpretation. The method comprises the following steps: dividing a plurality of feature points into a plurality of sets based on the horizontal coordinates of the plurality of feature points; determining a fault zone core feature curve, an up-dip feature curve and a down-dip feature curve for each set; determining a plurality of first reference points on the up-dip feature curve, determining a first normal line of the up-dip feature curve at each first reference point, determining a first intersection point of the first normal line and the fault zone core feature curve, determining a plurality of second reference points on the down-dip feature curve, determining a second normal line of the down-dip feature curve at each second reference point, and determining a second intersection point of the second normal line and the fault zone core feature curve; determining a reference point pair; determining a horizon value of the reference point pair; and determining a fault throw corresponding to the reference point pair. The application can improve the efficiency of determining fault throw.

Description

Technical Field

[0001] This application relates to the field of seismic data interpretation technology, and in particular to a method, apparatus, computer equipment, and storage medium for determining fault displacement. Background Technology

[0002] Faults are an important manifestation of stress release during crustal movement. The three elements describing a fault include the fault plane, fault line, and fault displacement. Among these, fault displacement reflects the intensity of crustal movement and is one of the most important evaluation parameters for describing faults. Accurate measurement of fault displacement is of great significance in geological research and mineral exploration.

[0003] Currently, fault displacement is mostly measured manually. Technicians observe the images corresponding to the fault zone planar data and, through experience, determine the corresponding feature points in the uplift and downlift blocks, as well as the core feature points of the fault zone, from a large number of feature points. Then, the fault displacement is calculated based on the stratigraphic values ​​of these points. This method of measuring fault displacement is inefficient, especially when the fault zone planar data is very complex, and it consumes a lot of time. Summary of the Invention

[0004] This application provides a method, apparatus, computer device, and storage medium for determining fault displacement, which can solve related problems. The technical solution is as follows:

[0005] In a first aspect, a method for determining fault displacement is provided, the method comprising:

[0006] Obtain the fault zone plane data of the target fault, wherein the fault zone plane data includes the horizontal coordinates and stratigraphic values ​​of multiple feature points of the target fault, and some of the feature points are core feature points of the fault zone.

[0007] Based on the horizontal coordinates of the multiple feature points, the multiple feature points are divided into multiple sets, wherein for any feature point, there is at least one other feature point in the set to which the feature point belongs that the distance between the feature point and the feature point is less than a distance threshold.

[0008] For each set, curve fitting is performed based on the core feature points of the fault zone in the set to obtain the core feature curve of the fault zone. Curve fitting is performed based on feature points other than the core feature points of the fault zone in the set to obtain the first feature curve. Based on the intersection of the first feature curve and the core feature curve of the fault zone, the first feature curve is divided into a second feature curve and a third feature curve. Based on the stratigraphic values ​​of each feature point corresponding to the second feature curve and the third feature curve, the uplifted block feature curve and the downlifted block feature curve are determined in the second feature curve and the third feature curve.

[0009] For each rising disk characteristic curve, multiple first reference points are determined on the rising disk characteristic curve at a first preset interval, and a first normal line of the rising disk characteristic curve at each first reference point is determined, and a first intersection point of each first normal line with the core characteristic curve of the fracture zone is determined.

[0010] For each descending disk characteristic curve, multiple second reference points are determined on the descending disk characteristic curve at the first preset interval, and a second normal line of the descending disk characteristic curve at each second reference point is determined, and a second intersection point of each second normal line with the core characteristic curve of the fracture zone is determined.

[0011] For each first intersection point, determine the second intersection point that is closest to the first intersection point and less than the second preset distance, and determine the first reference point corresponding to the first intersection point and the second reference point corresponding to the second intersection point as a reference point pair;

[0012] Based on the layer values ​​of the multiple feature points, determine the layer value corresponding to each first reference point, each second reference point, each first intersection point, and each second intersection point;

[0013] For each reference point pair, a first fault displacement is determined based on the stratigraphic value of the first reference point in the reference point pair and the stratigraphic value of the corresponding first intersection point. A second fault displacement is determined based on the stratigraphic value of the second reference point in the reference point pair and the stratigraphic value of the corresponding second intersection point. The sum of the first fault displacement and the second fault displacement is determined to obtain the fault displacement corresponding to the reference point pair. The fault displacements corresponding to all reference point pairs are determined as the fault displacements of the target fault.

[0014] In one possible implementation, prior to acquiring the fracture zone plane data of the target fault, the method further includes:

[0015] Acquire seismic data volume for a target area, wherein the seismic data volume includes the horizontal coordinates, layer values, fault zone core probability, and attribute values ​​of multiple fault attributes for each spatial point in the target area, wherein the fault zone core probability is used to represent the probability that the spatial point is located at the core of the fault zone;

[0016] For each fault attribute, based on the horizontal coordinates, layer values, and attribute values ​​of each spatial point in the seismic data volume, the fault corresponding to each spatial point is determined, and multiple target spatial points corresponding to the first surface of the target fault are determined. The horizontal coordinates and fault zone core probability of the target spatial points in the seismic data volume are determined as the initial plane data of the fault zone of the target fault, wherein the first surface is the upper or lower surface of the target fault.

[0017] In the initial plane data of the fault zone of the target fault corresponding to each fault attribute, the initial plane data of the first fault zone is selected;

[0018] In the initial planar data of the first fault zone, the target horizontal coordinates of the corresponding fault zone core probability within a preset interval are determined. The target spatial points with the target horizontal coordinates among multiple target spatial points are determined as feature points of the target fault. The horizontal coordinates and stratigraphic values ​​of the feature points of the target fault are determined as the fault zone planar data of the target fault.

[0019] In one possible implementation, the plurality of fault attributes includes at least two of the following: coherence attributes, curvature attributes, ant body attributes, energy attributes, and artificial intelligence fault prediction attributes.

[0020] In one possible implementation, selecting the first initial plane data of the fault zone from the initial plane data of the target fault corresponding to each attribute includes:

[0021] Obtain the initial planar data of the fracture zone of the target fault corresponding to each fault attribute;

[0022] Determine the data continuity value and signal-to-noise ratio corresponding to the initial planar data of each fault zone;

[0023] For each initial plane data of a fault zone, a weighted average of the data continuity value and the signal-to-noise ratio corresponding to the initial plane data of the fault zone is determined based on the first weighting coefficient corresponding to the data continuity value and the second weighting coefficient corresponding to the signal-to-noise ratio, and is used as the selection score corresponding to the initial plane data of the fault zone.

[0024] Among all the initial plane data of the fault zones, the first initial plane data of the fault zone with the highest selection score is selected.

[0025] In one possible implementation, the preset interval includes a fault zone core interval and a fault boundary interval, wherein the value range of the fault zone core interval is [0.9, 1], the value range of the fault boundary interval is [0.2, 0.3], and the target horizontal coordinate includes a first horizontal coordinate corresponding to the fault zone core probability within the fault zone core interval and a second horizontal coordinate corresponding to the fault zone core probability within the fault boundary interval.

[0026] The method further includes:

[0027] The target spatial point with the first horizontal coordinate among multiple target spatial points is determined as the core feature point of the fracture zone of the target fault.

[0028] In one possible implementation, determining the layer value corresponding to each first reference point, each second reference point, each first intersection point, and each second intersection point based on the layer values ​​of the plurality of feature points includes:

[0029] For any first reference point, if the first reference point coincides with the first feature point among the feature points corresponding to the first ascending disk feature curve where the first reference point is located, then the layer value of the first feature point is determined as the layer value corresponding to the first reference point. If the first reference point does not coincide with any of the feature points corresponding to the first ascending disk feature curve, then the layer value corresponding to the first reference point is determined based on the layer values ​​of the two feature points closest to the first reference point among the feature points corresponding to the first ascending disk feature curve. Here, each feature point corresponding to the first ascending disk feature curve is a feature point used to fit the first ascending disk feature curve.

[0030] For any second reference point, if the second reference point coincides with the second feature point among the feature points corresponding to the first descending disk feature curve where the second reference point is located, then the layer value of the second feature point is determined as the layer value corresponding to the second reference point. If the second reference point does not coincide with any of the feature points corresponding to the first descending disk feature curve, then the layer value corresponding to the second reference point is determined based on the layer values ​​of the two feature points closest to the second reference point among the feature points corresponding to the first descending disk feature curve. Here, each feature point corresponding to the first descending disk feature curve is a feature point used to fit the first descending disk feature curve.

[0031] For any first intersection point, if the first intersection point coincides with the third feature point among the feature points corresponding to the core feature curve of the first fault zone where the first intersection point is located, then the stratigraphic value of the third feature point is determined as the stratigraphic value corresponding to the first intersection point. If the first intersection point does not coincide with any of the feature points corresponding to the core feature curve of the first fault zone, then the stratigraphic value corresponding to the first intersection point is determined based on the stratigraphic values ​​of the two feature points closest to the first intersection point among the feature points corresponding to the core feature curve of the first fault zone. Here, each feature point corresponding to the core feature curve of the first fault zone is a feature point used to fit the core feature curve of the first fault zone.

[0032] For any second intersection point, if the second intersection point coincides with the fourth feature point among the feature points corresponding to the core feature curve of the first fault zone where the second intersection point is located, then the stratigraphic value of the fourth feature point is determined as the stratigraphic value corresponding to the second intersection point. If the second intersection point does not coincide with any of the feature points corresponding to the core feature curve of the first fault zone, then the stratigraphic value corresponding to the second intersection point is determined based on the stratigraphic values ​​of the two feature points closest to the second intersection point among the feature points corresponding to the core feature curve of the first fault zone.

[0033] In one possible implementation, determining the stratigraphic value corresponding to the first reference point based on the stratigraphic values ​​of the two feature points closest to the first reference point among all feature points corresponding to the first ascending disk feature curve includes:

[0034] Determine the two feature points that are closest to the first reference point among the feature points corresponding to the first rising plate feature curve, wherein the two feature points include the fifth feature point and the sixth feature point;

[0035] On the first rising plate characteristic curve, determine the first curve point closest to the fifth characteristic point and the second curve point closest to the sixth characteristic point;

[0036] Based on the first arc length between the first curve point and the first reference point, the second arc length between the second curve point and the first reference point, and the layer value of the fifth feature point and the layer value of the sixth feature point, the layer value corresponding to the first reference point is determined.

[0037] The step of determining the stratigraphic value corresponding to the second reference point based on the stratigraphic values ​​of the two feature points closest to the second reference point among the feature points corresponding to the first descending disk feature curve includes:

[0038] Determine the two feature points that are closest to the second reference point among the feature points corresponding to the first descending disk feature curve, wherein the two feature points include the seventh feature point and the eighth feature point;

[0039] On the first descending disk characteristic curve, determine the third curve point that is closest to the seventh characteristic point and the fourth curve point that is closest to the eighth characteristic point;

[0040] Based on the third arc length between the third curve point and the second reference point, the fourth arc length between the fourth curve point and the second reference point, and the layer value of the seventh feature point and the layer value of the eighth feature point, the layer value corresponding to the second reference point is determined.

[0041] The step of determining the stratigraphic value corresponding to the first intersection point based on the stratigraphic values ​​of the two feature points closest to the first intersection point among the feature points corresponding to the core feature curve of the first fault zone includes:

[0042] Identify the two feature points that are closest to the first intersection point among the feature points corresponding to the core feature curve of the first fault zone, wherein the two feature points include the ninth feature point and the tenth feature point;

[0043] On the core characteristic curve of the first fault zone, determine the fifth curve point that is closest to the ninth characteristic point and the sixth curve point that is closest to the tenth characteristic point;

[0044] Based on the fifth arc length between the fifth curve point and the first intersection point, the sixth arc length between the sixth curve point and the first intersection point, the layer value of the ninth feature point and the layer value of the tenth feature point, the layer value corresponding to the first intersection point is determined.

[0045] The step of determining the stratigraphic value corresponding to the second intersection point based on the stratigraphic values ​​of the two feature points closest to the second intersection point among the feature points corresponding to the core feature curve of the first fault zone includes:

[0046] Identify the two feature points that are closest to the second intersection point among the feature points corresponding to the core feature curve of the first fault zone, wherein the two feature points include the eleventh feature point and the twelfth feature point;

[0047] On the core characteristic curve of the first fault zone, determine the seventh curve point that is closest to the eleventh characteristic point and the eighth curve point that is closest to the twelfth characteristic point;

[0048] Based on the seventh arc length between the seventh curve point and the second intersection point, the eighth arc length between the eighth curve point and the second intersection point, and the layer value of the eleventh feature point and the layer value of the twelfth feature point, the layer value corresponding to the second intersection point is determined.

[0049] Secondly, an apparatus for determining fault displacement is provided, the apparatus comprising:

[0050] The acquisition module is used to acquire the fracture zone plane data of the target fault, wherein the fracture zone plane data includes the horizontal coordinates and stratigraphic values ​​of multiple feature points of the target fault, and some of the feature points are core feature points of the fracture zone.

[0051] The partitioning module is used to divide the multiple feature points into multiple sets based on the horizontal coordinates of the multiple feature points, wherein, for any feature point, there is at least one other feature point in the set to which the feature point belongs, and the distance between the feature point and the feature point is less than a distance threshold.

[0052] The fitting module is used to perform curve fitting based on the core feature points of the fault zone in each set to obtain the core feature curve of the fault zone, and to perform curve fitting based on feature points other than the core feature points of the fault zone in the set to obtain the first feature curve. Based on the intersection of the first feature curve and the core feature curve of the fault zone, the first feature curve is divided into a second feature curve and a third feature curve. Based on the stratigraphic values ​​of each feature point corresponding to the second feature curve and the third feature curve, the rising block feature curve and the falling block feature curve are determined in the second feature curve and the third feature curve.

[0053] The first determining module is used to determine multiple first reference points on each rising disk feature curve at a first preset interval, determine the first normal line of the rising disk feature curve at each first reference point, and determine the first intersection point of each first normal line with the core feature curve of the fracture zone.

[0054] For each descending disk characteristic curve, multiple second reference points are determined on the descending disk characteristic curve at the first preset interval, and a second normal line of the descending disk characteristic curve at each second reference point is determined, and a second intersection point of each second normal line with the core characteristic curve of the fracture zone is determined.

[0055] For each first intersection point, determine the second intersection point that is closest to the first intersection point and less than the second preset distance, and determine the first reference point corresponding to the first intersection point and the second reference point corresponding to the second intersection point as a reference point pair;

[0056] Based on the layer values ​​of the multiple feature points, determine the layer value corresponding to each first reference point, each second reference point, each first intersection point, and each second intersection point;

[0057] The second determining module is used to determine a first fault displacement for each reference point pair based on the stratigraphic value of the first reference point in the reference point pair and the stratigraphic value of the corresponding first intersection point; determine a second fault displacement based on the stratigraphic value of the second reference point in the reference point pair and the stratigraphic value of the corresponding second intersection point; determine the sum of the first fault displacement and the second fault displacement to obtain the fault displacement corresponding to the reference point pair; and determine the fault displacement corresponding to all reference point pairs as the fault displacement of the target fault.

[0058] In one possible implementation, the acquisition module is further configured to:

[0059] Acquire seismic data volume for a target area, wherein the seismic data volume includes the horizontal coordinates, layer values, fault zone core probability, and attribute values ​​of multiple fault attributes for each spatial point in the target area, wherein the fault zone core probability is used to represent the probability that the spatial point is located at the core of the fault zone;

[0060] The first determining module is used for:

[0061] For each fault attribute, based on the horizontal coordinates, layer values, and attribute values ​​of each spatial point in the seismic data volume, the fault corresponding to each spatial point is determined, and multiple target spatial points corresponding to the first surface of the target fault are determined. The horizontal coordinates and fault zone core probability of the target spatial points in the seismic data volume are determined as the initial plane data of the fault zone of the target fault, wherein the first surface is the upper or lower surface of the target fault.

[0062] In the initial plane data of the fault zone of the target fault corresponding to each fault attribute, the initial plane data of the first fault zone is selected;

[0063] In the initial planar data of the first fault zone, the target horizontal coordinates of the corresponding fault zone core probability within a preset interval are determined. The target spatial points with the target horizontal coordinates among multiple target spatial points are determined as feature points of the target fault. The horizontal coordinates and stratigraphic values ​​of the feature points of the target fault are determined as the fault zone planar data of the target fault.

[0064] In one possible implementation, the plurality of fault attributes includes at least two of the following: coherence attributes, curvature attributes, ant body attributes, energy attributes, and artificial intelligence fault prediction attributes.

[0065] In one possible implementation, the first determining module is configured to:

[0066] Obtain the initial planar data of the fracture zone of the target fault corresponding to each fault attribute;

[0067] Determine the data continuity value and signal-to-noise ratio corresponding to the initial planar data of each fault zone;

[0068] For each initial plane data of a fault zone, a weighted average of the data continuity value and the signal-to-noise ratio corresponding to the initial plane data of the fault zone is determined based on the first weighting coefficient corresponding to the data continuity value and the second weighting coefficient corresponding to the signal-to-noise ratio, and is used as the selection score corresponding to the initial plane data of the fault zone.

[0069] Among all the initial plane data of the fault zones, the first initial plane data of the fault zone with the highest selection score is selected.

[0070] In one possible implementation, the preset interval includes a fault zone core interval and a fault boundary interval, wherein the value range of the fault zone core interval is [0.9, 1], the value range of the fault boundary interval is [0.2, 0.3], and the target horizontal coordinate includes a first horizontal coordinate corresponding to the fault zone core probability within the fault zone core interval and a second horizontal coordinate corresponding to the fault zone core probability within the fault boundary interval.

[0071] The first determining module is further configured to:

[0072] The target spatial point with the first horizontal coordinate among multiple target spatial points is determined as the core feature point of the fracture zone of the target fault.

[0073] In one possible implementation, the second determining module is configured to:

[0074] For any first reference point, if the first reference point coincides with the first feature point among the feature points corresponding to the first ascending disk feature curve where the first reference point is located, then the layer value of the first feature point is determined as the layer value corresponding to the first reference point. If the first reference point does not coincide with any of the feature points corresponding to the first ascending disk feature curve, then the layer value corresponding to the first reference point is determined based on the layer values ​​of the two feature points closest to the first reference point among the feature points corresponding to the first ascending disk feature curve. Here, each feature point corresponding to the first ascending disk feature curve is a feature point used to fit the first ascending disk feature curve.

[0075] For any second reference point, if the second reference point coincides with the second feature point among the feature points corresponding to the first descending disk feature curve where the second reference point is located, then the layer value of the second feature point is determined as the layer value corresponding to the second reference point. If the second reference point does not coincide with any of the feature points corresponding to the first descending disk feature curve, then the layer value corresponding to the second reference point is determined based on the layer values ​​of the two feature points closest to the second reference point among the feature points corresponding to the first descending disk feature curve. Here, each feature point corresponding to the first descending disk feature curve is a feature point used to fit the first descending disk feature curve.

[0076] For any first intersection point, if the first intersection point coincides with the third feature point among the feature points corresponding to the core feature curve of the first fault zone where the first intersection point is located, then the stratigraphic value of the third feature point is determined as the stratigraphic value corresponding to the first intersection point. If the first intersection point does not coincide with any of the feature points corresponding to the core feature curve of the first fault zone, then the stratigraphic value corresponding to the first intersection point is determined based on the stratigraphic values ​​of the two feature points closest to the first intersection point among the feature points corresponding to the core feature curve of the first fault zone. Here, each feature point corresponding to the core feature curve of the first fault zone is a feature point used to fit the core feature curve of the first fault zone.

[0077] For any second intersection point, if the second intersection point coincides with the fourth feature point among the feature points corresponding to the core feature curve of the first fault zone where the second intersection point is located, then the stratigraphic value of the fourth feature point is determined as the stratigraphic value corresponding to the second intersection point. If the second intersection point does not coincide with any of the feature points corresponding to the core feature curve of the first fault zone, then the stratigraphic value corresponding to the second intersection point is determined based on the stratigraphic values ​​of the two feature points closest to the second intersection point among the feature points corresponding to the core feature curve of the first fault zone.

[0078] In one possible implementation, the second determining module is configured to:

[0079] Determine the two feature points that are closest to the first reference point among the feature points corresponding to the first rising plate feature curve, wherein the two feature points include the fifth feature point and the sixth feature point;

[0080] On the first rising plate characteristic curve, determine the first curve point closest to the fifth characteristic point and the second curve point closest to the sixth characteristic point;

[0081] Based on the first arc length between the first curve point and the first reference point, the second arc length between the second curve point and the first reference point, and the layer value of the fifth feature point and the layer value of the sixth feature point, the layer value corresponding to the first reference point is determined.

[0082] The step of determining the stratigraphic value corresponding to the second reference point based on the stratigraphic values ​​of the two feature points closest to the second reference point among the feature points corresponding to the first descending disk feature curve includes:

[0083] Determine the two feature points that are closest to the second reference point among the feature points corresponding to the first descending disk feature curve, wherein the two feature points include the seventh feature point and the eighth feature point;

[0084] On the first descending disk characteristic curve, determine the third curve point that is closest to the seventh characteristic point and the fourth curve point that is closest to the eighth characteristic point;

[0085] Based on the third arc length between the third curve point and the second reference point, the fourth arc length between the fourth curve point and the second reference point, and the layer value of the seventh feature point and the layer value of the eighth feature point, the layer value corresponding to the second reference point is determined.

[0086] The step of determining the stratigraphic value corresponding to the first intersection point based on the stratigraphic values ​​of the two feature points closest to the first intersection point among the feature points corresponding to the core feature curve of the first fault zone includes:

[0087] Identify the two feature points that are closest to the first intersection point among the feature points corresponding to the core feature curve of the first fault zone, wherein the two feature points include the ninth feature point and the tenth feature point;

[0088] On the core characteristic curve of the first fault zone, determine the fifth curve point that is closest to the ninth characteristic point and the sixth curve point that is closest to the tenth characteristic point;

[0089] Based on the fifth arc length between the fifth curve point and the first intersection point, the sixth arc length between the sixth curve point and the first intersection point, the layer value of the ninth feature point and the layer value of the tenth feature point, the layer value corresponding to the first intersection point is determined.

[0090] The step of determining the stratigraphic value corresponding to the second intersection point based on the stratigraphic values ​​of the two feature points closest to the second intersection point among the feature points corresponding to the core feature curve of the first fault zone includes:

[0091] Identify the two feature points that are closest to the second intersection point among the feature points corresponding to the core feature curve of the first fault zone, wherein the two feature points include the eleventh feature point and the twelfth feature point;

[0092] On the core characteristic curve of the first fault zone, determine the seventh curve point that is closest to the eleventh characteristic point and the eighth curve point that is closest to the twelfth characteristic point;

[0093] Based on the seventh arc length between the seventh curve point and the second intersection point, the eighth arc length between the eighth curve point and the second intersection point, and the layer value of the eleventh feature point and the layer value of the twelfth feature point, the layer value corresponding to the second intersection point is determined.

[0094] Thirdly, a computer device is provided, the computer device including a processor and a memory, the memory storing at least one instruction, the instruction being loaded and executed by the processor to implement the first aspect and its possible implementations.

[0095] Fourthly, a computer-readable storage medium is provided, wherein at least one instruction is stored therein, the instruction being loaded and executed by a processor to implement the first aspect and a method thereof.

[0096] Fifthly, a computer program product is provided, the computer program product including at least one instruction, the at least one instruction being loaded and executed by a processor to implement the first aspect and its possible implementations.

[0097] The beneficial effects of the technical solutions provided in this application are:

[0098] In this embodiment, the uplift and downlift characteristic curves are determined from the fault zone planar data of the target fault, and reference point pairs are identified within these curves. Based on the stratigraphic values ​​of the reference point pairs, the fault displacement corresponding to each pair is determined. This method of calculating the fault displacement allows for the rapid determination of corresponding points between the uplift and downlift blocks using computer equipment, eliminating the need for technicians to visually determine corresponding feature points on the uplift and downlift blocks. This improves the efficiency of calculating the fault displacement between the uplift and downlift blocks. Attached Figure Description

[0099] To more clearly illustrate the technical solutions in the embodiments of this application, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0100] Figure 1 This is a flowchart of a method for determining fault displacement provided in an embodiment of this application;

[0101] Figure 2 This is the initial planar data map of the fracture zone of the target fault corresponding to the artificial intelligence fault prediction attributes provided in the embodiments of this application;

[0102] Figure 3 This is the initial planar data map of the fracture zone of the target fault corresponding to the coherence attributes provided in the embodiments of this application;

[0103] Figure 4This is a schematic diagram of the initial feature points in the initial planar data of the first fault zone provided in this application embodiment, where the core of the fault zone is likely located in the fault boundary interval;

[0104] Figure 5 This is a schematic diagram of the initial feature points in the initial planar data of the first fault zone provided in this application embodiment, where the fault zone core probability is located in the fault zone core region;

[0105] Figure 6 This is a schematic diagram illustrating the principle of determining the two disks of a fault, provided in an embodiment of this application.

[0106] Figure 7 This is a schematic diagram illustrating the principle of determining fault displacement provided in an embodiment of this application;

[0107] Figure 8 This is a schematic diagram illustrating the principle of determining fault displacement provided in an embodiment of this application;

[0108] Figure 9 This is a schematic diagram of a device for determining fault displacement provided in an embodiment of this application;

[0109] Figure 10 This is a structural block diagram of a server provided in an embodiment of this application. Detailed Implementation

[0110] To make the objectives, technical solutions, and advantages of this application clearer, the embodiments of this application will be described in further detail below with reference to the accompanying drawings.

[0111] This application provides a method for determining fault displacement, which can be implemented by a computer device, which can be a server or a terminal. The server can be a single server or a server cluster consisting of multiple servers.

[0112] A server may include a processor, memory, communication components, etc., with the processor connected to the memory and communication components respectively.

[0113] The processor can be a CPU (Central Processing Unit). The processor can be used to read instructions and process data, for example, acquiring fault zone planar data of the target fault, dividing the multiple feature points into multiple sets based on their horizontal coordinates, performing curve fitting based on the core feature points of the fault zone within each set, determining multiple reference points in the uplift and downlift feature curves at a first preset interval, determining the normal at each reference point and its intersection with the core feature curve of the fault zone, determining reference point pairs in the uplift and downlift feature curves, and determining the fault displacement corresponding to each reference point pair based on its stratigraphic value, etc.

[0114] Memory can include ROM (Read-Only Memory), RAM (Random Access Memory), CD-ROM (Compact Disc Read-Only Memory), hard disks, optical data storage devices, etc. Memory can be used for data storage, such as storing the acquired fault zone planar data of the target fault, storing the horizontal coordinates and stratigraphic values ​​of multiple feature points of the target fault, storing the ensemble data after dividing multiple feature points into multiple sets, storing the fitted core feature curves of the fault zone, storing the feature curves of the uplifted and downlifted blocks, and so on.

[0115] Communication components can be wired network connectors, WiFi (Wireless Fidelity) modules, Bluetooth modules, cellular network communication modules, etc. Communication components can be used to receive and transmit signals.

[0116] Figure 1 This is a flowchart illustrating a method for determining fault displacement according to an embodiment of this application. See also... Figure 1 This embodiment includes:

[0117] 101. Obtain the plane data of the fault zone of the target fault.

[0118] 102. Based on the horizontal coordinates of multiple feature points, divide the multiple feature points into multiple sets.

[0119] 103. For each set, curve fitting is performed based on the core feature points of the fault zone in the set to obtain the core feature curve of the fault zone. Curve fitting is performed based on the feature points other than the core feature points of the fault zone in the set to obtain the first feature curve. Based on the intersection of the first feature curve and the core feature curve of the fault zone, the first feature curve is divided into the second feature curve and the third feature curve. Based on the stratigraphic values ​​of each feature point corresponding to the second feature curve and the third feature curve, the uplifted block feature curve and the downlifted block feature curve are determined in the second feature curve and the third feature curve.

[0120] 104. For each ascending disk characteristic curve, determine multiple first reference points on the ascending disk characteristic curve at a first preset interval. Determine the first normal line of the ascending disk characteristic curve at each first reference point, and determine the first intersection point of each first normal line with the core characteristic curve of the fault zone. For each descending disk characteristic curve, determine multiple second reference points on the descending disk characteristic curve at a first preset interval. Determine the second normal line of the descending disk characteristic curve at each second reference point, and determine the second intersection point of each second normal line with the core characteristic curve of the fault zone. For each first intersection point, determine the second intersection point that is closest to the first intersection point and less than the second preset interval. Define the first reference point corresponding to the first intersection point and the second reference point corresponding to the second intersection point as a reference point pair.

[0121] 105. Based on the layer values ​​of multiple feature points, determine the layer value corresponding to each first reference point, each second reference point, each first intersection point, and each second intersection point.

[0122] 106. For each reference point pair, determine the first fault displacement based on the stratigraphic value of the first reference point and the stratigraphic value of the corresponding first intersection point in the reference point pair; determine the second fault displacement based on the stratigraphic value of the second reference point and the stratigraphic value of the corresponding second intersection point in the reference point pair; determine the sum of the first fault displacement and the second fault displacement to obtain the fault displacement corresponding to the reference point pair; determine the fault displacement corresponding to all reference point pairs as the fault displacement of the target fault.

[0123] Before performing step 101, it is necessary to first obtain the seismic data volume of the target area, and then process the seismic data volume accordingly to obtain the fault zone plane data of the target fault.

[0124] The seismic data volume is processed accordingly as follows:

[0125] The first step is to acquire the seismic data volume of the target area.

[0126] The seismic data volume includes the horizontal coordinates, horizon values, fault core probability, and attribute values ​​of multiple fault properties for each spatial point in the target area. The horizon value represents the depth of the spatial point underground; a larger horizon value indicates a deeper location. The unit of horizon value can be distance, such as meters, or time, such as seconds. The fault core probability represents the probability that a spatial point is located at the core of a fault zone. The fault core can be understood as the surface formed by the midpoints of the lines connecting corresponding points (adjacent points before the fault) on each pair of fault upside-downside blocks. Fault properties can be coherence properties, curvature properties, ant-like properties, energy properties, or AI-generated fault prediction properties, etc. AI-generated fault prediction properties are obtained by processing fault properties using a trained machine learning model. The input data for the model is the seismic data volume, and the output data is the attribute values ​​of the AI-generated fault prediction properties. The seismic data volume can be a three-dimensional matrix, where each element of the spatial matrix corresponds one-to-one with a spatial point in the target area. Each element in the spatial matrix can be a multidimensional vector data, which includes the horizontal coordinates of the corresponding spatial point, the stratigraphic value, the core probability of the fault zone, and the attribute values ​​of multiple fault attributes.

[0127] In practice, geological exploration blasting can be carried out near the surface at multiple exploration wells within the target area to artificially generate seismic waves. Seismic wave data is then received by surface-based detectors, and data collection and processing are performed to obtain the seismic data volume for the target area. Each spatial point can be numbered; for example, if there are 1000 spatial points within the target area, they can be sequentially numbered 0001, 0002, ..., 1000. When the layer value is in time, such as seconds, time-depth conversion can be performed on the layer value, and the depth corresponding to the layer value of each spatial point can be calculated using the seismic wave velocity.

[0128] Optionally, when using ground-based detectors to receive seismic wave data and collect and process data from geological exploration blasting, edge protection or structural guidance filtering methods can be used to preprocess the seismic data volume, making the seismic data volume clearer and more continuous.

[0129] The second step involves determining the fault corresponding to each spatial point based on the horizontal coordinates, layer values, and fault attribute values ​​of each spatial point in the seismic data volume. This process identifies multiple target spatial points corresponding to the first surface of the target fault and uses the horizontal coordinates and fault zone core probability of the target spatial points in the seismic data volume as the initial plane data of the fault zone of the target fault.

[0130] The first surface can be the upper or lower surface of the target fault, or it can be the middle layer of the target fault. Multiple fault attributes are considered, including at least two of the following: coherence attribute, curvature attribute, ant-like body attribute, energy attribute, and AI-generated fault prediction attribute. The initial planar data of the target fault zone can have multiple initial feature points, each corresponding one-to-one with each target spatial point. The initial planar data of the target fault zone can be a two-dimensional matrix, where each element of the planar matrix can be a multi-dimensional vector, including the horizontal coordinates of the corresponding target spatial point and the core probability of the fault zone.

[0131] In practice, for each fault attribute, the fault attribute values ​​of adjacent spatial points within the same fault are continuous. Therefore, the difference in fault attribute values ​​between adjacent spatial points within the same fault is generally small. When there is a sudden change in the difference between the attribute value of a certain fault attribute at any spatial point and the attribute value of the same fault attribute at its adjacent spatial points, it indicates that there is a fault between these two adjacent spatial points. Therefore, the distribution of faults can be described based on fault attributes. Simultaneously, multiple faults may exist within the target area, and correspondingly, the seismic data volume includes multiple faults. In this case, the faults can be numbered. For example, if there are three faults within the target area, these three faults can be numbered D1, D2, and D3 respectively. After obtaining the seismic data volume of the target area, the faults corresponding to each spatial point can be determined based on the horizontal coordinates and horizon values ​​of each spatial point in the seismic data volume. The numbers of spatial points belonging to the same fault are stored, and a fault-spatial point correspondence table is established based on the correspondence between fault numbers and the numbers of spatial points belonging to the same fault.

[0132] Table 1

[0133]

[0134] Technicians can select a target fault from multiple faults based on actual needs; for example, D3 can be designated as the target fault. After determining the target fault, spatial points located on the upper and lower surfaces of the target fault can be identified based on the horizontal coordinates and stratigraphic values ​​of the spatial points belonging to the target fault. Then, either the upper or lower surface of the target fault is selected as the first surface, and the spatial points belonging to the first surface are designated as target spatial points. The horizontal coordinates of the target spatial points and the core probability of the fault zone are then used to determine the initial planar data of the fault zone of the target fault. In this way, initial planar data of the fault zone of the target fault corresponding to at least two different fault attributes can be obtained.

[0135] The third step is to select the initial plane data of the first fault zone from the initial plane data of the fault zone of the target fault corresponding to each fault attribute.

[0136] Optionally, after obtaining the initial plane data of the fault zone of the target fault corresponding to each fault attribute, the data continuity value and signal-to-noise ratio (SNR) of each initial plane data point can be determined first. Then, for each initial plane data point, based on the first weighting coefficient corresponding to the data continuity value and the second weighting coefficient corresponding to the SNR, a weighted average of the data continuity value and SNR is determined as the selection score for the initial plane data point. Finally, among all the initial plane data points, the first initial plane data point with the highest selection score is selected.

[0137] The first weighting coefficient and the second weighting coefficient can be pre-designed by technicians according to actual needs, and the sum of the first weighting coefficient and the second weighting coefficient is 1.

[0138] In practice, the detection radius and continuity threshold can be preset.

[0139] Based on the detection radius, all initial feature points in the initial planar data can be traversed to obtain the signal-to-noise ratio (SNR) of the initial planar data of the target fault zone corresponding to each fault attribute. The process is as follows: Taking any initial feature point as the center, when drawing a circle with the detection radius, if there are no other initial feature points within that circle, the initial feature point at the center can be considered a noise point. Taking any initial feature point as the center, when drawing a circle with the detection radius, if there are at least one other initial feature point within that circle, the initial feature point at the center can be considered a valid signal point. Based on the ratio of the number of valid signal points to the number of noise points in the initial planar data, the SNR corresponding to the initial planar data of each fault zone is determined.

[0140] Simultaneously, based on the detection radius and continuity threshold, all initial feature points in the initial plane data of the target fault's fault zone can be traversed to obtain the data continuity value of the initial plane data of the target fault's fault zone corresponding to each fault attribute. The process is as follows: all initial feature points other than noise points are divided into multiple sets. In each set, when a circle is drawn with any initial feature point as the center and the detection radius as the center, the circle must include at least one other initial feature point from the same set and exclude any initial feature points from other sets. After dividing all initial feature points into multiple sets, the number of initial feature points in each set is determined, and the number of initial feature points in each set is compared with the continuity threshold. If the number of initial feature points in any set is greater than or equal to the continuity threshold, all initial feature points in that set are considered continuous data; if the number of initial feature points in any set is less than the continuity threshold, all initial feature points in that set are considered discontinuous data. For example, when the continuity threshold is preset to 10, if the number of initial feature points in a set is 10, then the data corresponding to these 10 initial feature points is considered continuous data; if the number of initial feature points in a set is 9, then the data corresponding to these 9 initial feature points is considered discontinuous data. After completing the above steps, the number of continuous data points in the initial planar data of the fault zone of the target fault corresponding to each fault attribute can be counted. Based on the ratio between the number of continuous data points and the total number of initial feature points, the data continuity value corresponding to the initial planar data of each fault zone can be determined.

[0141] like Figure 2 and Figure 3 As shown, Figure 2 The image shows the initial planar data of the fault zone of the target fault corresponding to the artificial intelligence fault prediction attributes. Figure 3 The diagram shows the initial planar data of the fault zone of the target fault corresponding to the coherence attributes. This is achieved through comparison. Figure 2 and Figure 3 It can be seen that the initial plane data of the fault zone of the target fault corresponding to the artificial intelligence fault prediction attribute has a high signal-to-noise ratio and a large data continuity value. Therefore, the initial plane data of the fault zone of the target fault corresponding to the artificial intelligence fault prediction attribute is preferred as the initial plane data of the first fault zone.

[0142] The fourth step is to determine the target horizontal coordinates of the corresponding fault zone core probability within the preset interval in the initial planar data of the first fault zone, identify the target spatial points with target horizontal coordinates among multiple target spatial points as feature points of the target fault, and determine the horizontal coordinates and stratigraphic values ​​of the feature points of the target fault as the fault zone planar data of the target fault.

[0143] Among them, the fault zone core probability is used to represent the probability that the initial feature point is located at the core of the fault zone, and the probability of each initial feature point being located at the core of the fault zone is in the interval [0, 1].

[0144] In implementation, the core probability of the fault zone at all initial feature points in the initial planar data of the first fault zone can be compared with a preset interval. If the core probability of any initial feature point falls within the preset interval, the initial feature point is retained; if the core probability of any initial feature point falls outside the preset interval, the initial feature point is deleted. After comparing all initial feature points, the target horizontal coordinates of the corresponding fault zone core probabilities within the preset interval are determined. Target spatial points with target horizontal coordinates among multiple target spatial points are identified as feature points of the target fault. The horizontal coordinates and stratigraphic values ​​of the feature points of the target fault are then used to determine the fault zone planar data of the target fault.

[0145] Optionally, the preset interval includes the core interval of the fault zone and the fault boundary interval. The value range of the core interval of the fault zone is [0.9, 1], and the value range of the fault boundary interval is [0.2, 0.3]. The target horizontal coordinate includes the first horizontal coordinate of the corresponding fault zone core probability within the core interval of the fault zone and the second horizontal coordinate of the corresponding fault zone core probability within the fault boundary interval.

[0146] Optionally, the target spatial point with the first horizontal coordinate among multiple target spatial points can be determined as the core feature point of the fault zone of the target fault.

[0147] like Figure 4 and Figure 5 As shown, Figure 4 and Figure 5 and Figure 2 The fracture zone area within the dashed box corresponds to this. Figure 4 The diagram shows initial feature points in the initial planar data of the first fault zone where the core probability of the fault zone lies within the fault boundary interval [0.2, 0.3]. Figure 5 The figure shows the initial feature points in the initial planar data of the first fault zone where the core probability of some fault zones is located in the core interval [0.9, 1].

[0148] After obtaining the fault zone planar data of the target fault, the fault displacement can be calculated, and the corresponding processing is as follows:

[0149] Regarding step 101 above: Obtain the plane data of the fault zone of the target fault.

[0150] The fault zone planar data includes the horizontal coordinates and stratigraphic values ​​of multiple feature points of the target fault, some of which are core feature points of the fault zone. This planar data comprises multiple matrix data, each corresponding to a feature point in the target fault. The matrix data may include the horizontal coordinates and stratigraphic values ​​of the corresponding feature point, and may also include the fault zone core probability of the corresponding feature point.

[0151] In practice, the fault zone plane data of the target fault can be obtained. After obtaining the fault zone plane data of the target fault, the feature points corresponding to the fault zone core probability of each feature point in the fault zone plane data can be determined as the fault zone core feature points based on the fault zone core probability of each feature point in the fault zone plane data. Then, the following processing steps 102-106 are performed to obtain the fault displacement data of the target fault.

[0152] Regarding step 102 above: Based on the horizontal coordinates of multiple feature points, divide the multiple feature points into multiple sets.

[0153] Specifically, for any feature point, there exists at least one other feature point within the set to which that feature point belongs, whose distance to that feature point is less than a distance threshold. The horizontal coordinates can be the latitude and longitude coordinates of each feature point.

[0154] In practice, the target fault may have multiple fault zones, each consisting of an uplifted block and a downlifted block. In this case, a distance threshold can be preset, and then the horizontal distance between any two feature points is calculated based on their horizontal coordinates. If the horizontal distance between two feature points is less than the preset distance threshold, they are considered to belong to the same fault zone and are grouped into the same set. If the horizontal distance is greater than or equal to the preset distance threshold, they are considered not to belong to the same fault zone and are not grouped into the same set. By calculating the horizontal distance between any two feature points, all feature points in the fault zone planar data of the target fault can be divided into multiple sets, each set corresponding to one fault zone.

[0155] Regarding step 103 above: For each set, curve fitting is performed based on the core feature points of the fault zone in the set to obtain the core feature curve of the fault zone. Curve fitting is performed based on feature points other than the core feature points of the fault zone in the set to obtain the first feature curve. Based on the intersection of the first feature curve and the core feature curve of the fault zone, the first feature curve is divided into a second feature curve and a third feature curve. Based on the stratigraphic values ​​of each feature point corresponding to the second feature curve and the third feature curve, the uplifted block feature curve and the downlifted block feature curve are determined in the second feature curve and the third feature curve.

[0156] Among them, the characteristic curves of the ascending and descending disks correspond to the ascending and descending disks of the fault zone, respectively.

[0157] In implementation, such as Figure 6 As shown, Figure 6 This is a longitudinal section diagram of a fault, and also a schematic diagram for determining the uplift and downlift blocks of a fault. The stratigraphic dip angle refers to the angle of inclination between the entire fault and the horizontal plane, as shown by angle α in the diagram. The fault dip angle refers to the angle of inclination between the fault section and the horizontal plane, as shown by angle β in the diagram. The stratigraphic dip angle α is less than the fault dip angle β, because... Figure 6 It is known that the core of the first fault zone and the adjacent uplifted and downlifted blocks constitute one fault zone, and the core of the second fault zone and the adjacent uplifted and downlifted blocks constitute another fault zone. Between these two adjacent fault zones, the dip angle γ between the line connecting the cores of the first and second fault zones and the horizontal plane will always be between the fault dip angle and the stratigraphic dip angle. Therefore, the stratigraphic values ​​of each characteristic point on the characteristic curve of the fault zone core are always greater than those on the characteristic curve of the downlifted block, and the stratigraphic values ​​of each characteristic point on the characteristic curve of the fault zone core are always less than those on the characteristic curve of the uplifted block. Therefore, by comparing the stratigraphic values ​​of each characteristic point corresponding to the second and third characteristic curves, the characteristic curve corresponding to the characteristic point with the larger stratigraphic value can be identified as the uplifted block characteristic curve, and the characteristic curve corresponding to the characteristic point with the smaller stratigraphic value can be identified as the downlifted block characteristic curve.

[0158] After determining the characteristic curves of the uplift and downlift sides, the fault displacement can be determined using the scatter interpolation method. The principle is as follows: characteristic points on the uplift and downlift side characteristic curves are expanded outwards, reference point pairs are determined based on the expansion results, and the fault displacement is determined at the core of the fault zone. For details of the scatter interpolation method, please refer to steps 104-106.

[0159] Regarding step 104 above: For each ascending disk characteristic curve, multiple first reference points are determined on the ascending disk characteristic curve at a first preset interval. The first normal line of the ascending disk characteristic curve at each first reference point is determined, and the first intersection point of each first normal line with the core characteristic curve of the fault zone is determined. For each descending disk characteristic curve, multiple second reference points are determined on the descending disk characteristic curve at a first preset interval. The second normal line of the descending disk characteristic curve at each second reference point is determined, and the second intersection point of each second normal line with the core characteristic curve of the fault zone is determined. For each first intersection point, the second intersection point that is closest to the first intersection point and less than the second preset interval is determined. The first reference point corresponding to the first intersection point and the second reference point corresponding to the second intersection point are determined as a reference point pair.

[0160] The first preset spacing and the second preset spacing can both be preset curve lengths.

[0161] In implementation, such as Figure 7 and Figure 8 As shown, Figure 7 This is a planar data map of the fault zone of the target fault. Figure 8 This is a longitudinal section diagram of the fault zone. In the fault zone planar data of the target fault, if the first normal of any point on the characteristic curve of the upside-side fault intersects the second normal of any point on the characteristic curve of the downside-side fault at the same point on the core characteristic curve of the fault zone, then the intersection of the first normal with the characteristic curve of the upside-side fault and the intersection of the second normal with the characteristic curve of the downside-side fault are generally considered to be a pair of corresponding points on the two sides of the fault. Therefore, the intersection of the characteristic curves of the upside-side and downside-side faults can be determined as the first first reference point and the first second reference point. Then, according to the curve direction of the characteristic curves of the upside-side and downside-side faults, a first reference point and a second reference point can be determined every arc length of a first preset interval. In this way, multiple first reference points and multiple second reference points can be determined in the characteristic curves of the upside-side and downside-side faults respectively. Then, at each first reference point in the ascending disk characteristic curve, the first normal line of the ascending disk characteristic curve can be determined, and the first intersection point of each first normal line with the core characteristic curve of the fault zone can be determined. Similarly, at each second reference point in the descending disk characteristic curve, the second normal line of the descending disk characteristic curve can be determined, and the second intersection point of each second normal line with the core characteristic curve of the fault zone can be determined. This results in multiple first intersection points and multiple second intersection points on the core characteristic curve of the fault zone. Finally, for each first intersection point, the second intersection point that is closest to the first intersection point and less than a second preset distance is determined. The first reference point corresponding to the first intersection point and the second reference point corresponding to the second intersection point are then defined as a reference point pair.

[0162] For any first intersection point, if there are two second intersection points whose distances from the first intersection point are less than a second preset spacing, and these two second intersection points are equidistant from the first intersection point, then the second intersection point that was not previously identified as a reference point pair can be identified as the second intersection point corresponding to the first intersection point. The first reference point corresponding to the first intersection point and the second reference point corresponding to the second intersection point can then be identified as a reference point pair. The vertical distance between the first reference point and the second reference point (i.e., the difference in the layer values ​​of the reference point pair) is the break distance between the reference point pairs.

[0163] Regarding step 105 above: Based on the layer values ​​of multiple feature points, determine the layer value corresponding to each first reference point, each second reference point, each first intersection point, and each second intersection point.

[0164] The following provides the stratigraphic values ​​for each first reference point, each second reference point, each first intersection point, and each second intersection point:

[0165] (1) Regarding the stratigraphic values ​​corresponding to each first reference point

[0166] For any first reference point, if the first reference point coincides with the first feature point among the feature points corresponding to the first ascending disk feature curve where the first reference point is located, then the stratigraphic value of the first feature point is determined as the stratigraphic value corresponding to the first reference point. If the first reference point does not coincide with any of the feature points corresponding to the first ascending disk feature curve, then the stratigraphic value corresponding to the first reference point is determined based on the stratigraphic values ​​of the two feature points closest to the first reference point among the feature points corresponding to the first ascending disk feature curve.

[0167] Among them, the feature points corresponding to the first rising plate feature curve are the feature points used to fit the first rising plate feature curve.

[0168] In practice, based on the layer values ​​of each feature point corresponding to the second and third feature curves, after determining the ascending disk feature curve and the descending disk feature curve in the second and third feature curves, the feature points included in the ascending disk feature curve can be marked to indicate that these feature points are the feature points used to fit the ascending disk feature curve.

[0169] Optionally, if none of the feature points corresponding to the first reference point and the first ascending disk characteristic curve coincide, then the two feature points closest to the first reference point among the feature points corresponding to the first ascending disk characteristic curve can be determined, wherein the two feature points include the fifth feature point and the sixth feature point. On the first ascending disk characteristic curve, the first curve point closest to the fifth feature point and the second curve point closest to the sixth feature point are determined. Then, based on the first arc length between the first curve point and the first reference point, the second arc length between the second curve point and the first reference point, and the stratigraphic values ​​of the fifth and sixth feature points, the stratigraphic value corresponding to the first reference point is determined.

[0170] The first curve point and the second curve point are both located in the characteristic curve of the rising plate.

[0171] In implementation, since the rising plate characteristic curve is obtained by fitting multiple feature points, some feature points may not lie on the curve. When the fifth and sixth reference points, which are closest to the first reference point, are not on the rising plate characteristic curve, the points closest to them can be determined on the curve, and these two points can be designated as the first curve point and the second curve point, respectively. After determining the first and second curve points, the formula first arc length / second arc length = ( Stratigraphic value of the first reference point - Stratigraphic value of the fifth feature point ) / ( Stratigraphic value of the first reference point - Stratigraphic value of the sixth feature point ), determine the stratigraphic value of the first reference point.

[0172] (2) Regarding the stratigraphic values ​​corresponding to each second reference point

[0173] For any second reference point, if the second feature point among the feature points corresponding to the first descending disk feature curve of the second reference point coincides with the second feature point, then the stratigraphic value of the second feature point is determined as the stratigraphic value corresponding to the second reference point. If the second reference point and the feature points corresponding to the first descending disk feature curve do not coincide, then the stratigraphic value corresponding to the second reference point is determined based on the stratigraphic values ​​of the two feature points closest to the second reference point among the feature points corresponding to the first descending disk feature curve.

[0174] Among them, the feature points corresponding to the first descending disk feature curve are the feature points used to fit the first descending disk feature curve.

[0175] In practice, based on the layer values ​​of each feature point corresponding to the second and third feature curves, after determining the ascending disk feature curve and the descending disk feature curve in the second and third feature curves, the feature points included in the descending disk feature curve can be marked to indicate that these feature points are the feature points used to fit the descending disk feature curve.

[0176] Optionally, if none of the feature points corresponding to the first reference point and the first ascending disk feature curve coincide, then the two feature points closest to the second reference point among the feature points corresponding to the first descending disk feature curve can be determined, where the two feature points include the seventh feature point and the eighth feature point. On the first descending disk feature curve, the third curve point closest to the seventh feature point and the fourth curve point closest to the eighth feature point are determined. Then, based on the third arc length between the third curve point and the second reference point, the fourth arc length between the fourth curve point and the second reference point, and the stratigraphic values ​​of the seventh and eighth feature points, the stratigraphic value corresponding to the second reference point is determined.

[0177] The third and fourth curve points are both located within the descending disk characteristic curve.

[0178] The process of determining the stratigraphic value corresponding to the second reference point is the same as the process of determining the stratigraphic value corresponding to the first reference point described above, and will not be repeated here.

[0179] (3) Regarding the stratigraphic value corresponding to each first intersection point

[0180] For any first intersection point, if the first intersection point coincides with the third feature point among the feature points corresponding to the core feature curve of the first fault zone where the first intersection point is located, then the stratigraphic value of the third feature point is determined as the stratigraphic value corresponding to the first intersection point. If the first intersection point does not coincide with any of the feature points corresponding to the core feature curve of the first fault zone, then the stratigraphic value corresponding to the first intersection point is determined based on the stratigraphic values ​​of the two feature points closest to the first intersection point among the feature points corresponding to the core feature curve of the first fault zone.

[0181] Among them, each feature point corresponding to the core feature curve of the first fault zone is a feature point used to fit the core feature curve of the first fault zone.

[0182] Optionally, if none of the feature points corresponding to the first intersection point and the core feature curve of the first fault zone coincide, then the two feature points closest to the first intersection point among the feature points corresponding to the core feature curve of the first fault zone can be determined, where the two feature points include the ninth feature point and the tenth feature point. Then, on the core feature curve of the first fault zone, the fifth curve point closest to the ninth feature point and the sixth curve point closest to the tenth feature point are determined. Then, based on the fifth arc length between the fifth curve point and the first intersection point, the sixth arc length between the sixth curve point and the first intersection point, and the stratigraphic values ​​of the ninth and tenth feature points, the stratigraphic value corresponding to the first intersection point is determined.

[0183] The fifth and sixth curve points are both located within the core characteristic curves of the fault zone.

[0184] The process of determining the stratigraphic value corresponding to the first intersection point is the same as the process of determining the stratigraphic value corresponding to the first reference point described above, and will not be repeated here.

[0185] (4) Regarding the stratigraphic values ​​corresponding to each second intersection point

[0186] For any second intersection point, if the fourth feature point among the feature points corresponding to the core feature curve of the first fault zone where the second intersection point is located coincides with the fourth feature point, then the stratigraphic value of the fourth feature point is determined as the stratigraphic value corresponding to the second intersection point. If the second intersection point and the feature points corresponding to the core feature curve of the first fault zone do not coincide, then the stratigraphic value corresponding to the second intersection point is determined based on the stratigraphic values ​​of the two feature points closest to the second intersection point among the feature points corresponding to the core feature curve of the first fault zone.

[0187] Optionally, if none of the feature points corresponding to the first intersection point and the core feature curve of the first fault zone coincide, then the two feature points closest to the second intersection point among the feature points corresponding to the core feature curve of the first fault zone can be determined, where the two feature points include the eleventh and twelfth feature points. On the core feature curve of the first fault zone, the seventh curve point closest to the eleventh feature point and the eighth curve point closest to the twelfth feature point are determined. Based on the seventh arc length between the seventh curve point and the second intersection point, the eighth arc length between the eighth curve point and the second intersection point, and the stratigraphic values ​​of the eleventh and twelfth feature points, the stratigraphic value corresponding to the second intersection point is determined.

[0188] The seventh and eighth curve points are both located within the core characteristic curves of the fault zone.

[0189] The process of determining the stratigraphic value corresponding to the second intersection point is the same as the process of determining the stratigraphic value corresponding to the first reference point, and will not be described again here.

[0190] Regarding step 106 above: For each reference point pair, based on the stratigraphic value of the first reference point in the reference point pair and the stratigraphic value of the corresponding first intersection point, determine the first fault displacement; based on the stratigraphic value of the second reference point in the reference point pair and the stratigraphic value of the corresponding second intersection point, determine the second fault displacement; determine the sum of the first fault displacement and the second fault displacement to obtain the fault displacement corresponding to the reference point pair; determine the fault displacement corresponding to all reference point pairs as the fault displacement of the target fault.

[0191] The first fault displacement can be used to represent the vertical distance between each point in the uplifted block of the fault and the core of the fault zone, while the second fault displacement can be used to represent the vertical distance between each point in the downlifted block of the fault and the core of the fault zone.

[0192] Optionally, after completing step 104, the stratigraphic values ​​corresponding to each first reference point and each second reference point can be determined based on the stratigraphic values ​​of multiple feature points. The stratigraphic value corresponding to the first reference point is assigned to the first intersection point corresponding to the first reference point, and the stratigraphic value corresponding to the second reference point is assigned to the second intersection point corresponding to the second reference point. For each first intersection point, the second intersection point that is closest to the first intersection point and less than a second preset distance can be determined. The first reference point corresponding to the first intersection point and the second reference point corresponding to the second intersection point are defined as a reference point pair. For each reference point pair, the difference between the stratigraphic values ​​of the first and second intersection points corresponding to each reference point pair is directly determined at the core of the fault zone to obtain the fault displacement corresponding to each reference point pair. Finally, the fault displacements corresponding to all reference point pairs are determined as the fault displacement of the target fault.

[0193] In this embodiment, the uplift and downlift characteristic curves are determined from the fault zone planar data of the target fault, and reference point pairs are identified within these curves. Based on the stratigraphic values ​​of the reference point pairs, the fault displacement corresponding to each pair is determined. This method of calculating the fault displacement allows for the rapid determination of corresponding points between the uplift and downlift blocks using computer equipment, eliminating the need for technicians to visually determine corresponding feature points on the uplift and downlift blocks. This improves the efficiency of calculating the fault displacement between the uplift and downlift blocks.

[0194] This application provides an apparatus for determining fault displacement. This apparatus can be the computer device described in the above embodiments, such as... Figure 9 As shown, the device includes:

[0195] The acquisition module 910 is used to acquire the fracture zone plane data of the target fault, wherein the fracture zone plane data includes the horizontal coordinates and stratigraphic values ​​of multiple feature points of the target fault, and some of the feature points are core feature points of the fracture zone.

[0196] The partitioning module 920 is used to partition the multiple feature points into multiple sets based on the horizontal coordinates of the multiple feature points, wherein, for any feature point, there is at least one other feature point in the set to which the feature point belongs, and the distance between the feature point and the feature point is less than a distance threshold.

[0197] The fitting module 930 is used to perform curve fitting based on the core feature points of the fault zone in each set to obtain the core feature curve of the fault zone, and to perform curve fitting based on feature points other than the core feature points of the fault zone in the set to obtain a first feature curve. Based on the intersection of the first feature curve and the core feature curve of the fault zone, the first feature curve is divided into a second feature curve and a third feature curve. Based on the stratigraphic values ​​of each feature point corresponding to the second feature curve and the third feature curve, the rising disk feature curve and the falling disk feature curve are determined in the second feature curve and the third feature curve.

[0198] The first determining module 940 is used to determine multiple first reference points on each rising disk feature curve at a first preset interval, determine the first normal of the rising disk feature curve at each first reference point, and determine the first intersection point of each first normal with the core feature curve of the fracture zone.

[0199] For each descending disk characteristic curve, multiple second reference points are determined on the descending disk characteristic curve at the first preset interval, and a second normal line of the descending disk characteristic curve at each second reference point is determined, and a second intersection point of each second normal line with the core characteristic curve of the fracture zone is determined.

[0200] For each first intersection point, determine the second intersection point that is closest to the first intersection point and less than the second preset distance, and determine the first reference point corresponding to the first intersection point and the second reference point corresponding to the second intersection point as a reference point pair;

[0201] Based on the layer values ​​of the multiple feature points, determine the layer value corresponding to each first reference point, each second reference point, each first intersection point, and each second intersection point;

[0202] The second determining module 950 is used to determine a first discontinuity for each reference point pair based on the layer value of the first reference point in the reference point pair and the layer value of the corresponding first intersection point, determine a second discontinuity based on the layer value of the second reference point in the reference point pair and the layer value of the corresponding second intersection point, and determine the sum of the first discontinuity and the second discontinuity to obtain the discontinuity corresponding to the reference point pair.

[0203] In one possible implementation, the acquisition module 910 is further configured to:

[0204] Acquire seismic data volume for a target area, wherein the seismic data volume includes the horizontal coordinates, layer values, fault zone core probability, and attribute values ​​of multiple fault attributes for each spatial point in the target area, wherein the fault zone core probability is used to represent the probability that the spatial point is located at the core of the fault zone;

[0205] The first determining module 940 is used for:

[0206] For each fault attribute, based on the horizontal coordinates, layer values, and attribute values ​​of each spatial point in the seismic data volume, the fault corresponding to each spatial point is determined, and multiple target spatial points corresponding to the first surface of the target fault are determined. The horizontal coordinates and fault zone core probability of the target spatial points in the seismic data volume are determined as the initial plane data of the fault zone of the target fault, wherein the first surface is the upper or lower surface of the target fault.

[0207] In the initial plane data of the fault zone of the target fault corresponding to each fault attribute, the initial plane data of the first fault zone is selected;

[0208] In the initial planar data of the first fault zone, the target horizontal coordinates of the corresponding fault zone core probability within a preset interval are determined. The target spatial points with the target horizontal coordinates among multiple target spatial points are determined as feature points of the target fault. The horizontal coordinates and stratigraphic values ​​of the feature points of the target fault are determined as the fault zone planar data of the target fault.

[0209] In one possible implementation, the plurality of fault attributes includes at least two of the following: coherence attributes, curvature attributes, ant body attributes, energy attributes, and artificial intelligence fault prediction attributes.

[0210] In one possible implementation, the first determining module 940 is configured to:

[0211] Obtain the initial planar data of the fracture zone of the target fault corresponding to each fault attribute;

[0212] Determine the data continuity value and signal-to-noise ratio corresponding to the initial planar data of each fault zone;

[0213] For each initial plane data of a fault zone, a weighted average of the data continuity value and the signal-to-noise ratio corresponding to the initial plane data of the fault zone is determined based on the first weighting coefficient corresponding to the data continuity value and the second weighting coefficient corresponding to the signal-to-noise ratio, and is used as the selection score corresponding to the initial plane data of the fault zone.

[0214] Among all the initial plane data of the fault zones, the first initial plane data of the fault zone with the highest selection score is selected.

[0215] In one possible implementation, the preset interval includes a fault zone core interval and a fault boundary interval, wherein the value range of the fault zone core interval is [0.9, 1], the value range of the fault boundary interval is [0.2, 0.3], and the target horizontal coordinate includes a first horizontal coordinate corresponding to the fault zone core probability within the fault zone core interval and a second horizontal coordinate corresponding to the fault zone core probability within the fault boundary interval.

[0216] The first determining module is further configured to:

[0217] The target spatial point with the first horizontal coordinate among multiple target spatial points is determined as the core feature point of the fracture zone of the target fault.

[0218] In one possible implementation, the second determining module 950 is configured to:

[0219] For any first reference point, if the first reference point coincides with the first feature point among the feature points corresponding to the first ascending disk feature curve where the first reference point is located, then the layer value of the first feature point is determined as the layer value corresponding to the first reference point. If the first reference point does not coincide with any of the feature points corresponding to the first ascending disk feature curve, then the layer value corresponding to the first reference point is determined based on the layer values ​​of the two feature points closest to the first reference point among the feature points corresponding to the first ascending disk feature curve. Here, each feature point corresponding to the first ascending disk feature curve is a feature point used to fit the first ascending disk feature curve.

[0220] For any second reference point, if the second reference point coincides with the second feature point among the feature points corresponding to the first descending disk feature curve where the second reference point is located, then the layer value of the second feature point is determined as the layer value corresponding to the second reference point. If the second reference point does not coincide with any of the feature points corresponding to the first descending disk feature curve, then the layer value corresponding to the second reference point is determined based on the layer values ​​of the two feature points closest to the second reference point among the feature points corresponding to the first descending disk feature curve. Here, each feature point corresponding to the first descending disk feature curve is a feature point used to fit the first descending disk feature curve.

[0221] For any first intersection point, if the first intersection point coincides with the third feature point among the feature points corresponding to the core feature curve of the first fault zone where the first intersection point is located, then the stratigraphic value of the third feature point is determined as the stratigraphic value corresponding to the first intersection point. If the first intersection point does not coincide with any of the feature points corresponding to the core feature curve of the first fault zone, then the stratigraphic value corresponding to the first intersection point is determined based on the stratigraphic values ​​of the two feature points closest to the first intersection point among the feature points corresponding to the core feature curve of the first fault zone. Here, each feature point corresponding to the core feature curve of the first fault zone is a feature point used to fit the core feature curve of the first fault zone.

[0222] For any second intersection point, if the second intersection point coincides with the fourth feature point among the feature points corresponding to the core feature curve of the first fault zone where the second intersection point is located, then the stratigraphic value of the fourth feature point is determined as the stratigraphic value corresponding to the second intersection point. If the second intersection point does not coincide with any of the feature points corresponding to the core feature curve of the first fault zone, then the stratigraphic value corresponding to the second intersection point is determined based on the stratigraphic values ​​of the two feature points closest to the second intersection point among the feature points corresponding to the core feature curve of the first fault zone.

[0223] In one possible implementation, the second determining module 950 is configured to:

[0224] Determine the two feature points that are closest to the first reference point among the feature points corresponding to the first rising plate feature curve, wherein the two feature points include the fifth feature point and the sixth feature point;

[0225] On the first rising plate characteristic curve, determine the first curve point closest to the fifth characteristic point and the second curve point closest to the sixth characteristic point;

[0226] Based on the first arc length between the first curve point and the first reference point, the second arc length between the second curve point and the first reference point, and the layer value of the fifth feature point and the layer value of the sixth feature point, the layer value corresponding to the first reference point is determined.

[0227] The step of determining the stratigraphic value corresponding to the second reference point based on the stratigraphic values ​​of the two feature points closest to the second reference point among the feature points corresponding to the first descending disk feature curve includes:

[0228] Determine the two feature points that are closest to the second reference point among the feature points corresponding to the first descending disk feature curve, wherein the two feature points include the seventh feature point and the eighth feature point;

[0229] On the first descending disk characteristic curve, determine the third curve point that is closest to the seventh characteristic point and the fourth curve point that is closest to the eighth characteristic point;

[0230] Based on the third arc length between the third curve point and the second reference point, the fourth arc length between the fourth curve point and the second reference point, and the layer value of the seventh feature point and the layer value of the eighth feature point, the layer value corresponding to the second reference point is determined.

[0231] The step of determining the stratigraphic value corresponding to the first intersection point based on the stratigraphic values ​​of the two feature points closest to the first intersection point among the feature points corresponding to the core feature curve of the first fault zone includes:

[0232] Identify the two feature points that are closest to the first intersection point among the feature points corresponding to the core feature curve of the first fault zone, wherein the two feature points include the ninth feature point and the tenth feature point;

[0233] On the core characteristic curve of the first fault zone, determine the fifth curve point that is closest to the ninth characteristic point and the sixth curve point that is closest to the tenth characteristic point;

[0234] Based on the fifth arc length between the fifth curve point and the first intersection point, the sixth arc length between the sixth curve point and the first intersection point, the layer value of the ninth feature point and the layer value of the tenth feature point, the layer value corresponding to the first intersection point is determined.

[0235] The step of determining the stratigraphic value corresponding to the second intersection point based on the stratigraphic values ​​of the two feature points closest to the second intersection point among the feature points corresponding to the core feature curve of the first fault zone includes:

[0236] Identify the two feature points that are closest to the second intersection point among the feature points corresponding to the core feature curve of the first fault zone, wherein the two feature points include the eleventh feature point and the twelfth feature point;

[0237] On the core characteristic curve of the first fault zone, determine the seventh curve point that is closest to the eleventh characteristic point and the eighth curve point that is closest to the twelfth characteristic point;

[0238] Based on the seventh arc length between the seventh curve point and the second intersection point, the eighth arc length between the eighth curve point and the second intersection point, and the layer value of the eleventh feature point and the layer value of the twelfth feature point, the layer value corresponding to the second intersection point is determined.

[0239] It should be noted that the above embodiments of the device for determining fault displacement are only illustrative examples of the division of the functional modules described above. In practical applications, the functions described above can be assigned to different functional modules as needed, that is, the internal structure of the device can be divided into different functional modules to complete all or part of the functions described above. Furthermore, the device for determining fault displacement and the method embodiments for determining fault displacement provided in the above embodiments belong to the same concept, and their specific implementation process is detailed in the method embodiments, which will not be repeated here.

[0240] Figure 10 This is a schematic diagram of the structure of a computer device provided in an embodiment of this application. The computer device can be a server or a terminal. The computer device 1000 can vary significantly due to differences in configuration or performance, and may include one or more processors 1001 and one or more memories 1002. The memories 1002 store at least one instruction, which is loaded and executed by the processors 1001 to implement the methods provided in the above-described method embodiments. Of course, the server may also have wired or wireless network interfaces, a keyboard, and input / output interfaces for input and output. The server may also include other components for implementing device functions, which will not be elaborated upon here.

[0241] In an exemplary embodiment, a computer-readable storage medium is also provided, such as a memory including instructions executable by a processor in a terminal to perform the method for identifying seismic multiples in the above embodiments. This computer-readable storage medium may be non-transitory. For example, the computer-readable storage medium may be a ROM, RAM, CD-ROM, magnetic tape, floppy disk, or optical data storage device, etc.

[0242] Those skilled in the art will understand that all or part of the steps of the above embodiments can be implemented by hardware or by a program instructing related hardware. The program can be stored in a computer-readable storage medium, such as a read-only memory, a disk, or an optical disk.

[0243] The above description is merely an optional embodiment of this application and is not intended to limit this application. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this application should be included within the protection scope of this application.

Claims

1. A method for determining fault displacement, characterized in that, The method includes: Obtain the fault zone plane data of the target fault, wherein the fault zone plane data includes the horizontal coordinates and stratigraphic values ​​of multiple feature points of the target fault, and some of the feature points are core feature points of the fault zone. Based on the horizontal coordinates of the multiple feature points, the multiple feature points are divided into multiple sets, wherein for any feature point, there is at least one other feature point in the set to which the feature point belongs that the distance between the feature point and the feature point is less than a distance threshold. For each set, curve fitting is performed based on the core feature points of the fault zone in the set to obtain the core feature curve of the fault zone. Curve fitting is performed based on feature points other than the core feature points of the fault zone in the set to obtain the first feature curve. Based on the intersection of the first feature curve and the core feature curve of the fault zone, the first feature curve is divided into a second feature curve and a third feature curve. Based on the stratigraphic values ​​of each feature point corresponding to the second feature curve and the third feature curve, the uplifted block feature curve and the downlifted block feature curve are determined in the second feature curve and the third feature curve. For each rising disk characteristic curve, multiple first reference points are determined on the rising disk characteristic curve at a first preset interval, and a first normal line of the rising disk characteristic curve at each first reference point is determined, and a first intersection point of each first normal line with the core characteristic curve of the fracture zone is determined. For each descending disk characteristic curve, multiple second reference points are determined on the descending disk characteristic curve at the first preset interval, and a second normal line of the descending disk characteristic curve at each second reference point is determined, and a second intersection point of each second normal line with the core characteristic curve of the fracture zone is determined. For each first intersection point, determine the second intersection point that is closest to the first intersection point and less than the second preset distance, and determine the first reference point corresponding to the first intersection point and the second reference point corresponding to the second intersection point as a reference point pair; Based on the layer values ​​of the multiple feature points, determine the layer value corresponding to each first reference point, each second reference point, each first intersection point, and each second intersection point; For each reference point pair, a first dislocation is determined based on the stratigraphic value of the first reference point in the reference point pair and the stratigraphic value of the corresponding first intersection point. A second dislocation is determined based on the stratigraphic value of the second reference point in the reference point pair and the stratigraphic value of the corresponding second intersection point. The sum of the first dislocation and the second dislocation is determined to obtain the dislocation corresponding to the reference point pair. The fault displacements corresponding to all reference points are determined as the fault displacements of the target fault.

2. The method according to claim 1, characterized in that, Before acquiring the fault zone plane data of the target fault, the method further includes: Acquire seismic data volume for a target area, wherein the seismic data volume includes the horizontal coordinates, layer values, fault zone core probability, and attribute values ​​of multiple fault attributes for each spatial point in the target area, wherein the fault zone core probability is used to represent the probability that the spatial point is located at the core of the fault zone; For each fault attribute, based on the horizontal coordinates, layer values, and attribute values ​​of each spatial point in the seismic data volume, the fault corresponding to each spatial point is determined, and multiple target spatial points corresponding to the first surface of the target fault are determined. The horizontal coordinates and fault zone core probability of the target spatial points in the seismic data volume are determined as the initial plane data of the fault zone of the target fault, wherein the first surface is the upper or lower surface of the target fault. In the initial plane data of the fault zone of the target fault corresponding to each fault attribute, the initial plane data of the first fault zone is selected; In the initial planar data of the first fault zone, the target horizontal coordinates of the corresponding fault zone core probability within a preset interval are determined. The target spatial points with the target horizontal coordinates among multiple target spatial points are determined as feature points of the target fault. The horizontal coordinates and stratigraphic values ​​of the feature points of the target fault are determined as the fault zone planar data of the target fault.

3. The method according to claim 2, characterized in that, The aforementioned multiple fault attributes include at least two of the following: coherence attributes, curvature attributes, ant body attributes, energy attributes, and artificial intelligence fault prediction attributes.

4. The method according to claim 2, characterized in that, The step of selecting the first initial plane data of the fault zone from the initial plane data of the target fault corresponding to each attribute includes: Obtain the initial planar data of the fracture zone of the target fault corresponding to each fault attribute; Determine the data continuity value and signal-to-noise ratio corresponding to the initial planar data of each fault zone; For each initial plane data of a fault zone, a weighted average of the data continuity value and the signal-to-noise ratio corresponding to the initial plane data of the fault zone is determined based on the first weighting coefficient corresponding to the data continuity value and the second weighting coefficient corresponding to the signal-to-noise ratio, and is used as the selection score corresponding to the initial plane data of the fault zone. Among all the initial plane data of the fault zones, the first initial plane data of the fault zone with the highest selection score is selected.

5. The method according to claim 2, characterized in that, The preset interval includes a fault zone core interval and a fault boundary interval. The value range of the fault zone core interval is [0.9, 1], and the value range of the fault boundary interval is [0.2, 0.3]. The target horizontal coordinate includes the first horizontal coordinate of the corresponding fault zone core probability within the fault zone core interval and the second horizontal coordinate of the corresponding fault zone core probability within the fault boundary interval. The method further includes: The target spatial point with the first horizontal coordinate among multiple target spatial points is determined as the core feature point of the fracture zone of the target fault.

6. The method according to claim 1, characterized in that, The determination of the layer value corresponding to each first reference point, each second reference point, each first intersection point, and each second intersection point based on the layer values ​​of the multiple feature points includes: For any first reference point, if the first reference point coincides with the first feature point among the feature points corresponding to the first ascending disk feature curve where the first reference point is located, then the layer value of the first feature point is determined as the layer value corresponding to the first reference point. If the first reference point does not coincide with any of the feature points corresponding to the first ascending disk feature curve, then the layer value corresponding to the first reference point is determined based on the layer values ​​of the two feature points closest to the first reference point among the feature points corresponding to the first ascending disk feature curve. Here, each feature point corresponding to the first ascending disk feature curve is a feature point used to fit the first ascending disk feature curve. For any second reference point, if the second reference point coincides with the second feature point among the feature points corresponding to the first descending disk feature curve where the second reference point is located, then the layer value of the second feature point is determined as the layer value corresponding to the second reference point. If the second reference point does not coincide with any of the feature points corresponding to the first descending disk feature curve, then the layer value corresponding to the second reference point is determined based on the layer values ​​of the two feature points closest to the second reference point among the feature points corresponding to the first descending disk feature curve. Here, each feature point corresponding to the first descending disk feature curve is a feature point used to fit the first descending disk feature curve. For any first intersection point, if the first intersection point coincides with the third feature point among the feature points corresponding to the core feature curve of the first fault zone where the first intersection point is located, then the stratigraphic value of the third feature point is determined as the stratigraphic value corresponding to the first intersection point. If the first intersection point does not coincide with any of the feature points corresponding to the core feature curve of the first fault zone, then the stratigraphic value corresponding to the first intersection point is determined based on the stratigraphic values ​​of the two feature points closest to the first intersection point among the feature points corresponding to the core feature curve of the first fault zone. Here, each feature point corresponding to the core feature curve of the first fault zone is a feature point used to fit the core feature curve of the first fault zone. For any second intersection point, if the second intersection point coincides with the fourth feature point among the feature points corresponding to the core feature curve of the first fault zone where the second intersection point is located, then the stratigraphic value of the fourth feature point is determined as the stratigraphic value corresponding to the second intersection point. If the second intersection point does not coincide with any of the feature points corresponding to the core feature curve of the first fault zone, then the stratigraphic value corresponding to the second intersection point is determined based on the stratigraphic values ​​of the two feature points closest to the second intersection point among the feature points corresponding to the core feature curve of the first fault zone.

7. The method according to claim 6, characterized in that, The step of determining the stratigraphic value corresponding to the first reference point based on the stratigraphic values ​​of the two feature points closest to the first reference point among the feature points corresponding to the first ascending disc feature curve includes: Determine the two feature points that are closest to the first reference point among the feature points corresponding to the first rising plate feature curve, wherein the two feature points include the fifth feature point and the sixth feature point; On the first rising plate characteristic curve, determine the first curve point closest to the fifth characteristic point and the second curve point closest to the sixth characteristic point; Based on the first arc length between the first curve point and the first reference point, the second arc length between the second curve point and the first reference point, and the layer value of the fifth feature point and the layer value of the sixth feature point, the layer value corresponding to the first reference point is determined. The step of determining the stratigraphic value corresponding to the second reference point based on the stratigraphic values ​​of the two feature points closest to the second reference point among the feature points corresponding to the first descending disk feature curve includes: Determine the two feature points that are closest to the second reference point among the feature points corresponding to the first descending disk feature curve, wherein the two feature points include the seventh feature point and the eighth feature point; On the first descending disk characteristic curve, determine the third curve point that is closest to the seventh characteristic point and the fourth curve point that is closest to the eighth characteristic point; Based on the third arc length between the third curve point and the second reference point, the fourth arc length between the fourth curve point and the second reference point, and the layer value of the seventh feature point and the layer value of the eighth feature point, the layer value corresponding to the second reference point is determined. The step of determining the stratigraphic value corresponding to the first intersection point based on the stratigraphic values ​​of the two feature points closest to the first intersection point among the feature points corresponding to the core feature curve of the first fault zone includes: Identify the two feature points that are closest to the first intersection point among the feature points corresponding to the core feature curve of the first fault zone, wherein the two feature points include the ninth feature point and the tenth feature point; On the core characteristic curve of the first fault zone, determine the fifth curve point that is closest to the ninth characteristic point and the sixth curve point that is closest to the tenth characteristic point; Based on the fifth arc length between the fifth curve point and the first intersection point, the sixth arc length between the sixth curve point and the first intersection point, the layer value of the ninth feature point and the layer value of the tenth feature point, the layer value corresponding to the first intersection point is determined. The step of determining the stratigraphic value corresponding to the second intersection point based on the stratigraphic values ​​of the two feature points closest to the second intersection point among the feature points corresponding to the core feature curve of the first fault zone includes: Identify the two feature points that are closest to the second intersection point among the feature points corresponding to the core feature curve of the first fault zone, wherein the two feature points include the eleventh feature point and the twelfth feature point; On the core characteristic curve of the first fault zone, determine the seventh curve point that is closest to the eleventh characteristic point and the eighth curve point that is closest to the twelfth characteristic point; Based on the seventh arc length between the seventh curve point and the second intersection point, the eighth arc length between the eighth curve point and the second intersection point, and the layer value of the eleventh feature point and the layer value of the twelfth feature point, the layer value corresponding to the second intersection point is determined.

8. A device for determining fault displacement, characterized in that, The device includes: The acquisition module is used to acquire the fracture zone plane data of the target fault, wherein the fracture zone plane data includes the horizontal coordinates and stratigraphic values ​​of multiple feature points of the target fault, and some of the feature points are core feature points of the fracture zone. The partitioning module is used to divide the multiple feature points into multiple sets based on the horizontal coordinates of the multiple feature points, wherein, for any feature point, there is at least one other feature point in the set to which the feature point belongs, and the distance between the feature point and the feature point is less than a distance threshold. The fitting module is used to perform curve fitting based on the core feature points of the fault zone in each set to obtain the core feature curve of the fault zone, and to perform curve fitting based on feature points other than the core feature points of the fault zone in the set to obtain the first feature curve. Based on the intersection of the first feature curve and the core feature curve of the fault zone, the first feature curve is divided into a second feature curve and a third feature curve. Based on the stratigraphic values ​​of each feature point corresponding to the second feature curve and the third feature curve, the rising block feature curve and the falling block feature curve are determined in the second feature curve and the third feature curve. The first determining module is used to determine multiple first reference points on each rising disk feature curve at a first preset interval, determine the first normal line of the rising disk feature curve at each first reference point, and determine the first intersection point of each first normal line with the core feature curve of the fracture zone. For each descending disk characteristic curve, multiple second reference points are determined on the descending disk characteristic curve at the first preset interval, and a second normal line of the descending disk characteristic curve at each second reference point is determined, and a second intersection point of each second normal line with the core characteristic curve of the fracture zone is determined. For each first intersection point, determine the second intersection point that is closest to the first intersection point and less than the second preset distance, and determine the first reference point corresponding to the first intersection point and the second reference point corresponding to the second intersection point as a reference point pair; Based on the layer values ​​of the multiple feature points, determine the layer value corresponding to each first reference point, each second reference point, each first intersection point, and each second intersection point; The second determining module is used to determine a first discontinuity for each reference point pair based on the layer value of the first reference point in the reference point pair and the layer value of the corresponding first intersection point, determine a second discontinuity based on the layer value of the second reference point in the reference point pair and the layer value of the corresponding second intersection point, and determine the sum of the first discontinuity and the second discontinuity to obtain the discontinuity corresponding to the reference point pair.

9. A computer device, characterized in that, The computer device includes a processor and a memory, the memory storing at least one instruction, which is loaded and executed by the processor to perform the operations performed by the method for determining fault displacement as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, The storage medium stores at least one instruction, which is loaded and executed by a processor to perform the operations of the method for determining fault displacement as described in any one of claims 1 to 7.

Images