Fault plane analysis method and device based on tendency true curvature, equipment and medium

Abstract

The invention discloses a fault plane analysis method, device and equipment based on tendency true curvature, and a medium. The method comprises the following steps: acquiring seismic horizon data of a target reservoir on a depth domain; determining the tendency true curvature of a local curved surface corresponding to each common center point in the seismic horizon data; according to the tendency true curvature of the local curved surface corresponding to each common center point, a first curvature combination, a second curvature combination, and a common center point corresponding to the first curvature combination and a common center point corresponding to the second curvature combination are determined, the first curvature combination comprises at least two positive tendency true curvature extreme points, and the second curvature combination comprises at least two positive tendency true curvature extreme points; the second curvature combination comprises at least two negative tendency true curvature extreme points; and determining fault plane distribution data in the target reservoir according to the common center point corresponding to the first curvature combination and the common center point corresponding to the second curvature combination. According to the technical scheme provided by the embodiment of the invention, the technical effect of automatically and accurately determining the fault surface distribution data can be achieved.

Description

Technical Field

[0001] The invention relates to the technical field of petroleum exploration, and in particular to a fault plane analysis method, device, equipment and medium based on dip true curvature. Background Art

[0002] In petroleum exploration, the important issue that affects the exploration and development of oil and gas fields is fault structure. The existing technology often uses conventional 3D seismic interpretation. Conventional 3D seismic interpretation is basically carried out on vertical seismic profiles. It requires interpreting each survey line in the work area one by one to determine the plane distribution of faults and stratigraphic features, so as to use manual analysis, judgment and interpretation on the plane. However, this method is cumbersome and requires manual analysis, and it is impossible to automatically and accurately determine the fault surface distribution data. Summary of the invention

[0003] The present invention provides a fault plane analysis method, device, equipment and medium based on dip true curvature to solve the problem that the prior art cannot automatically and accurately determine fault plane distribution data.

[0004] According to one aspect of the present invention, a fault plane analysis method based on dip true curvature is provided, the method comprising:

[0005] Acquire seismic layer data of the target reservoir in the depth domain;

[0006] Determine the true curvature of the local surface corresponding to each common center point in the seismic horizon data;

[0007] According to the inclination true curvature of the local curved surface corresponding to each of the common center points, determine a first curvature combination, a second curvature combination, and a common center point corresponding to the first curvature combination and a common center point corresponding to the second curvature combination, wherein the first curvature combination includes at least two positive inclination true curvature extreme value points, and the second curvature combination includes at least two negative inclination true curvature extreme value points;

[0008] The fault plane distribution data in the target reservoir is determined according to the common center point corresponding to the first curvature combination and the common center point corresponding to the second curvature combination.

[0009] According to another aspect of the present invention, a fault plane analysis device based on dip true curvature is provided, the device comprising:

[0010] A depth domain layer acquisition module is used to obtain seismic layer data of the target reservoir in the depth domain;

[0011] A curvature calculation module, used to determine the true curvature of the local surface corresponding to each common center point in the seismic layer data;

[0012] A curvature combination module, for determining a first curvature combination, a second curvature combination, and a common center point corresponding to the first curvature combination and a common center point corresponding to the second curvature combination according to the inclination true curvature of the local surface corresponding to each of the common center points, wherein the first curvature combination includes at least two positive inclination true curvature extreme value points, and the second curvature combination includes at least two negative inclination true curvature extreme value points;

[0013] A fault plane analysis module is used to determine the fault plane distribution data in the target reservoir according to the common center point corresponding to the first curvature combination and the common center point corresponding to the second curvature combination.

[0014] According to another aspect of the present invention, an electronic device is provided, the electronic device comprising:

[0015] at least one processor; and

[0016] a memory communicatively connected to the at least one processor; wherein,

[0017] The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the fault plane analysis method based on dip true curvature described in any embodiment of the present invention.

[0018] According to another aspect of the present invention, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the fault plane analysis method based on the dip true curvature described in any embodiment of the present invention when executed.

