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

Abstract

The application discloses a fault plane analysis method and device based on a tendency true curvature, equipment and a medium, and the method comprises the following steps: acquiring seismic horizon data of a target reservoir in a depth domain; determining the tendency true curvature of a local surface corresponding to each common midpoint in the seismic horizon data; determining a first curvature combination, a second curvature combination, and common midpoints corresponding to the first curvature combination and common midpoints corresponding to the second curvature combination according to the tendency true curvature of the local surface corresponding to each common midpoint, wherein the first curvature combination comprises at least two positive tendency true curvature extreme points, and 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 midpoints corresponding to the first curvature combination and the common midpoints corresponding to the second curvature combination. The technical scheme provided by the embodiment of the application can automatically and accurately determine the fault plane distribution data.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of oil exploration, and in particular to a fault plane analysis method and device based on true dip curvature, equipment and a medium. BACKGROUND

[0002] In oil exploration, the important problem affecting the exploration and development of oil and gas fields is fault structure. The existing technology often uses conventional three-dimensional seismic interpretation. Conventional three-dimensional seismic interpretation is basically carried out on a vertical seismic profile, which requires each survey line in the work area to be interpreted one by one to determine the planar distribution of fault and stratum characteristics, so as to analyze, judge and interpret on the plane by using artificial methods. However, this method is complicated and requires manual analysis, and cannot automatically and accurately determine the fault plane distribution data. SUMMARY

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

[0004] According to an aspect of the present application, a fault plane analysis method based on true dip curvature is provided, which comprises:

[0005] obtaining seismic horizon data of a target reservoir in a depth domain;

[0006] determining the true dip curvature of the local surface corresponding to each common midpoint in the seismic horizon data;

[0007] determining a first curvature combination, a second curvature combination, and the common midpoint corresponding to the first curvature combination and the common midpoint corresponding to the second curvature combination according to the true dip curvature of the local surface corresponding to each common midpoint, the first curvature combination comprising at least two positive true dip curvature extreme points, and the second curvature combination comprising at least two negative true dip curvature extreme points;

[0008] determining the fault plane distribution data in the target reservoir according to the common midpoint corresponding to the first curvature combination and the common midpoint corresponding to the second curvature combination.

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

[0010] a depth domain horizon acquisition module for obtaining seismic horizon data of a target reservoir in a depth domain;

[0011] a curvature calculation module for determining the true dip curvature of the local surface corresponding to each common midpoint in the seismic horizon data;

[0012] a curvature combination module configured to determine a first curvature combination, a second curvature combination, and a common point corresponding to the first curvature combination and a common point corresponding to the second curvature combination according to the local surface tendency true curvatures of each of the common points, the first curvature combination comprising at least two positive tendency true curvature extreme points, and the second curvature combination comprising at least two negative tendency true curvature extreme points;

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

[0014] According to another aspect of the present application, an electronic device is provided, which comprises:

[0015] at least one processor; and

[0016] a memory connected to the at least one processor in communication; 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 to enable the at least one processor to perform the tendency true curvature based fault plane analysis method according to any one of the embodiments of the present application.

[0018] According to another aspect of the present application, a computer readable storage medium is provided, which stores computer instructions for enabling a processor to perform the tendency true curvature based fault plane analysis method according to any one of the embodiments of the present application.

[0019] The tendency true curvature based fault plane analysis method provided by the embodiments of the present application comprises the following steps: obtaining seismic horizon data of a target reservoir in a depth domain; determining tendency true curvatures of local surfaces corresponding to each common point in the seismic horizon data; determining a first curvature combination, a second curvature combination, and a common point corresponding to the first curvature combination and a common point corresponding to the second curvature combination according to the tendency true curvatures of the local surfaces corresponding to each of the common points, the first curvature combination comprising at least two positive tendency true curvature extreme points, and the second curvature combination comprising at least two negative tendency true curvature extreme points; and determining fault plane distribution data in the target reservoir according to the common point corresponding to the first curvature combination and the common point corresponding to the second curvature combination. The technical effect of determining the fault plane distribution data in the target reservoir according to the tendency true curvatures of the local surfaces corresponding to each of the common points in the seismic horizon data is achieved, and the method is simple and accurate.

