CN114519784A

CN114519784A - Three-dimensional fault modeling method and device

Info

Publication number: CN114519784A
Application number: CN202210150784.0A
Authority: CN
Inventors: 严志刚; 芦辰; 刘畅; 袁金龙; 姜明; 王鹏; 何鸿鹏
Original assignee: General Coal Research Institute Co Ltd
Current assignee: General Coal Research Institute Co Ltd
Priority date: 2022-02-14
Filing date: 2022-02-14
Publication date: 2022-05-20
Anticipated expiration: 2042-02-14
Also published as: CN114519784B

Abstract

The disclosure provides a modeling method and device for a three-dimensional fault, and relates to the technical field of three-dimensional geological modeling. The method comprises the following steps: acquiring the type of a fault and the fault distance corresponding to each marking point on a trace of the fault; determining an upper coiling line and a lower coiling line of the fault according to the trace and the corresponding fault distance of each marking point; determining a fault boundary line according to the upper coil line, the lower coil line and a preset offset distance; and determining a fault plane according to the type of the fault, the upper coil line, the lower coil line and the fault boundary line. Therefore, the upper coil line and the lower coil line of the fault are determined according to the fault distance corresponding to each marking point on the trace of the fault, and the fault plane is formed according to the upper coil line and the lower coil line, so that the determined fault plane is more accurate, and the workload is reduced.

Description

Three-dimensional fault modeling method and device

Technical Field

The disclosure relates to the technical field of three-dimensional geological modeling, in particular to a modeling method and device of a three-dimensional fault.

Background

The three-dimensional geological modeling plays an important role in geological exploration understanding, geological evolution process analysis and the like, and the establishment of the three-dimensional geological model is an important means in the field of geoscience and relates to evaluation of natural resources, disaster simulation prediction and the like. Because faults cut rock units and play an important role in natural phenomena in the subsurface, the construction of fault structures plays a crucial role in three-dimensional geological modeling.

Disclosure of Invention

The present disclosure is directed to solving, at least in part, one of the technical problems in the related art.

The embodiment of the first aspect of the disclosure provides a method for modeling a three-dimensional fault, which includes:

acquiring the type of a fault and the fault distance corresponding to each marking point on a trace of the fault;

determining an upper coiling line and a lower coiling line of the fault according to the trace and the fault distance corresponding to each marking point;

determining a fault boundary line according to the upper coil line, the lower coil line and a preset offset distance;

and determining a fault plane according to the type of the fault, the upper coil line, the lower coil line and the fault boundary line.

Optionally, the acquiring the type of the fault and the fault distance corresponding to each marking point on the trace of the fault includes:

determining a maximum fault-section marking point on the trace and a corresponding maximum fault-section;

acquiring the distance between each mark point and the maximum fault-section mark point;

and determining the fault distance corresponding to each marking point on the trace according to the maximum fault distance and the distance between each marking point and the maximum fault distance marking point based on a distance inverse ratio method.

Optionally, the determining an upper coil and a lower coil of the fault according to the trace and the fault distance corresponding to each mark point includes:

Determining a plurality of reference points corresponding to the upper coil wire and a plurality of reference points corresponding to the lower coil wire according to the trace and the fault distance corresponding to each marking point;

generating the upper winding according to the positions of a plurality of reference points corresponding to the upper winding;

and generating the lower winding line according to the positions of the reference points corresponding to the lower winding line.

Optionally, the determining a fault boundary line according to the upper and lower coil lines and a preset offset distance includes:

determining that the distance between the upper plate line and the upper plate line is the preset offset distance, and the boundary line far away from the trace line is the upper plate fault boundary line;

determining that the distance between the lower disc line and the lower disc line is the preset offset distance, and the boundary line far away from the trace line is the lower disc fault boundary line;

and determining the fault boundary line according to the upper disc fault boundary line and the lower disc fault boundary line.

Optionally, the determining a fault plane according to the type of the fault, the upper disc line, the lower disc line, and the fault boundary line includes:

in response to the fact that the fault type is a positive fault, determining a plane where the upper coil line and the upper coil fault boundary line are located as an upper disc surface, and determining a plane where the lower coil line and the lower coil fault boundary line are located as a lower disc surface;

And determining the fault plane according to the upper disc surface and the lower disc surface.

Optionally, the determining a fault plane according to the type of the fault, the upper disc line, the lower disc line and the fault boundary line includes:

in response to the fault type being a reverse fault, determining a plane where the upper coil line and the lower coil fault boundary line are located as an upper disc surface, and determining a plane where the lower coil line and the upper coil fault boundary line are located as a lower disc surface;

Optionally, the method further includes:

in response to the intersection of the traces of any two faults, determining the fault plane corresponding to the two faults according to the traces of any two faults.

Optionally, the method further includes:

respectively acquiring a first lower disc surface corresponding to a first fault and a second upper disc surface corresponding to a second fault in response to the fact that the traces of any two faults do not intersect and the distance between the two traces is smaller than a threshold value;

and determining a transition surface according to the first lower disk surface and the second upper disk surface.

Optionally, the method further includes:

acquiring a ground area boundary, the type of each fault contained in the ground area boundary, a trace of each fault, and an upper coil line and a lower coil line corresponding to each fault;

Determining a plurality of sub-regions contained in the ground region boundary according to the type of each fault, the trace of each fault, and the upper coil and the lower coil corresponding to each fault;

and performing Kriging interpolation on the plurality of sub-regions to determine a continuous region consisting of the plurality of sub-regions.

The embodiment of the second aspect of the present disclosure provides a three-dimensional fault modeling apparatus, including:

the first acquisition module is used for acquiring the type of a fault and the fault distance corresponding to each marking point on a trace of the fault;

the first determining module is used for determining an upper coil line and a lower coil line of the fault according to the trace and the fault distance corresponding to each marking point;

the second determining module is used for determining a fault boundary line according to the upper coil line, the lower coil line and a preset offset distance;

and the third determining module is used for determining the fault plane according to the type of the fault, the upper coil line, the lower coil line and the fault boundary line.

Optionally, the first obtaining module is specifically configured to:

Optionally, the first determining module is specifically configured to:

generating the upper wire winding according to the positions of a plurality of reference points corresponding to the upper wire winding;

and generating the lower winding line according to the positions of the plurality of reference points corresponding to the lower winding line.

Optionally, the second determining module is specifically configured to:

Optionally, the third determining module is specifically configured to:

Optionally, the system further includes a fourth determining module, specifically configured to:

in response to the intersection of the traces of any two faults, determining the fault plane corresponding to the two faults according to the traces of the any two faults.

Optionally, the system further includes a fifth determining module, specifically configured to:

Optionally, the method further includes:

the second acquisition module is used for acquiring a ground area boundary, the type of each fault contained in the ground area boundary, a trace of each fault, and an upper coil and a lower coil corresponding to each fault;

A sixth determining module, configured to determine, according to the type of each fault, the trace of each fault, and an upper tray line and a lower tray line corresponding to each fault, a plurality of sub-regions included in the ground region boundary;

and the seventh determining module is used for performing Kriging interpolation on the plurality of sub-areas so as to determine a continuous area formed by the plurality of sub-areas.

An embodiment of a third aspect of the present disclosure provides an electronic device, including: the three-dimensional fault modeling method comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor executes the program to realize the three-dimensional fault modeling method provided by the embodiment of the first aspect of the disclosure.

A fourth aspect of the present disclosure provides a computer-readable storage medium, which stores a computer program, and when the computer program is executed by a processor, the computer program implements the method for modeling a three-dimensional fault as set forth in the first aspect of the present disclosure.

A fifth aspect of the present disclosure provides a computer program product, wherein when the instructions of the computer program product are executed by a processor, the method for modeling a three-dimensional fault as set forth in the first aspect of the present disclosure is implemented.

The modeling method and the device for the three-dimensional fault have the following beneficial effects:

in the embodiment of the disclosure, the type of the fault and the fault distance corresponding to each marking point on the trace of the fault are obtained, then the upper coil and the lower coil of the fault are determined according to the trace and the fault distance corresponding to each marking point, the fault boundary line is determined according to the upper coil, the lower coil and the preset offset distance, and finally the fault plane is determined according to the type of the fault, the upper coil, the lower coil and the fault boundary line. Therefore, the upper coil line and the lower coil line of the fault are determined according to the fault distance corresponding to each marking point on the trace of the fault, and the fault plane is formed according to the upper coil line and the lower coil line, so that the determined fault plane is more accurate, and the workload is reduced.

Additional aspects and advantages of the disclosure will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the disclosure.

Drawings

The foregoing and/or additional aspects and advantages of the present disclosure will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings of which:

fig. 1 is a schematic flow chart of a method for modeling a three-dimensional fault according to an embodiment of the present disclosure;

FIG. 2 is a schematic view of a normal fault provided by an embodiment of the present disclosure;

FIG. 3 is a schematic view of a fault boundary line provided in an embodiment of the present disclosure;

FIG. 4 is a schematic flow chart diagram illustrating a method for modeling a three-dimensional fault according to another embodiment of the present disclosure;

FIG. 5 is a schematic view of another fault boundary line provided in an embodiment of the present disclosure;

FIG. 6 is a schematic view of a fault provided by an embodiment of the present disclosure;

FIG. 7 is a schematic diagram of a display interface for creating faults provided by embodiments of the present disclosure;

FIG. 8 is a schematic view of a two intersecting fault-forming fault boundary line provided by an embodiment of the present disclosure;

FIG. 9 is a schematic diagram of two non-intersecting positive fault forming transition surfaces provided by an embodiment of the present disclosure;

FIG. 10 is a schematic flow chart diagram illustrating a method for modeling a three-dimensional fault according to another embodiment of the present disclosure;

FIG. 11a is a schematic diagram of a plurality of faults included within a boundary of a ground area according to an embodiment of the present disclosure;

fig. 11b is a schematic diagram of a small region segmented within a boundary of a ground region according to an embodiment of the present disclosure;