[0019] The technical solution of the fault plane analysis method based on dip true curvature provided by the embodiment of the present invention obtains the seismic horizon data of the target reservoir in the depth domain; determines the dip true curvature of the local surface corresponding to each common center point in the seismic horizon data; determines the first curvature combination, the second curvature combination, and the common center point corresponding to the first curvature combination and the common center point corresponding to the second curvature combination according to the dip true curvature of the local surface corresponding to each common center point, wherein the first curvature combination includes at least two positive dip true curvature extreme points, and the second curvature combination includes at least two negative dip true curvature extreme points; determines the fault plane distribution data in the target reservoir according to the common center point corresponding to the first curvature combination and the common center point corresponding to the second curvature combination. The technical effect of determining the fault plane distribution data in the target reservoir according to the dip true curvature corresponding to each common center point in the seismic horizon data is achieved, which is simple and accurate.

[0020] It should be understood that the contents described in this section are not intended to identify the key or important features of the embodiments of the present invention, nor are they intended to limit the scope of the present invention. Other features of the present invention will become easily understood through the following description. BRIEF DESCRIPTION OF THE DRAWINGS

[0021] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0022] Figure 1 is a flow chart of a fault plane analysis method based on dip true curvature provided in accordance with an embodiment of the present invention;

[0023] Figure 2 is a schematic structural diagram of a target reservoir provided according to an embodiment of the present invention;

[0024] Figure 3 is a schematic diagram of a facet provided according to an embodiment of the present invention;

[0025] Figure 4 is a schematic diagram of a local curved surface provided according to an embodiment of the present invention;

[0026] Figure 5 is a schematic diagram of the inclination and direction of the local curved surface corresponding to the common center point B at the common center point B provided by an embodiment of the present invention;

[0027] Figure 6 is a schematic diagram of a trend line and a tendency line provided according to an embodiment of the present invention;

[0028] Figure 7 is a schematic diagram of the true curvature of the target reservoir provided according to an embodiment of the present invention;

[0029] Figure 8 is a partially enlarged schematic diagram of the inclination true curvature provided according to an embodiment of the present invention;

[0030] Fig. 9 is a schematic diagram of the distribution of fault boundaries corresponding to the first curvature combination and the fault boundaries corresponding to the second curvature combination provided according to an embodiment of the present invention;

[0031] Fig.10 is a schematic diagram of fault centerline segment distribution provided according to an embodiment of the present invention;

[0032] Fig.11 is a three-dimensional schematic diagram of a target reservoir including multiple fault planes provided according to an embodiment of the present invention;

[0033] Fig.12 is a structural schematic diagram of a fault plane analysis device based on dip true curvature provided according to an embodiment of the present invention;

[0034] Fig.13 It is a schematic diagram of the structure of an electronic device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0035] In order to enable those skilled in the art to better understand the scheme of the present invention, the technical scheme in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present invention.

[0036] It should be noted that the terms "purpose", "goal", etc. in the specification and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units that are clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0037] Figure 1 The present invention provides a flowchart of a fault plane analysis method based on dip true curvature. This embodiment is applicable to the case where fault plane distribution data is automatically determined based on seismic layer data in the depth domain. The method can be performed by a fault plane analysis device based on dip true curvature. The fault plane analysis device based on dip true curvature can be implemented in the form of hardware and / or software. The fault plane analysis device based on dip true curvature can be configured in a processor of an electronic device. Figure 1 As shown, the method includes:

[0038] S110, obtaining seismic layer data of the target reservoir in the depth domain.

[0039] Figure 2 The schematic diagram of the structure of the target reservoir provided by the embodiment of the present invention is shown in FIG.

[0040] The seismic horizon data of the target reservoir in the depth domain can be determined based on the seismic horizon data of the target reservoir in the time domain. In one embodiment, the seismic horizon data of the target reservoir in the time domain is obtained; based on the seismic horizon data of the target reservoir in the time domain and the longitudinal wave velocity of the target reservoir, the seismic horizon data of the target reservoir in the depth domain is determined. It can be understood that the seismic horizon data of the target reservoir in the depth domain is three-dimensional data, and each coordinate point (common center point) includes a row number, a channel number, and a horizon depth data. Therefore, the seismic horizon data of the target reservoir in the depth domain is mapped to a three-dimensional coordinate system to obtain a three-dimensional model of the target reservoir.