[0020] It should be understood that the matters described in this detailed description are intended to be illustrative only and are not intended to limit the scope of the present application. Other features of the present application will be apparent from the following description. BRIEF DESCRIPTION OF DRAWINGS

[0021] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings needed to be used in the embodiments will be briefly introduced as follows. Obviously, the drawings in the following description are only some embodiments of the present application, and for those skilled in the art, other drawings can also be obtained from these drawings without creative effort.

[0022] Figure 1 is a flow chart of a fault plane analysis method based on true dip curvature according to an embodiment of the present application;

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

[0024] Figure 3 is a bin schematic diagram according to an embodiment of the present application;

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

[0026] Figure 5 is a schematic diagram of a dip and strike of a local surface corresponding to a common depth point B at the common depth point B according to an embodiment of the present application;

[0027] Figure 6 is a strike line and dip line schematic diagram according to an embodiment of the present application;

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

[0029] Figure 8 is a local zoom-in schematic diagram of true dip curvature according to an embodiment of the present application;

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

[0031] Figure 10 is a fault center line segment distribution schematic diagram according to an embodiment of the present application;

[0032] Figure 11 is a three-dimensional schematic diagram of a target reservoir including a plurality of fault planes according to an embodiment of the present application;

[0033] Figure 12 is a structural schematic diagram of a fault plane analysis device based on true curvature of tendency provided by an embodiment of the present application;

[0034] Figure 13 is a structural schematic diagram of an electronic device provided by an embodiment of the present application. DETAILED DESCRIPTION

[0035] In order to make the personnel in the technical field better understand the present application scheme, the technical scheme in the embodiment of the present application will be described clearly and completely in combination with the drawings in the embodiment of the present application. Obviously, the described embodiment is only a part of the embodiment of the present application, not all. Based on the embodiment in the present application, all other embodiments obtained by the person skilled in the art without creative labor should belong to the protection scope of the present application.

[0036] It should be noted that the terms "purpose", "target" and the like in the specification and claims of the present application and the above-mentioned drawings are used to distinguish similar objects, and do not necessarily describe a specific order or sequence. It should be understood that the data thus used can be interchanged under appropriate circumstances, so that the embodiments of the present application described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "include" and "have" and any variations thereof are intended to cover non-exclusive inclusion, for example, a process, method, system, product or device including a series of steps or units does not necessarily limit to those steps or units clearly listed, but can include other steps or units not clearly listed or inherent to these processes, methods, products or devices.

[0037] Figure 1 The flowchart of the fault plane analysis method based on true curvature of tendency provided by the embodiment of the present application can be applied to the case of automatically determining fault plane distribution data based on seismic horizon data in depth domain. The method can be executed by a fault plane analysis device based on true curvature of tendency, which can be realized in the form of hardware and / or software, and can be configured in the processor of an electronic device. As shown in the figure, the method comprises: Figure 1

[0038] S110, acquiring seismic horizon data of a target reservoir in depth domain.

[0039] Figure 2 The structural schematic diagram of the target reservoir provided by the embodiment of the present application is shown in the figure. The target reservoir refers to the target horizon, i.e. the horizon of interest, which includes the fault plane.