FIG. 12 is a schematic structural diagram of a three-dimensional fault modeling apparatus provided in an embodiment of the present disclosure;

FIG. 13 illustrates a block diagram of an exemplary electronic device suitable for use in implementing embodiments of the present disclosure.

Detailed Description

Reference will now be made in detail to the embodiments of the present disclosure, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the drawings are exemplary and intended to be illustrative of the present disclosure, and should not be construed as limiting the present disclosure.

A three-dimensional fault modeling method and apparatus according to an embodiment of the present disclosure will be described below with reference to the drawings.

Fig. 1 is a schematic flow chart of a three-dimensional fault modeling method according to an embodiment of the present disclosure.

The embodiment of the disclosure exemplifies that the modeling method of the three-dimensional fault is configured in a modeling apparatus of the three-dimensional fault, which can be applied to any electronic device, so that the electronic device can perform a modeling function of the three-dimensional fault.

The electronic device may be a Personal Computer (PC), a cloud device, a mobile device, and the like, and the mobile device may be a hardware device having various operating systems, touch screens, and/or display screens, such as a mobile phone, a tablet Computer, a Personal digital assistant, a wearable device, and an in-vehicle device.

As shown in fig. 1, the method for modeling a three-dimensional fault may include the steps of:

step 101, acquiring the type of the fault and the fault distance corresponding to each marking point on the trace of the fault.

Fig. 2 is a schematic diagram of a normal fault according to an embodiment of the disclosure. As shown in FIG. 2, the fault distance is the horizontal distance between the upper tray and the lower tray, the inclination angle is the included angle between the fault plane and the horizontal plane, and the fall is the vertical distance between the upper tray and the lower tray.

The fault is a structure in which the earth crust is stressed to break and rock masses on two sides of the broken surface are subjected to obvious relative displacement.

Among them, the types of faults may include: normal faults, reverse faults, and the like. The present disclosure is not limited thereto.

Wherein the fault distance refers to the relative distance between the corresponding rock formations of the faulted rock formation on the two discs. The pitch values are different in the cross-sections of different orientations.

Additionally, the trace of the fault may be a line in the middle of the fault plane. The present disclosure is not limited thereto.

It should be noted that, since the rock near the fault line is easily weathered and eroded, which may cause that the mark points on the collected fault line are not easily measured, in the embodiment of the present disclosure, the upper winding line and the lower winding line may be determined by acquiring the fault distance corresponding to the mark point on the trace.

The mark points on the trace can be cosine marked points, and the distance between every two adjacent mark points can be the same or different. The present disclosure is not limited thereto.

Alternatively, the fault distance corresponding to each marking point on the trace of the fault may be measured in advance. And determining the corresponding fault distance of each marking point by using a distance inverse comparison method. The present disclosure is not limited thereto.

And 102, determining an upper coil line and a lower coil line of the fault according to the trace and the fault distance corresponding to each marking point.

It can be understood that, in the embodiment of the present disclosure, for a completely disconnected fault, the upper coil and the lower coil of the fault can be directly obtained by translating along the inclination direction. For faults that are not completely broken, the upper and lower coil lines may be determined from the fault distance corresponding to each marker point on the trace.

It should be noted that, since the trace is the line at the middle of the fault plane, the distance between the trace and the upper disc line is equal to the distance between the trace and the lower disc line. Therefore, after the break distance corresponding to each marking point on the trace is determined, a plurality of reference points corresponding to the upper wire rod and a plurality of reference points corresponding to the lower wire rod can be determined according to the trace and the break distance corresponding to each marking point, and then the upper wire rod is generated according to the positions of the reference points corresponding to the upper wire rod; and generating a lower coil according to the positions of the reference points corresponding to the lower coil.

And 103, determining a fault boundary line according to the upper coil line, the lower coil line and a preset offset distance.

For example, the preset offset distance may be 200m,300m, and so on. The present disclosure is not limited thereto.

The fault boundary line may be a range of faults, and may affect the shape of a fault plane.

It is understood that, since the fault plane is easily eroded by the weather, after the upper and lower routes are determined, the upper and lower routes may be shifted in a direction away from the trace line based on a preset shift distance to determine the fault boundary line. Wherein a schematic view of the fault boundary line may be as shown in fig. 3.

And step 104, determining a fault plane according to the type of the fault, the upper coil line, the lower coil line and the fault boundary line.

Optionally, the upper tray surface and the lower tray surface of the fault may be determined according to the type of the fault, the upper tray line, the lower tray line and the fault boundary line, and then the fault surface may be determined according to the upper tray surface and the lower tray surface.

It should be noted that, for the normal fault and the reverse fault, different regions are divided according to the fault boundary line and the characteristics (such as elevation, range, etc.) of the adjacent region to form a trend surface, that is, region analysis interpolation is performed on both sides of the fault.