[0041] Specifically, the depth domain horizon of the target reservoir is determined by the following formula:

[0042]

[0043] Among them, time is the time horizon, and the one-way propagation time of seismic waves is v is the P-wave velocity. According to the above formula, the depth domain layer z of the target reservoir is determined.

[0044] S120, determining the true curvature of the local surface corresponding to each common center point in the seismic layer data.

[0045] In one embodiment, the following steps are used to determine the true curvature of the local surface corresponding to each common center point in the seismic horizon data, including:

[0046] Step a1, determining the local surface expression corresponding to each common center point in the seismic layer data.

[0047] In one embodiment, a two-variable binomial equation is used to construct the local surface expression corresponding to each common center point. Specifically, Figure 2 As shown, for each common center point, such as common center point A, common center point B and common center point C; determine the surface element where each common center point is located, and the surface element includes 5×5 common center points including the corresponding common center point, and the corresponding common center point is in the middle of the 5×5 common center points (see Figure 3 ), input the data of the 25 common center points into the binary binomial equation, and obtain the specific values ​​of the polynomial coefficients a, b, c, d, e, and f in the equation, so as to obtain the local surface expression corresponding to the current common center point, that is, the local surface corresponding to the current common center point. Figure 4 The local surfaces corresponding to the common center point A, the common center point B and the common center point C are shown respectively.

[0048] The two-variable binomial equation is as follows:

[0049] z(x,y)=ax 2 +by 2+cxy+dx+ey+f

[0050] Step a2: determining the direction and inclination of the local surface corresponding to each of the common center points based on the local surface expressions corresponding to each of the common center points.

[0051] For each common center point, after the local surface corresponding to the common center point is determined, the gradient vector of the local surface corresponding to the common center point at the common center point is determined, and the modulus of the gradient is used as the formation gradient corresponding to the common center point. The projection of the formation gradient vector corresponding to the common center point on the horizontal plane is used as the dip corresponding to the common center point; the direction perpendicular to the formation dip in the horizontal plane is used as the formation strike. Among them, the strike can be used to determine the extension direction of the fault plane. Specifically:

[0052] The gradient vector of the local surface corresponding to each common center point at the corresponding common center point is as follows:

[0053]

[0054] in are the first-order partial derivatives of the local surface with respect to x and y, respectively, which can be expressed as:

[0055]

[0056]

[0057] The modulus of the gradient vector is defined as the formation gradient grad, which can be expressed as:

[0058]

[0059] The direction of the gradient vector is taken as the formation dip φ dip , which can be expressed as:

[0060]

[0061] Where atan2 is the four-quadrant inverse tangent function.

[0062] The direction perpendicular to the formation dip in the horizontal plane is taken as the formation strike, which can be expressed as:

[0063]

[0064] Figure 5 A schematic diagram of the inclination and direction of the local surface corresponding to the common center point provided by the embodiment of the present invention at the common center point. Figure 5Take the local surface corresponding to the common center point B as an example. Assuming that the formation gradient of the local surface corresponding to the common center point B is equal to 2 (that is, the directional derivative is the maximum value and equal to 2) and the azimuth is 12° (that is, α = 12°), the direction of the formation gradient (which can be understood as the maximum slope direction of the local surface corresponding to the common center point B) projected on the horizontal plane (x, y) is determined as the inclination of the local surface corresponding to the common center point B (such as Figure 5 The s in the equation is the inclination), and the direction perpendicular to the inclination of the local surface is determined as the strike direction of the local surface (e.g. Figure 5 The t in the figure is the strike direction, which can be used to determine the extension direction of the fault plane).