[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, and the seismic horizon data of the target reservoir in the depth domain is determined according to the seismic horizon data of the target reservoir in the time domain and the P-wave velocity of the target reservoir. It can be understood that the seismic horizon data of the target reservoir in the depth domain is three-dimensional data, each coordinate point (common midpoint) includes a row number, a trace number, and horizon depth data, and therefore the three-dimensional model of the target reservoir can be obtained by mapping the seismic horizon data of the target reservoir in the depth domain to a three-dimensional coordinate system.

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

[0042]

[0043] wherein time is the time horizon, and the seismic wave one-way propagation time is v is the P-wave velocity. The horizon z in the depth domain of the target reservoir is determined according to the above formula.

[0044] S120, determining the tendency true curvature of the local surface corresponding to each common midpoint in the seismic horizon data.

[0045] In one embodiment, the tendency true curvature of the local surface corresponding to each common midpoint in the seismic horizon data is determined by the following steps, including:

[0046] Step a1, determining the local surface expression corresponding to each common midpoint in the seismic horizon data.

[0047] In one embodiment, the local surface expression corresponding to each common midpoint is constructed by using a binary binomial equation. Specifically, as shown in Figure 2 for each common midpoint, such as common midpoint A, common midpoint B, and common midpoint C, the bin of each common midpoint is determined, which includes 5*5 common midpoints corresponding to the common midpoint, and the corresponding common midpoint is in the middle of the 5*5 common midpoints (see Figure 3 ), the data of the 25 common midpoints is input into the binary binomial equation, and the specific values of the polynomial coefficients a, b, c, d, e, and f in the equation are obtained, so as to obtain the local surface expression corresponding to the current common midpoint, that is, the local surface corresponding to the current common midpoint. Figure 4 The local surfaces corresponding to common midpoint A, common midpoint B, and common midpoint C are shown.

[0048] wherein the binary binomial equation is as follows:

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

[0050] Step a2, determining the strike and dip of the local surface corresponding to each said common midpoint based on the local surface expression corresponding to each said common midpoint.

[0051] For each common midpoint, after determining the local surface corresponding to the common midpoint, determining the gradient vector of the local surface corresponding to the common midpoint at the common midpoint, taking the module length of the gradient as the formation gradient corresponding to the common midpoint. Taking the projection of the formation gradient vector corresponding to the common midpoint on the horizontal plane as the dip corresponding to the common midpoint; taking the direction perpendicular to the formation dip in the horizontal plane as the strike. Wherein, 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 midpoint at the corresponding common midpoint is as follows:

[0053]

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

[0055]

[0056]

[0057] The module length 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] Wherein atan2 is the four quadrant arctangent function.

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

[0063]

[0064] Figure 5 The schematic diagram of the dip and strike of the local surface corresponding to the common midpoint at the common midpoint provided by the embodiment of the present application is shown in FIG. 1. Taking Figure 5Let's take the local surface corresponding to the common center point B as an example. Assume the stratigraphic gradient of the local surface corresponding to the common center point B is equal to 2 (i.e., the directional derivative is at its maximum value and equal to 2) and its azimuth is 12° (i.e., α = 12°). Then, the direction of the stratigraphic gradient (which can be understood as the direction of the maximum slope of the local surface corresponding to the common center point B) projected onto the horizontal plane (x, y) is determined as the dip direction of the local surface corresponding to the common center point B (e.g., ...). Figure 5 In this context, 's' represents the inclination, and the direction perpendicular to the inclination of the local surface is defined as the orientation of the local surface (e.g., ...). Figure 5 In this context, 't' represents the orientation, which can be used to determine the direction of the fault plane's extension.

[0065] refer to Figure 5 Let's take the local surface corresponding to the common center point B as an example for illustration. The inclination s of this local surface maps to the vector s′ on the two-dimensional coordinate plane, and the angle between vector s′ and the x-coordinate axis... Based on the orientation t of the local surface, determine the vector t′ that maps the orientation of each local surface to on the horizontal plane in the three-dimensional coordinate system, and the angle between this vector t′ and the x-coordinate axis. The included angle of the local surface can be calculated using existing formulas, i.e. Where atan2 represents the arctangent of the four quadrants, and d and e represent the coefficients of the polynomial of each local surface, as detailed in the section on determining the expression of the local surface above.

[0066] The above method can be used to calculate... It is 12°. It is 102°.

[0067] Furthermore, for each local surface corresponding to a common center point, after determining the direction and inclination of the local surface corresponding to the current common center point; a target plane is determined based on the direction and inclination; the target plane is used to cut the local surface corresponding to the current common center point to obtain the transverse section corresponding to the common center point; the curvature of the transverse section corresponding to each common center point is used 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 the diagram is the transverse section of the local surface at the common center point B. Figure 6 The dashed lines s1s2 (tendency lines s1s2) are the transverse sections of the local surface corresponding to the common center point A at that point. The method for determining the transverse sections 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 whole surface, the surface is disassembled into multiple local surfaces, and the curvature of the cross section of each local surface is solved, that is, the local surface is inclined to the true curvature, which can improve the accuracy of the curvature determination of the common point. The method for solving the cross section of each local surface is as follows:

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

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