In the embodiment of the disclosure, the type of the fault and the fault distance corresponding to each marking point on the trace of the fault are acquired, then the upper coil and the lower coil of the fault are determined according to the trace and the fault distance corresponding to each marking point, the fault boundary line is determined according to the upper coil, the lower coil and the preset offset distance, and finally the fault plane is determined according to the type of the fault, the upper coil, the lower coil and the fault boundary line. Therefore, the upper coil line and the lower coil line of the fault are determined according to the fault distance corresponding to each marking point on the trace of the fault, and the fault plane is formed according to the upper coil line and the lower coil line, so that the determined fault plane is more accurate, and the workload is reduced.

Fig. 4 is a schematic flowchart of a method for modeling a three-dimensional fault according to an embodiment of the present disclosure, and as shown in fig. 4, the method for modeling a three-dimensional fault may include the following steps:

step 401, determine the maximum distance marking point on the trace and the corresponding maximum distance.

And the maximum fault-section mark point is the mark point with the maximum fault-section on the trace.

Optionally, the fall corresponding to the maximum fault-distance marking point and the inclination corresponding to the fault can be measured first, and then the maximum fault-distance is determined according to the fall and the inclination.

Step 402, obtaining the distance between each mark point and the maximum fault-distance mark point.

The distance between each marking point and the maximum fault-section marking point can be a straight line distance.

And step 403, based on a distance inverse comparison method, determining the fault distance corresponding to each mark point on the trace according to the maximum fault distance and the distance between each mark point and the maximum fault distance mark point.

The distance inverse proportion method can be that the larger the distance between the mark point and the maximum fault-section mark point is, the smaller the fault-section corresponding to the mark point is; the smaller the distance between the mark point and the maximum fault-section mark point is, the larger the fault-section corresponding to the mark point is.

In the embodiment of the present disclosure, the determined distance between each marked point on the trace may be as shown in table 1 by using an inverse distance comparison method.

TABLE 1

Wherein, the distance from the mark point 2 to the maximum fault-section mark point is 0, which indicates that the mark point 2 is the maximum fault-section mark point; the distance between the mark point 0 and the mark point 4 is farthest from the mark point with the largest breaking distance, and the corresponding breaking distance is 0, namely the upper disc and the lower disc are not disconnected at the mark point 0 and the mark point 4; and determining the corresponding fault distances of the marking

points

2 and 4 by using a distance inverse comparison method.

It is understood that each element and each corresponding relationship in table 1 exist independently; these elements, correspondences, are exemplarily listed in the same table, but do not represent that all elements, correspondences, in the table must coexist according to the representation in table 1. Wherein the value of each element and each correspondence is independent of any other element value or correspondence in table 1. Therefore, as can be understood by those skilled in the art, the value of each element and each corresponding relationship in table 1 are independent embodiments.

In the embodiment of the disclosure, the fault distance corresponding to each marking point on the trace of the fault is determined based on a distance inverse comparison method, so that the problem that data such as fall or fault distance corresponding to each marking point need to be measured when a worker measures in the field is avoided, and the workload of the worker is reduced.

And step 404, determining an upper coil line and a lower coil line of the fault according to the trace and the fault distance corresponding to each marking point.

The specific implementation form of step 404 may refer to the detailed description in other embodiments of the present disclosure, and is not described in detail here.

Optionally, after the upper coil line and the lower coil line of the fault are determined, parameters such as the upper coil line, the lower coil line and the fault distance may be stored in the database, so that the data can be saved and called conveniently.

Step 405, determining a fault boundary line according to the upper coil line, the lower coil line and a preset offset distance.

Optionally, determining that the distance between the upper disc line and the upper disc line is a preset offset distance, and a boundary line far away from the trace line is an upper disc fault boundary line; determining that the distance between the lower disc line and the lower disc line is a preset offset distance, and the boundary line far away from the trace line is the lower disc fault boundary line; and determining a fault boundary line according to the upper disc fault boundary line and the lower disc fault boundary line.

Fig. 5 is a schematic view of another fault boundary line provided in an embodiment of the present disclosure. As shown in fig. 5, a line after shifting the upper disc line in a direction away from the trace by a preset shift distance is an upper disc fault boundary line; and a line obtained after the lower disc line is shifted by a preset offset distance in a direction away from the trace line is taken as a lower disc fault boundary line, and the upper disc fault boundary line and the lower disc fault boundary line are connected to obtain a fault boundary line.

And step 406, determining a fault plane according to the fault type, the upper coil line, the lower coil line and the fault boundary line.

Optionally, in response to that the fault type is a positive fault, determining a plane where the upper coil line and the upper coil fault boundary line are located as an upper disc surface, determining a plane where the lower coil line and the lower coil fault boundary line are located as a lower disc surface, and then determining a fault surface according to the upper disc surface and the lower disc surface.

Optionally, in response to that the fault type is a reverse fault, determining a plane where the upper coil line and the lower coil fault boundary line are located as an upper disc surface, determining a plane where the lower coil line and the upper coil fault boundary line are located as a lower disc surface, and then determining a fault surface according to the upper disc surface and the lower disc surface.

Fig. 6 is a schematic diagram of a fault according to an embodiment of the disclosure. As shown in fig. 6, wherein 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 are boundary points that can generate area trends. For the normal fault, the upper disc line forming area is an area enclosed by 1, 2, 3, 4, 5 and 6, and the surface of the area can be the upper disc surface; the area formed by the lower coiling line can be the area enclosed by 4, 7, 6, 8, 9 and 10, and the surface of the area can be the lower plate surface.