[0065] refer to Figure 5 , still taking the local surface corresponding to the common center point B as an example. The inclination s corresponding to the local surface maps the vector s′ on the two-dimensional coordinate plane, and the angle between the vector s′ and the x-axis is According to the direction t of the local surface, determine the vector t' of each local surface's direction mapped on the horizontal plane in the three-dimensional coordinate system, and the angle between the vector t' and the x-axis Among them, the existing calculation formula can be used to solve the angle of the local surface, that is, Among them, atan2 represents the four-quadrant inverse tangent, d and e represent the coefficients of the polynomial of each local surface. For details, please refer to the above-mentioned part on determining the local surface expression.

[0066] According to the above method, we can find is 12°, It is 102°.

[0067] Furthermore, for the local surfaces corresponding to each common center point, after determining the direction and inclination of the local surface corresponding to the current common center point; determine the target plane based on the direction and the inclination; use the target plane to cut the local surface corresponding to the current common center point to obtain the cross-section corresponding to the common center point; and use the curvature of the cross-section corresponding to each common center point as the true curvature of the inclination of the corresponding local surface at the corresponding common center point.

[0068] Figure 5 The solid line PQ in is the cross section of the local surface with the common center point B. Figure 6 The dotted line s1s2 (inclination line s1s2) in is the cross-section line of the local surface corresponding to the common center point A at the common center point A. The method for determining the cross-section line of the local surface corresponding to the common center point C is similar and will not be repeated here.

[0069] Compared with solving the curvature of the entire surface, decomposing the surface into multiple local surfaces and calculating the curvature of the cross section of each local surface, that is, the inclination true curvature of the local surface, can improve the accuracy of determining the curvature corresponding to the common center point. The specific method for solving the cross section of each local surface is as follows:

[0070] By rotating the x and y plane coordinates into the coordinate system represented by t and s, the local surface function degenerates into a parabolic function. The specific formula is derived as follows:

[0071] x=scosφ dip ,y=ssinφ dip .

[0072]

[0073]

[0074] Therefore, it can be seen from the above function that the local surface function is transformed into a two-variable linear equation, which degenerates into a parabolic function. For the specific determination process of a, b, c, d, e, and f, refer to S120.

[0075] After finding the function of the cross section of each local surface, the curvature is calculated using the existing technology, and the true curvature of the cross section of the local surface is further solved in combination with the above formula, as follows:

[0076]

[0077] The above method can solve the curvature of the cross-section of each local surface at the corresponding common center point, and compared with the prior art that directly solves the curvature of the surface, the embodiment of the present invention converts it into solving the curvature of the cross-section of the local surface, which can improve the accuracy of determining the true curvature of the surface, and lays the foundation for the subsequent accurate determination of the boundary of the fault plane.

[0078] S130. Determine a first curvature combination, a second curvature combination, and a common center point corresponding to the first curvature combination and a common center point corresponding to the second curvature combination according to the inclination true curvature of the local surface corresponding to each of the common center points, wherein the first curvature combination includes at least two extreme points of true curvature with a positive inclination, and the second curvature combination includes at least two extreme points of true curvature with a negative inclination.

[0079] The first curvature combination corresponds to the second curvature combination. The first curvature combination includes the positively inclined true curvature extreme value points in the curvature pair, and the second curvature combination includes the negatively inclined true curvature extreme value points in the curvature pair.

[0080] In one embodiment, the true curvature of the inclination of the target reservoir is determined according to the true curvature of the local surface corresponding to each common center point, such as Figure 7 As shown in the figure, the darker the color, the higher the absolute value of the true curvature of the local surface corresponding to the position.