[0072]

[0073]

[0074] Therefore, from the above function, the local surface function is converted into a binary linear equation, which degenerates into a parabolic function. Among them, The specific determination process of a, b, c, d, e and f is described in S120.

[0075] After the function of the cross section of each local surface is solved, the curvature is calculated by using the prior art, and the true curvature of the cross section of the local surface is further solved by combining 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 point, and compared with the prior art of directly solving the curvature of the surface, the curvature of the cross section of the local surface is converted by the embodiment of the application, which can improve the accuracy of the determination of the inclined true curvature of the surface, and lays a foundation for the subsequent accurate determination of the boundary of the fault surface.

[0078] S130, according to the inclined true curvature of the local surface corresponding to each common point, determine the first curvature combination, the second curvature combination, and the common point corresponding to the first curvature combination and the common point corresponding to the second curvature combination, the first curvature combination includes at least two positive inclined true curvature extreme points, and the second curvature combination includes at least two negative inclined true curvature extreme points.

[0079] The first curvature combination and the second curvature combination correspond. The first curvature combination includes the positive inclined true curvature extreme point in the curvature pair, and the second curvature combination includes the negative inclined true curvature extreme point in the curvature pair.

[0080] In one embodiment, according to the inclined true curvature of the local surface corresponding to each common point, the inclined true curvature plan of the target reservoir is determined, as shown inFigure 7 The darker the color in the figure, the higher the absolute value of the tendency true curvature of the local surface corresponding to the position.

[0081] In one embodiment, a curvature image is determined according to the tendency true curvature of each said co-center point corresponding local surface at the corresponding co-center point (see FIG. 2B). Figure 8 The determination of the curvature image includes adding all the positive tendency true curvature extreme points in all the curvature pairs to a first curvature combination (corresponding to the black area on the left side of FIG. 2B), and adding all the negative tendency true curvature extreme points in all the curvature pairs to a second curvature combination (corresponding to the black area on the right side of FIG. 2B). Wherein, the curvature image can be understood as a projection image obtained by projecting the tendency true curvature of each co-center point corresponding local surface to a plane coordinate system. Figure 8 Figure 8 S140, determining the fault plane distribution data in the target reservoir according to the co-center points corresponding to the first curvature combination and the co-center points corresponding to the second curvature combination.

[0082] According to the co-center points corresponding to the first curvature combination, the hanging wall boundary distribution data of the fault is determined; according to the co-center points corresponding to the second curvature combination, the footwall boundary distribution data of the fault is determined; the hanging wall boundary corresponding to the hanging wall boundary distribution data of the fault is opposite to the footwall boundary corresponding to the footwall boundary distribution data of the fault.

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

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

[0085] ​In one embodiment, the combination of the common points corresponding to the fault center line segment is determined according to the combination of the common points corresponding to the first curvature and the combination of the common points corresponding to the second curvature. The fault center line segment is determined according to the combination of the common points corresponding to the fault center line segment. It can be understood that the combination of the common points corresponding to the fault center line segment is located in the middle of the black solid line A and the black solid line C in FIG. 1 and includes the common point B. Figure 9

[0086] Figure 10 The distribution diagram of the fault center line segment provided by the embodiment of the present application is a top view of the fault plane in FIG. 1. Different colors of the curved surface represent different fault levels. Figure 11

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

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

[0089] The technical scheme of the fault plane analysis method based on the preferred true curvature provided by the embodiment of the present application comprises the following steps: obtaining the seismic horizon data of a target reservoir in a depth domain; determining the preferred true curvature of the local curved surface corresponding to each common point in the seismic horizon data; determining a first curvature combination, a second curvature combination, a common point corresponding to the first curvature combination, and a common point corresponding to the second curvature combination according to the preferred true curvature of the local curved surface corresponding to each common point, the first curvature combination including at least two positive preferred true curvature extreme points, and the second curvature combination including at least two negative preferred true curvature extreme points; and determining the distribution data of the fault plane in the target reservoir according to the common point corresponding to the first curvature combination and the common point corresponding to the second curvature combination. The technical effect of determining the distribution data of the fault plane in the target reservoir according to the preferred true curvature of the local curved surface corresponding to each common point in the seismic horizon data is achieved, and the method is simple and accurate.