For the reverse fault, the area formed by the upper disc line is the area enclosed by 1, 2, 3, 4, 7 and 6, and the surface of the area can be the lower disc surface; the areas formed by the lower coiling lines are 6, 5, 4, 10, 9 and 8, and the surface of the area can be the upper coiling surface.

It can be understood that due to the particularity of the geological fault, after the fault distance corresponding to each marking point on the trace is determined by the over-distance inverse ratio method, the determined upper trace and the lower trace may have a difference from the actual situation according to the fault distance corresponding to each marking point. Therefore, the distance between the two lines can be determined by the inverse distance comparison method, and then the user can manually adjust the distance between the two lines, so as to adjust the upper line and the lower line.

Fig. 7 is a schematic diagram of a display interface for creating a fault according to an embodiment of the present disclosure. As shown in fig. 7, in creating a fault display interface, parameters such as a fault distance and a fall can be manually adjusted, so that an upper coil and a lower coil are adjusted.

It should be noted that, if two or more faults are located closer to each other, the fault boundary is very complicated. Fault boundaries formed in an offset manner interfere with each other and a conforming actual fault plane cannot be formed. Accordingly, embodiments of the present disclosure propose methods of forming fault planes from multiple faults together, such that complex fault boundaries may be formed.

In the embodiment of the present disclosure, two cases, that is, two traces intersect with each other and two traces do not intersect with each other, can be divided for two adjacent faults.

Optionally, in a case where the traces of any two faults intersect, the fault plane corresponding to the two faults may be determined according to the traces of any two faults.

Fig. 8 is a schematic diagram of a boundary line of two intersecting faults in accordance with an embodiment of the present disclosure. As shown in fig. 8, the trace 1 and the trace 2 are two intersecting traces, and the two faults may be multiplied together by a fault boundary line and then divided into different regions along the intersection to form a fault plane.

Optionally, under the condition that the traces of any two faults do not intersect and the distance between the two traces is smaller than a threshold, a first lower disk surface corresponding to the first fault and a second upper disk surface corresponding to the second fault are respectively obtained; and determining a transition surface according to the first lower disk surface and the second upper disk surface.

Fig. 9 is a schematic diagram of a transition surface formed by two non-intersecting positive faults according to an embodiment of the present disclosure. As shown in fig. 9, the first fault and the second fault are both positive faults, the upper tray surface of the first fault and the lower tray surface of the second fault are formed separately, the first lower tray surface of the first fault and the second upper tray surface of the second fault form transition surfaces together, and the region enclosed by 1, 2, 3, 4, 5, 6, 7 and 8 shown in fig. 9 is the transition surface.

In the embodiment of the disclosure, the maximum fault-section marking point and the corresponding maximum fault-section distance on the trace are determined, then the distance between each marking point and the maximum fault-section distance marking point is obtained, the fault-section distance corresponding to each marking point on the trace is determined according to the maximum fault-section distance and the distance between each marking point and the maximum fault-section distance marking point based on a distance inverse comparison method, then the upper coil and the lower coil of the fault are determined according to the trace and the fault-section distance corresponding to each marking point, the fault boundary line is determined according to the upper coil, the lower coil and the preset offset distance, and finally the fault plane is determined according to the type of the fault, the upper coil, the lower coil and the fault boundary line. Therefore, the fault distance corresponding to each marking point on the trace is determined according to the distance inverse proportion method, and the upper coil and the lower coil of the fault are determined according to the fault distance corresponding to each marking point, so that the determined fault plane can be more accurate, and the workload is further reduced.

Fig. 10 is a schematic flowchart of a method for modeling a three-dimensional fault according to another embodiment of the present disclosure, and as shown in fig. 10, the method for modeling a three-dimensional fault may include the following steps:

step 1001, a ground area boundary, a type of each fault included in the ground area boundary, a trace of each fault, and an upper coil and a lower coil corresponding to each fault are acquired.

Wherein, the ground area boundary can be a preset area boundary. The size of the ground area boundary is not limited in the embodiments of the present disclosure. For example, the land area boundaries may be 1000m × 1000m,1000m × 500m, and so on.

It should be noted that a plurality of faults may be included within the boundary of the ground area. Alternatively, the upper and lower coil lines corresponding to each fault may be determined according to the trace of each fault. Alternatively, if the upper and lower coils of any of the plurality of faults are stored in the database, the upper and lower coils corresponding to any fault may be directly acquired from the database.

FIG. 11a is a schematic diagram of a ground area including multiple faults within a boundary thereof according to an embodiment of the present disclosure. As shown in FIG. 11a, fault F1 is the reverse fault, and fault F2, fault F3, and fault F4 are the normal faults. Wherein fault F2 and fault F3 intersect.

Step 1002, determining a plurality of sub-regions included in the ground region boundary according to the type of each fault, the trace of each fault, and the upper coil and the lower coil corresponding to each fault.