[0081] In one embodiment, the curvature image is determined according to the true curvature of the local surface corresponding to each of the common center points at the corresponding common center points (see Figure 8 ). Determine all the curvature pairs included in the curvature image, and add the positively inclined true curvature extreme points in all the curvature pairs to the first curvature combination (corresponding to Figure 8 The black area on the left side of B in Figure 2 adds the negatively inclined true curvature extreme points in all curvature pairs to the second curvature combination (corresponding to Figure 8 The curvature image can be understood as a projection image obtained by projecting the true curvature of the inclination corresponding to each common center point onto a plane coordinate system.

[0082] S140. Determine fault plane distribution data in the target reservoir according to the common center point corresponding to the first curvature combination and the common center point corresponding to the second curvature combination.

[0083] The fault hanging plate boundary distribution data is determined according to the common center point corresponding to the first curvature combination; the fault footwall boundary distribution data is determined according to the common center point corresponding to the second curvature combination; the fault hanging plate boundary corresponding to the fault hanging plate boundary distribution data is opposite to the fault footwall boundary corresponding to the fault footwall boundary distribution data.

[0084] Specifically, assuming that the true curvature of the inclination corresponding to the common center point A is 0.247, multiple three-dimensional coordinate points can be determined along the direction t of the common center point A, and their corresponding true curvature of the inclination is around 0.247; the true curvature of the inclination corresponding to the common center point B is -0.038, which is approximately 0 (its cross section is approximately a straight line). Similarly, multiple three-dimensional coordinate points can be determined along the direction t of the common center point B, and their corresponding true curvature of the inclination is approximately 0; the true curvature of the inclination corresponding to the common center point C is -0.232, and multiple three-dimensional coordinate points can be determined along the direction t of the common center point C, and their corresponding true curvature of the inclination is approximately -0.232. If 0.247 is the positive dip true curvature extreme point and -0.232 is the negative dip true curvature extreme point, then the fault hanging wall boundary distribution data is determined based on a plurality of three-dimensional coordinate points corresponding to the positive dip true curvature extreme point, and the fault footwall boundary distribution data is determined based on a plurality of three-dimensional coordinate points corresponding to the negative dip true curvature extreme point. It can be understood that the fault hanging wall boundary corresponding to the fault hanging wall boundary distribution data is opposite to the fault footwall boundary corresponding to the fault footwall boundary distribution data, and both are distributed on both sides of the fault centerline segment.

[0085] In one embodiment, a common center point combination corresponding to the fault center line segment is determined based on the common center point corresponding to the first curvature combination and the common center point corresponding to the second curvature combination; and a fault center line segment is determined based on the common center point combination corresponding to the fault center line segment. It can be understood that the common center point combination corresponding to the fault center line segment is located at Fig. 9 The middle of the black solid line A and the black solid line C in the figure, and including the common center point B.

[0086] Fig.10 A schematic diagram of the distribution of fault centerline segments provided in an embodiment of the present invention is specifically as follows: Fig.11 A top view of the fault plane in Figure 1. Different surface colors indicate different fault levels.

[0087] In one embodiment, the fault centerline segment corresponding to the fault centerline segment distribution data, that is, the fault centerline segment distribution schematic diagram is displayed in a visual interactive interface, so that the user can intuitively determine the distribution of the fault plane.

[0088] In one embodiment, a two-dimensional top view of the fault plane and / or a three-dimensional stereogram of the fault plane corresponding to the fault plane distribution data are displayed in a visual interactive interface.

[0089] The technical solution of the fault plane analysis method based on dip true curvature provided by the embodiment of the present invention obtains the seismic horizon data of the target reservoir in the depth domain; determines the dip true curvature of the local surface corresponding to each common center point in the seismic horizon data; determines the first curvature combination, the second curvature combination, and the common center point corresponding to the first curvature combination and the common center point corresponding to the second curvature combination according to the dip true curvature of the local surface corresponding to each common center point, wherein the first curvature combination includes at least two positive dip true curvature extreme points, and the second curvature combination includes at least two negative dip true curvature extreme points; determines the fault plane distribution data in the target reservoir according to the common center point corresponding to the first curvature combination and the common center point corresponding to the second curvature combination. The technical effect of determining the fault plane distribution data in the target reservoir according to the dip true curvature corresponding to each common center point in the seismic horizon data is achieved, which is simple and accurate.

[0090] Fig.12 The schematic diagram of the structure of the fault plane analysis device based on the true curvature of the dip provided by the embodiment of the present invention. Fig.12 As shown, the device comprises:

[0091] A depth domain layer acquisition module 110 is used to acquire seismic layer data of a target reservoir in the depth domain;

[0092] The curvature calculation module 120 is used to determine the true curvature of the local surface corresponding to each common center point in the seismic layer data;

[0093] The curvature combination module 130 is used to determine a first curvature combination, a second curvature combination, and a common center point corresponding to the first curvature combination and a common center point corresponding to the second curvature combination according to the inclination true curvature of the local surface corresponding to each of the common center points, wherein the first curvature combination includes at least two positive inclination true curvature extreme value points, and the second curvature combination includes at least two negative inclination true curvature extreme value points;

[0094] The fault plane analysis module 140 is used to determine the fault plane distribution data in the target reservoir according to the common center point corresponding to the first curvature combination and the common center point corresponding to the second curvature combination.

[0095] In one embodiment, the curvature calculation module 120 includes:

[0096] An expression unit, used to determine a local surface expression corresponding to each common center point in the seismic layer data;

[0097] A tendency unit, used for determining the direction and tendency of the local surface corresponding to each of the common center points based on the local surface expression corresponding to each of the common center points;

[0098] The curvature unit is used to determine the true curvature of the inclination of the corresponding local surface at the corresponding common center point based on the direction and inclination of the local surface corresponding to each of the common center points.

[0099] In one embodiment, the curvature unit is used to:

[0100] For the local surfaces corresponding to each common center point, determine the direction and inclination of the local surface corresponding to the current common center point;

[0101] determining a target plane based on the direction and the inclination;

[0102] Using the target plane to cut the local surface corresponding to the current common center point to obtain a cross-section corresponding to the common center point;

[0103] The curvature of the cross-section line corresponding to each of the common center points is taken as the inclination true curvature of the corresponding local surface at the corresponding common center point.

[0104] In one embodiment, the local surface expression is determined based on a two-variable quadratic polynomial equation.

[0105] In one embodiment, the fault plane analysis module 140 is used to:

[0106] Determine the fault hanging wall boundary distribution data according to the common center point corresponding to the first curvature combination;

[0107] Determine the fault footwall boundary distribution data according to the common center point corresponding to the second curvature combination;

[0108] The fault hanging wall boundary corresponding to the fault hanging wall boundary distribution data is opposite to the fault footing wall boundary corresponding to the fault footing wall boundary distribution data.

[0109] In one embodiment, the fault plane distribution data is determined based on the fault hanging wall boundary distribution data and the fault hanging wall boundary distribution data.

[0110] In one embodiment, the device further comprises a display module, wherein the display module is used to:

[0111] The two-dimensional top view of the fault plane and / or the three-dimensional stereogram of the fault plane corresponding to the fault plane distribution data are displayed in the visual interactive interface.

[0112] In one embodiment, the depth domain layer acquisition module 110 is used to:

[0113] Obtain seismic horizon data of the target reservoir in the time domain;

[0114] The seismic layer data of the target reservoir in the depth domain are determined according to the seismic layer data of the target reservoir in the time domain and the P-wave velocity curve in the seismic layer data.

[0115] The technical solution of the fault plane analysis device based on dip true curvature provided by the embodiment of the present invention obtains the seismic horizon data of the target reservoir in the depth domain; determines the dip true curvature of the local surface corresponding to each common center point in the seismic horizon data; determines the first curvature combination, the second curvature combination, and the common center point corresponding to the first curvature combination and the common center point corresponding to the second curvature combination according to the dip true curvature of the local surface corresponding to each common center point, wherein the first curvature combination includes at least two positive dip true curvature extreme points, and the second curvature combination includes at least two negative dip true curvature extreme points; determines the fault plane distribution data in the target reservoir according to the common center point corresponding to the first curvature combination and the common center point corresponding to the second curvature combination. The technical effect of determining the fault plane distribution data in the target reservoir according to the dip true curvature corresponding to each common center point in the seismic horizon data is achieved, which is simple and accurate.

[0116] The fault plane analysis device based on the inclination true curvature provided in the embodiment of the present invention can execute the fault plane analysis method based on the inclination true curvature provided in any embodiment of the present invention, and has the corresponding functional modules and beneficial effects of the execution method.

[0117] Fig.13A schematic diagram of the structure of an electronic device 10 that can be used to implement an embodiment of the present invention is shown. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as personal digital processing, cellular phones, smart phones, wearable devices (such as helmets, glasses, watches, etc.) and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the present invention described and / or required herein.

[0118] like Fig.13 As shown, the electronic device 10 includes at least one processor 11, and a memory connected to the at least one processor 11, such as a read-only memory (ROM) 12, a random access memory (RAM) 13, etc., wherein the memory stores a computer program that can be executed by at least one processor, and the processor 11 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 12 or the computer program loaded from the storage unit 18 to the random access memory (RAM) 13. In the RAM 13, various programs and data required for the operation of the electronic device 10 can also be stored. The processor 11, the ROM 12, and the RAM 13 are connected to each other through a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

[0119] A number of components in the electronic device 10 are connected to the I / O interface 15, including: an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, speakers, etc.; a storage unit 18, such as a disk, an optical disk, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices through a computer network such as the Internet and / or various telecommunication networks.

[0120] The processor 11 may be a variety of general and / or special processing components with processing and computing capabilities. Some examples of the processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any appropriate processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as a fault plane analysis method based on the true curvature of the inclination.

[0121] In some embodiments, the method for parsing a tomographic plane based on the inclination true curvature may be implemented as a computer program, which is tangibly contained in a computer-readable storage medium, such as a storage unit 18. In some embodiments, part or all of the computer program may be loaded and / or installed on the electronic device 10 via the ROM 12 and / or the communication unit 19. When the computer program is loaded into the RAM 13 and executed by the processor 11, one or more steps of determining the tomographic plane based on the inclination true curvature described above may be performed. Alternatively, in other embodiments, the processor 11 may be configured to perform the method for parsing a tomographic plane based on the inclination true curvature in any other appropriate manner (e.g., by means of firmware).

[0122] Various implementations of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), systems on chips (SOCs), load programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various implementations can include: being implemented in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which can be a special purpose or general purpose programmable processor that can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit data and instructions to the storage system, the at least one input device, and the at least one output device.

[0123] Computer programs for implementing the methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, so that when the computer program is executed by the processor, the functions / operations specified in the flow chart and / or block diagram are implemented. The computer program may be executed entirely on the machine, partially on the machine, partially on the machine and partially on a remote machine as a stand-alone software package, or entirely on a remote machine or server.

[0124] In the context of the present invention, a computer-readable storage medium may be a tangible medium that may contain or store a computer program for use by or in combination with an instruction execution system, device or equipment. A computer-readable storage medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, devices or equipment, or any suitable combination of the foregoing. Alternatively, a computer-readable storage medium may be a machine-readable signal medium. A more specific example of a machine-readable storage medium may include an electrical connection based on one or more lines, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0125] To provide interaction with a user, the systems and techniques described herein may be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and a pointing device (e.g., a mouse or trackball) through which the user can provide input to the electronic device. Other types of devices may also be used to provide interaction with the user; for example, the feedback provided to the user may be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user may be received in any form (including acoustic input, voice input, or tactile input).