[0090] Figure 12 The structure diagram of the fault plane analysis device based on the preferred true curvature provided by the embodiment of the present application is shown in FIG. 2. As shown in FIG. 2, the device comprises: Figure 12 A depth domain horizon acquisition module 110 is configured to acquire the seismic horizon data of a target reservoir in a depth domain.

[0091] A curvature calculation module 120 is configured to determine the preferred true curvature of the local curved surface corresponding to each common point in the seismic horizon data.

[0092]

[0093] ​​​The curvature combination module 130 is configured to determine a first curvature combination and a second curvature combination according to the local surface of each of the common points, and the common point corresponding to the first curvature combination and the common point corresponding to the second curvature combination, the first curvature combination including at least two positive inclination true curvature extreme points, and the second curvature combination including at least two negative inclination true curvature extreme points.

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

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

[0096] An expression unit is configured to determine the local surface expression corresponding to each of the common points in the seismic horizon data.

[0097] An inclination unit is configured to determine the trend and inclination of the local surface corresponding to each of the common points based on the local surface expression corresponding to each of the common points.

[0098] A curvature unit is configured to determine the inclination true curvature of the local surface at the corresponding common point based on the trend and inclination of the local surface corresponding to each of the common points.

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

[0100] For the local surface corresponding to each of the common points, determine the trend and inclination of the local surface corresponding to the current common point.

[0101] Determine a target plane based on the trend and the inclination.

[0102] Cut the local surface corresponding to the current common point by using the target plane to obtain a cross-section line corresponding to the common point.

[0103] Determine the curvature of each of the cross-section lines corresponding to the common points as the inclination true curvature of the local surface at the corresponding common point.

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

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

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

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

[0108] The upper wall boundary corresponding to the upper wall boundary distribution data is opposite to the lower wall boundary corresponding to the lower wall boundary distribution data.

[0109] In one embodiment, the fault plane distribution data is determined according to the upper wall boundary distribution data and the lower wall boundary distribution data.

[0110] In one embodiment, the apparatus further comprises a display module configured to:

[0111] Display a two-dimensional plan view and / or a three-dimensional perspective view of the fault plane corresponding to the fault plane distribution data in the visualization interface.

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

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

[0114] Determine seismic horizon data of the target reservoir in the depth domain according to the seismic horizon data of the target reservoir in the time domain and a P-wave velocity curve in the seismic horizon data.

[0115] The technical scheme of the fault plane analysis device based on the inclined true curvature provided by the embodiments of the present application acquires seismic horizon data of the target reservoir in the depth domain, determines the inclined true curvature of each common midpoint corresponding to a local surface in the seismic horizon data, determines a first curvature combination, a second curvature combination, and a common midpoint corresponding to the first curvature combination and a common midpoint corresponding to the second curvature combination according to the inclined true curvature of each common midpoint corresponding to the local surface, the first curvature combination includes at least two positive inclined true curvature extreme points, and the second curvature combination includes at least two negative inclined true curvature extreme points, and determines fault plane distribution data in the target reservoir according to the common midpoint corresponding to the first curvature combination and the common midpoint corresponding to the second curvature combination. The technical effect of determining the fault plane distribution data in the target reservoir according to the inclined true curvature of each common midpoint in the seismic horizon data is achieved, and the method is simple and accurate.

[0116] The fault plane analysis device based on the inclined true curvature provided by the embodiments of the present application can execute the fault plane analysis method based on the inclined true curvature provided by any of the embodiments of the present application, has the corresponding function modules and beneficial effects of executing the method.

[0117] Figure 13A schematic diagram of an electronic device 10 that can be used to implement embodiments 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 processors, cellular phones, smartphones, wearable devices (e.g., helmets, glasses, watches, etc.), and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely illustrative and are not intended to limit the implementation of the invention described and / or claimed herein.