Fig. 11b is a schematic diagram of a small region segmented within a boundary of a ground region according to an embodiment of the present disclosure. As shown in fig. 11b, scanning is performed from left to right in sequence from the boundary of the ground area, and F1 is a reverse fault, so that a first area is divided along the upper disc line of the fault F1 and the extension line of the trace, and the first area is an area surrounded by 1, 2, 3, 4, 5 and 6 in fig. 11 b; then, starting from F1, searching along the scanning direction, searching a fault F2 and a fault F3, because F2 and F3 are crossed, dividing regions jointly, namely extending the trace of the fault F2 to the lower disc line of the fault F1, extending the trace of the fault F3 to the boundary of a ground region, and dividing a region two and a region three, wherein the region two is the region surrounded by 6, 7, 9, 13, 12 and 11 in FIG. 11b, and the region three is the region surrounded by 4, 8, 9, 14, 15 and 10 in FIG. 11 b; continuing the scan direction search from slice F3, slice F4 is found, and according to the trace extension of slice F4, area four and area five are divided, area four being enclosed by 15, 16, 18, 25, 20, 22, 23, 24 in fig. 11b, and area five being enclosed by 12, 17, 18, 19, 20, 21 in fig. 11 b.

The first area, the second area, the third area, the fourth area and the fifth area are a plurality of sub-areas contained in the ground area boundary.

Step 1003, performing kriging interpolation on the multiple sub-areas to determine a continuous area composed of the multiple sub-areas.

It will be appreciated that after determining that a plurality of sub-regions are contained within the boundary of the ground area, kriging interpolation may be performed on the plurality of sub-regions to connect the sub-regions into a large area, so that the fault plane contained within the boundary of the ground area may be clearly determined.

In the embodiment of the disclosure, a ground area boundary, a type of each fault contained in the ground area boundary, a trace of each fault, and an upper coil and a lower coil corresponding to each fault are obtained, then a plurality of sub-areas contained in the ground area boundary are determined according to the type of each fault, the trace of each fault, the upper coil and the lower coil corresponding to each fault, and finally kriging interpolation is performed on the plurality of sub-areas to determine a continuous area formed by the plurality of sub-areas. Therefore, the ground area boundary can be divided into a plurality of small areas according to the type of each fault in the ground area boundary, the trace line, the upper coil line and the lower coil line of each fault, and then interpolation is carried out on each small area, so that the fault plane contained in the ground area boundary can be accurately determined.

In order to realize the embodiment, the disclosure further provides a three-dimensional fault modeling device.

Fig. 12 is a schematic structural diagram of a three-dimensional fault modeling apparatus according to an embodiment of the present disclosure.

As shown in fig. 12, the modeling apparatus 1200 for a three-dimensional fault may include: a first obtaining module 1210, a first determining module 1220, a second determining module 1230, and a third determining module 1240.

The first acquisition module 1210 is configured to acquire a type of a fault and a fault distance corresponding to each marking point on a trace of the fault;

the first determining module 1220 is configured to determine an upper coil line and a lower coil line of a fault according to the trace and the fault distance corresponding to each marking point;

the second determining module 1230 is configured to determine a fault boundary line according to the upper coil line, the lower coil line and a preset offset distance;

and a third determining module 1240 for determining the fault plane according to the fault type, the upper coil, the lower coil and the fault boundary line.

Optionally, the first obtaining module 1210 is specifically configured to:

and based on a distance inverse ratio method, determining the fault distance corresponding to each marking point on the trace according to the maximum fault distance and the distance between each marking point and the maximum fault distance marking point.

Optionally, the first determining module 1220 is specifically configured to:

determining a plurality of reference points corresponding to the upper wire rod and a plurality of reference points corresponding to the lower wire rod according to the trace and the fault distance corresponding to each marking point;

generating an upper coil according to the positions of a plurality of reference points corresponding to the upper coil;

and generating a lower coil according to the positions of the reference points corresponding to the lower coil.

Optionally, the second determining module 1230 is specifically configured to:

determining that the distance between the upper disc line and the upper disc line is a preset offset distance, and the boundary line far away from the trace line is the upper disc fault boundary line;

determining that the distance between the lower disc line and the lower disc line is a preset offset distance, and the boundary line far away from the trace line is the lower disc fault boundary line;

and determining a fault boundary line according to the upper disc fault boundary line and the lower disc fault boundary line.

Optionally, the third determining module 1240 is specifically configured to:

in response to the fact that the fault type is a positive fault, determining a plane where an upper coil line and an upper coil fault boundary line are located as an upper disc surface, and determining a plane where a lower coil line and a lower coil fault boundary line are located as a lower disc surface;

Optionally, the third determining module 1240 is specifically configured to:

and in response to the intersection of the traces of any two faults, determining a fault plane corresponding to the two faults according to the traces of any two faults.

Optionally, the method further includes:

the second acquisition module is used for acquiring the ground area boundary, the type of each fault contained in the ground area boundary, the trace of each fault, and an upper coil and a lower coil corresponding to each fault;