[0126] The systems and techniques described herein may be implemented in a computing system that includes backend components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes frontend components (e.g., a user computer with a graphical user interface or a web browser through which a user can interact with implementations of the systems and techniques described herein), or a computing system that includes any combination of such backend components, middleware components, or frontend components. The components of the system may be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: a local area network (LAN), a wide area network (WAN), a blockchain network, and the Internet.

[0127] A computing system may include a client and a server. The client and the server are generally remote from each other and usually interact through a communication network. The client and server relationship is generated by computer programs running on the corresponding computers and having a client-server relationship with each other. The server may be a cloud server, also known as a cloud computing server or cloud host, which is a host product in the cloud computing service system to solve the defects of difficult management and weak business scalability in traditional physical hosts and VPS services.

[0128] It should be understood that the various forms of processes shown above can be used to reorder, add or delete steps. For example, the steps described in the present invention can be executed in parallel, sequentially or in different orders, as long as the desired results of the technical solution of the present invention can be achieved, and this document does not limit this.

[0129] The above specific implementations do not constitute a limitation on the protection scope of the present invention. It should be understood by those skilled in the art that various modifications, combinations, sub-combinations and substitutions can be made according to design requirements and other factors. Any modification, equivalent substitution and improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims

1. A fault plane analysis method based on dip true curvature, characterized in that: include: Obtain seismic layer data of the target reservoir in the depth domain; Determine the true curvature of the local surface corresponding to each common center point in the seismic horizon data; According to the inclination true curvature of the local curved surface corresponding to each of the common center points, determine a first curvature combination, a second curvature combination, and a common center point corresponding to the first curvature combination and a common center point corresponding to the second curvature combination, wherein the first curvature combination includes at least two positive inclination true curvature extreme value points, and the second curvature combination includes at least two negative inclination true curvature extreme value points; The fault plane distribution data in the target reservoir is determined according to the common center point corresponding to the first curvature combination and the common center point corresponding to the second curvature combination.