[0118] like Figure 13 As shown, the electronic device 10 includes at least one processor 11 and a memory, such as a read-only memory (ROM) 12 or a random access memory (RAM) 13, communicatively connected to the at least one processor 11. The memory stores computer programs executable by the at least one processor. The processor 11 can perform various appropriate actions and processes based on the computer program stored in the ROM 12 or loaded from storage unit 18 into the RAM 13. The RAM 13 may also store various programs and data required for the operation of the electronic device 10. The processor 11, ROM 12, and RAM 13 are interconnected via a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

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

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

[0121] In some embodiments, the true-dip curvature based fault plane resolution method can be implemented as a computer program tangibly embodied in a computer readable storage medium, e.g., storage unit 18. In some embodiments, parts or all of the computer program can be loaded and / or installed onto electronic device 10 via, e.g., ROM 12 and / or communication unit 19. When the computer program is loaded onto RAM 13 and executed by processor 11, one or more of the steps of the true-dip curvature based fault plane determination described above can be performed. Alternatively, in other embodiments, processor 11 can be configured to perform the true-dip curvature based fault plane resolution method by other means, e.g., via firmware.

[0122] Various implementations of the systems and techniques described above can be realized in digital electronic circuitry, integrated circuitry, a field programmable gate array (FPGA), an application specific integrated circuit (ASIC), a system on a chip (SOC), a programmable logic device (PLD), a computer hardware, firmware, software, and / or combinations thereof. These various implementations can include implementation in one or more computer programs that are executable and / or interpretable on a programmable system including at least one programmable processor, which can be special or general purpose, coupled to receive data and instructions from, and to transmit data and instructions to, a storage system, at least one input device, and at least one output device.

[0123] Computer programs used to implement the processes of the present application can be written in any combination of one or more programming languages. These computer programs can be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the computer program, when executed by the processor of the machine, implements the functions / acts specified in the flowcharts and / or block diagrams. The computer program can be executed entirely on a machine, partially on a machine, partially on a machine as a stand-alone software package, partially on a machine and partially on a remote machine or entirely on a remote machine or server.

[0124] In the context of the present application, a computer-readable storage medium can be a tangible medium that can contain or store a computer program for use by or in connection with an instruction execution system, apparatus, or device. A computer-readable storage medium can include, but is not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any suitable combination of the foregoing. Alternatively, a computer-readable storage medium can be a machine-readable signal medium. More specific examples of a machine-readable storage medium will include one or more lines of a program of instructions in a transitory signal, a portable computer diskette, 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 disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0125] To provide for interaction with a user, the systems and techniques described here can 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 a trackball) by which the user can provide input to the electronic device. Other kinds of devices can be used to provide for interaction with a user as well; for example, feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form, including acoustic, speech, or tactile input.

[0126] The systems and techniques described here can be implemented in a computing system that includes a back end component (e.g., as a data server), or that includes a middleware component (e.g., an application server), or that includes a front end component (e.g., a user computer having a graphical user interface or a Web browser through which a user can interact with an implementation of the systems and techniques described here), or any combination of such back end, middleware, or front end components. The components of the system can 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] The computing system can include clients and servers. A client and server are generally remote from each other and typically interact through a communication network. The relationship of client and server arises by virtue of computer programs running on the respective computers and having a client-server relationship to each other. The server can 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 large management difficulty and weak business scalability in traditional physical host and VPS service.

[0128] It should be understood that the various forms of flow shown above can be used to reorder, add or delete steps. For example, each step described in the present application can be executed in parallel, sequentially or in a different order, as long as the desired results of the technical solutions of the present application can be achieved, which is not limited herein.

[0129] The above detailed description does not constitute a limitation on the scope of protection of the present application. Those skilled in the art should understand that various modifications, combinations, sub-combinations and substitutions can be made according to design requirements and other factors. Any modifications, equivalent replacements and improvements made within the spirit and principles of the present application shall be included in the scope of protection of the present application.

Claims