the sixth determining module is used for determining a plurality of sub-regions contained in the ground region boundary according to the type of each fault, the trace of each fault, and the upper coil and the lower coil corresponding to each fault;

and the seventh determining module is used for performing kriging interpolation on the plurality of sub-regions to determine a continuous region formed by the plurality of sub-regions.

The functions and specific implementation principles of the above modules in the embodiments of the present disclosure may refer to the above method embodiments, which are not described herein again.

The modeling device for the three-dimensional fault of the embodiment of the disclosure firstly acquires the type of the fault and the fault distance corresponding to each marking point on the trace of the fault, then determines the upper coil and the lower coil of the fault according to the trace and the fault distance corresponding to each marking point, determines the fault boundary line according to the upper coil, the lower coil and the preset offset distance, and finally determines the fault plane according to the type of the fault, the upper coil, the lower coil and the fault boundary line. Therefore, the upper coil line and the lower coil line of the fault are determined according to the fault distance corresponding to each marking point on the trace of the fault, and the fault plane is formed according to the upper coil line and the lower coil line, so that the determined fault plane is more accurate, and the workload is reduced.

In order to implement the above embodiment, the present disclosure further provides an electronic device, including: the three-dimensional fault modeling method comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein when the processor executes the program, the three-dimensional fault modeling method provided by the previous embodiment of the disclosure is realized.

In order to achieve the above embodiments, the present disclosure further provides a computer readable storage medium storing a computer program, which when executed by a processor, implements the modeling method of a three-dimensional fault as set forth in the foregoing embodiments of the present disclosure.

In order to implement the above embodiments, the present disclosure also proposes a computer program product, wherein when instructions in the computer program product are executed by a processor, the modeling method of the three-dimensional fault proposed by the foregoing embodiments of the present disclosure is implemented.

FIG. 13 illustrates a block diagram of an exemplary electronic device suitable for use in implementing embodiments of the present disclosure. The electronic device 12 shown in fig. 13 is only an example, and should not bring any limitation to the functions and the scope of use of the embodiments of the present disclosure.

As shown in fig. 13, electronic device 12 is embodied in the form of a general purpose computing device. The components of computer device 12 may include, but are not limited to: one or more processors or processing units 16, a system memory 28, and a bus 18 that couples various system components including the system memory 28 and the processing unit 16.

Bus 18 represents one or more of any of several types of bus structures, including a memory bus or memory controller, a peripheral bus, an accelerated graphics port, and a processor or local bus using any of a variety of bus architectures. These architectures include, but are not limited to, Industry Standard Architecture (ISA) bus, Micro Channel Architecture (MAC) bus, enhanced ISA bus, Video Electronics Standards Association (VESA) local bus, and Peripheral Component Interconnect (PCI) bus, to name a few.

Computer device 12 typically includes a variety of computer system readable media. Such media can be any available media that is accessible by computer device 12 and includes both volatile and nonvolatile media, removable and non-removable media.

Memory 28 may include computer system readable media in the form of volatile Memory, such as Random Access Memory (RAM) 30 and/or cache Memory 32. The computer device 12 may further include other removable/non-removable, volatile/nonvolatile computer system storage media. By way of example only, storage system 34 may be used to read from and write to non-removable, nonvolatile magnetic media (not shown in FIG. 13 and commonly referred to as a "hard drive"). Although not shown in FIG. 13, a disk drive for reading from and writing to a removable, nonvolatile magnetic disk (e.g., a "floppy disk") and an optical disk drive for reading from or writing to a removable, nonvolatile optical disk (e.g., a Compact disk Read Only Memory (CD-ROM), a Digital versatile disk Read Only Memory (DVD-ROM), or other optical media) may be provided. In these cases, each drive may be connected to bus 18 by one or more data media interfaces. Memory 28 may include at least one program product having a set (e.g., at least one) of program modules that are configured to carry out the functions of embodiments of the disclosure.

A program/utility 40 having a set (at least one) of program modules 42 may be stored, for example, in memory 28, such program modules 42 including, but not limited to, an operating system, one or more application programs, other program modules, and program data, each of which examples or some combination thereof may comprise an implementation of a network environment. Program modules 42 generally perform the functions and/or methodologies of the embodiments described in this disclosure.

Computer device 12 may also communicate with one or more external devices 14 (e.g., keyboard, pointing device, display 24, etc.), with one or more devices that enable a user to interact with computer device 12, and/or with any devices (e.g., network card, modem, etc.) that enable computer device 12 to communicate with one or more other computing devices. Such communication may be through an input/output (I/O) interface 22. Moreover, computer device 12 may also communicate with one or more networks (e.g., a Local Area Network (LAN), a Wide Area Network (WAN), and/or a public Network such as the Internet) via Network adapter 20. As shown, network adapter 20 communicates with the other modules of computer device 12 via bus 18. It should be understood that although not shown in the figures, other hardware and/or software modules may be used in conjunction with computer device 12, including but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, and data backup storage systems, among others.

The processing unit 16 executes various functional applications and data processing, for example, implementing the methods mentioned in the foregoing embodiments, by running a program stored in the system memory 28.