2. The method according to claim 1, characterized in that The determining of the true curvature of the local surface corresponding to each common center point in the seismic horizon data comprises: Determine the local surface expression corresponding to each common center point in the seismic layer data; Determine the direction and inclination of the local surface corresponding to each of the common center points based on the local surface expressions corresponding to each of the common center points; Based on the direction and inclination of the local curved surface corresponding to each of the common center points, the inclination true curvature of the corresponding local curved surface at the corresponding common center point is determined.

3. The method according to claim 2, characterized in that The determining the inclination true curvature of the corresponding local curved surface at the corresponding common center point based on the direction and inclination of the local curved surface corresponding to each of the common center points comprises: For the local surfaces corresponding to each common center point, determine the direction and inclination of the local surface corresponding to the current common center point; determining a target plane based on the direction and the inclination; Using the target plane to cut the local surface corresponding to the current common center point to obtain a cross-section corresponding to the common center point; The curvature of the cross-section line corresponding to each of the common center points is taken as the inclination true curvature of the corresponding local surface at the corresponding common center point.

4. The method according to claim 2, characterized in that: The local surface expression is determined based on a two-variable quadratic polynomial equation.

5. The method according to claim 1, characterized in that Determining fault plane distribution data in the target reservoir according to the common center point corresponding to the first curvature combination and the common center point corresponding to the second curvature combination includes: Determine the fault hanging wall boundary distribution data according to the common center point corresponding to the first curvature combination; Determine the fault footwall boundary distribution data according to the common center point corresponding to the second curvature combination; The fault hanging wall boundary corresponding to the fault hanging wall boundary distribution data is opposite to the fault footing wall boundary corresponding to the fault footing wall boundary distribution data.

6. The method according to claim 5, characterized in that Also includes: The fault plane distribution data is determined according to the fault hanging wall boundary distribution data and the fault footwall boundary distribution data.

7. The method according to claim 6, characterized in that Also includes: The two-dimensional top view of the fault plane and / or the three-dimensional stereogram of the fault plane corresponding to the fault plane distribution data are displayed in the visual interactive interface.

8. A fault plane analysis device based on dip true curvature, characterized in that: include: A depth domain layer acquisition module is used to obtain seismic layer data of the target reservoir in the depth domain; A curvature calculation module, used to determine the true curvature of the local surface corresponding to each common center point in the seismic layer data; A curvature combination module, for determining a first curvature combination, a second curvature combination, and a common center point corresponding to the first curvature combination and a common center point corresponding to the second curvature combination according to the inclination true curvature of the local surface corresponding to each of the common center points, wherein the first curvature combination includes at least two positive inclination true curvature extreme value points, and the second curvature combination includes at least two negative inclination true curvature extreme value points; A fault plane analysis module is used to determine the fault plane distribution data in the target reservoir according to the common center point corresponding to the first curvature combination and the common center point corresponding to the second curvature combination.

9. An electronic device, characterized in that: The electronic device comprises: at least one processor; and a memory communicatively connected to the at least one processor; wherein, The memory stores a computer program executable by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can perform the fault plane analysis method based on dip true curvature according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a processor to implement the fault plane analysis method based on dip true curvature according to any one of claims 1 to 7 when executed.

Images