1. A fault plane analysis method based on dip true curvature, characterized by, The method comprises the following steps: acquiring seismic horizon data of a target reservoir in a depth domain; determining a local surface corresponding to each common midpoint in the seismic horizon data; determining a first curvature combination and a second curvature combination according to the local surface corresponding to each common midpoint, wherein the first curvature combination comprises at least two positive inclination true curvature extreme points, and the second curvature combination comprises at least two negative inclination true curvature extreme points; determining fault surface distribution data in the target reservoir according to the common midpoint corresponding to the first curvature combination and the common midpoint corresponding to the second curvature combination. The method further comprises the following steps: determining a local surface expression corresponding to each common midpoint in the seismic horizon data; determining the trend and inclination of the local surface corresponding to each common midpoint based on the local surface expression corresponding to each common midpoint; determining the inclination true curvature of the local surface at the corresponding common midpoint based on the trend and inclination of the local surface corresponding to each common midpoint. The method further comprises the following steps: determining the trend and inclination of the local surface corresponding to each common midpoint; determining a target plane based on the trend and inclination; cutting the local surface corresponding to the current common midpoint by using the target plane to obtain a cross section corresponding to the common midpoint; taking the curvature of the cross section corresponding to each common midpoint as the inclination true curvature of the local surface at the corresponding common midpoint. The method further comprises the following steps: determining fault hanging wall boundary distribution data according to the common midpoint corresponding to the first curvature combination; determining fault footwall boundary distribution data according to the common midpoint 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 footwall boundary corresponding to the fault footwall boundary distribution data.

2. The method of claim 1, wherein, The local surface expression is determined based on a binary quadratic polynomial equation.

3. The method of claim 1, wherein, The method further comprises the following steps: determining fault surface distribution data according to the fault hanging wall boundary distribution data and the fault footwall boundary distribution data.

4. The method of claim 3, wherein, The method further comprises the following steps: displaying a two-dimensional overhead view and / or a three-dimensional perspective view of the fault surface corresponding to the fault surface distribution data in a visual interactive interface.

5. A fault plane analysis apparatus based on true curvature of tendency, characterized by, The method comprises the following steps: a depth domain horizon acquisition module, configured to acquire seismic horizon data of a target reservoir in a depth domain; a curvature calculation module, configured to determine the inclination true curvature of a local surface corresponding to each common midpoint in the seismic horizon data; The curvature combination module is configured to determine a first curvature combination and a second curvature combination according to the principal curvatures of the local surface corresponding to each of the common points, and the common point corresponding to the first curvature combination and the common point corresponding to the second curvature combination, the first curvature combination including at least two positive principal curvature extreme points, and the second curvature combination including at least two negative principal curvature extreme points. The fault plane analysis module is configured to determine fault plane distribution data in the target reservoir according to the common point corresponding to the first curvature combination and the common point corresponding to the second curvature combination. The determination of the principal curvatures of the local surface corresponding to each of the common points in the seismic horizon data includes: determining a local surface expression corresponding to each of the common points in the seismic horizon data; determining a strike and a dip of the local surface corresponding to each of the common points based on the local surface expression corresponding to each of the common points; determining the principal curvature of the local surface corresponding to each of the common points at the corresponding common point based on the strike and the dip of the local surface corresponding to each of the common points; The determination of the principal curvature of the local surface corresponding to each of the common points at the corresponding common point based on the strike and the dip of the local surface corresponding to each of the common points includes: determining the strike and the dip of the local surface corresponding to the current common point for the local surface corresponding to each of the common points; determining a target plane based on the strike and the dip; cutting the local surface corresponding to the current common point by using the target plane to obtain a cross-section line corresponding to the common point; taking the curvature of the cross-section line corresponding to each of the common points as the principal curvature of the local surface corresponding to each of the common points at the corresponding common point. The determination of the fault plane distribution data in the target reservoir according to the common point corresponding to the first curvature combination and the common point corresponding to the second curvature combination includes: determining fault hanging wall boundary distribution data according to the common point corresponding to the first curvature combination; determining fault footwall boundary distribution data according to the common 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 footwall boundary corresponding to the fault footwall boundary distribution data.

6. An electronic device, comprising: The electronic device includes: at least one processor; and a memory connected with the at least one processor in communication; 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 to enable the at least one processor to execute the method for analyzing a fault plane based on a principal curvature according to any one of claims 1-4.

7. A computer-readable storage medium, characterized in that, The computer readable storage medium stores computer instructions for enabling a processor to execute the method for analyzing a fault plane based on a principal curvature according to any one of claims 1-4 when executed by the processor.

Images