According to the technical scheme, the type of the fault and the fault distance corresponding to each marking point on the trace of the fault are obtained, then the upper coil and the lower coil of the fault are determined according to the trace and the fault distance corresponding to each marking point, the fault boundary line is determined according to the upper coil, the lower coil and the preset offset distance, and finally the fault surface is determined according to the type of the fault, the upper coil, the lower coil and the fault boundary line. Therefore, the upper coil line and the lower coil line of the fault are determined according to the fault distance corresponding to each marking point on the trace of the fault, and the fault plane is formed according to the upper coil line and the lower coil line, so that the determined fault plane is more accurate, and the workload is reduced.

In the description herein, references to the description of the term "one embodiment," "some embodiments," "an example," "a specific example," or "some examples," etc., mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present disclosure. In this specification, the schematic representations of the terms used above are not necessarily intended to refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Moreover, various embodiments or examples and features of various embodiments or examples described in this specification can be combined and combined by one skilled in the art without being mutually inconsistent.

Furthermore, the terms "first", "second" and "first" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one of the feature. In the description of the present disclosure, "a plurality" means at least two, e.g., two, three, etc., unless explicitly specifically limited otherwise.

Any process or method descriptions in flow charts or otherwise described herein may be understood as representing modules, segments, or portions of code which include one or more executable instructions for implementing steps of a custom logic function or process, and alternate implementations are included within the scope of the preferred embodiment of the present disclosure in which functions may be executed out of order from that shown or discussed, including substantially concurrently or in reverse order, depending on the functionality involved, as would be understood by those reasonably skilled in the art of the embodiments of the present disclosure.

The logic and/or steps represented in the flowcharts or otherwise described herein, e.g., an ordered listing of executable instructions that can be considered to implement logical functions, can be embodied in any computer-readable medium for use by or in connection with an instruction execution system, apparatus, or device, such as a computer-based system, processor-containing system, or other system that can fetch the instructions from the instruction execution system, apparatus, or device and execute the instructions. For the purposes of this description, a "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device. More specific examples (a non-exhaustive list) of the computer-readable medium would include the following: an electrical connection (electronic device) having one or more wires, a portable computer diskette (magnetic device), a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disc read-only memory (CDROM). Additionally, the computer-readable medium could even be paper or another suitable medium upon which the program is printed, as the program can be electronically captured, via for instance optical scanning of the paper or other medium, then compiled, interpreted or otherwise processed in a suitable manner if necessary, and then stored in a computer memory.

It should be understood that portions of the present disclosure may be implemented in hardware, software, firmware, or a combination thereof. In the above embodiments, the various steps or methods may be implemented in software or firmware stored in memory and executed by a suitable instruction execution system. If implemented in hardware, as in another embodiment, any one or combination of the following techniques, which are known in the art, may be used: a discrete logic circuit having a logic gate circuit for implementing a logic function on a data signal, an application specific integrated circuit having an appropriate combinational logic gate circuit, a Programmable Gate Array (PGA), a Field Programmable Gate Array (FPGA), or the like.

It will be understood by those skilled in the art that all or part of the steps carried by the method for implementing the above embodiments may be implemented by hardware related to instructions of a program, which may be stored in a computer readable storage medium, and when the program is executed, the program includes one or a combination of the steps of the method embodiments.

In addition, functional units in the embodiments of the present disclosure may be integrated into one processing module, or each unit may exist alone physically, or two or more units are integrated into one module. The integrated module can be realized in a hardware mode, and can also be realized in a software functional module mode. The integrated module, if implemented in the form of a software functional module and sold or used as a stand-alone product, may also be stored in a computer readable storage medium.

The storage medium mentioned above may be a read-only memory, a magnetic or optical disk, etc. While embodiments of the present disclosure have been shown and described above, it will be understood that the above embodiments are exemplary and not to be construed as limiting the present disclosure, and that changes, modifications, substitutions and alterations may be made to the above embodiments by those of ordinary skill in the art within the scope of the present disclosure.

Claims

1. A method of modeling a three-dimensional fault, comprising:

2. The method of claim 1, wherein the acquiring the type of fault and the fault distance corresponding to each marker point on the trace of the fault comprises:

3. The method of claim 1, wherein said determining an upper trajectory and a lower trajectory of said fault from said trace and a corresponding fault distance of said each marker point comprises:

4. The method of claim 1, wherein determining fault boundary lines based on the upper and lower trajectories and a predetermined offset distance comprises:

5. The method of claim 4, wherein determining a fault plane based on the type of fault, the upper disc line, the lower disc line, and the fault boundary line comprises:

6. The method of claim 4, wherein determining a fault plane based on the type of fault, the upper disc line, the lower disc line, and the fault boundary line comprises:

7. The method of any of claims 1-6, further comprising:

8. The method of any of claims 1-6, further comprising:

9. The method of any of claims 1-6, further comprising:

performing Krigin interpolation on the plurality of sub-regions to determine a continuous region composed of the plurality of sub-regions.

10. A modeling apparatus for a three-dimensional fault, comprising: