JP2006009262A - System for creating and displaying geologic structure model, stratum modeling method, and program for creating and displaying geologic structure model - Google Patents

Abstract

<P>PROBLEM TO BE SOLVED: To provide a system, a method and a program for creating and displaying a geologic structure model, which enable the removal of an unnecessary stratum section without designation and selection by a user. <P>SOLUTION: This system for creating and displaying the geologic structure model comprises: a stratum boundary surface generating and computing part which determines the depth distribution of a stratum boundary surface in the above area by prescribed computing from a plurality of coordinate values in an evaluative object area to be input, and geologic data for each predetermined depth distance at each of the coordinate values; an intersecting point computing part for determining a nodal line of the two stratum boundary surfaces, when the stratum boundary surfaces cross each other; an area computing and splitting part for splitting each of the crossing stratum boundary surfaces as different split boundary surfaces by the nodal line; and a correlation determining and computing part for extracting the split boundary surface incapable of existing, with respect to relevant information indicating a relationship between the upper and lower positions of the crossing strata. <P>COPYRIGHT: (C)2006,JPO&NCIPI

Description

  According to the present invention, a geological structure model is created in which a predetermined calculation is performed by estimating and extracting a formation boundary surface from geological data obtained by boring excavation, and this model is displayed three-dimensionally on a computer screen. The present invention relates to a display system, a geological modeling method, and a geological structure model creation / display program.

When building a building, knowing the stratum of the land to be built is a very important condition for increasing the strength of the building.
For this reason, geological structure modeling is conducted by boring excavation at a predetermined density in a target region, and geological information (related to the position of discrete specific stratum boundaries) at a certain depth unit for each excavation position. Information: coordinate values and depth direction geological data), and the shape of the formation is estimated and extracted.
The above-described estimation needs to reconstruct and display many measurement points and three-dimensional measurement data in the depth direction in the shape of the formation.
For this reason, a system that estimates the formation shape by calculation by inputting the data and displays it three-dimensionally (see, for example, Non-Patent Document 1).
3D underground information management system, vulcan (registered trademark), [online], [March 18, 2004 search], Internet, <URL: http://www.vulcan3d.com>

However, when the geological structure model creation / display system shown in Non-Patent Document 1 generates a geological boundary surface indicating a geological structure, when showing the interrelationship between the surface and fault inconsistency, which is a phenomenon peculiar to the geological field, The user needs to edit the image of the inconsistent portion.
Here, it is assumed that the right side of the formations T1, T2, and T3 has risen due to crustal movement, for example, as shown in FIG.
Then, after a certain period of time, as shown in FIG. 11 (c), a part of the formations T1 and T2 is removed by erosion due to rain or wind, and as a result, the ground surface B (the formation boundary surface B) is formed. Suppose.

Using the conventional geological structure model creation / display system (CAD system for estimating and displaying the stratum structure), when the stratum boundary surface formed by the various layers T1 to T3 described above is estimated and calculated, In the calculation process of the structure with C, it is not possible to determine that the formation boundary surface C is missing at the portion where the formation boundary surface B intersects and is only the formation boundary surface B. Therefore, as shown in FIG. A three-dimensional stratum in which a stratum boundary C that does not exist is formed on the top of the stratum boundary B, and the two stratum boundaries intersect, even though there cannot be a stratum boundary above the stratum boundary B The modeled image data will be displayed.
In other words, in the conventional formation CAD system, a portion of the formation boundary surface C that does not actually exist is formed below the formation boundary surface C.

For this reason, in the conventional geological CAD system, on the display screen, the user selects (i) two planes for editing the hierarchical relation of the geological boundary plane (ii) extraction of the intersection of the two planes (iii) two geological formations There is a disadvantage that an unnecessary part of the boundary surface is designated and selected (iv) The above unnecessary part removal operation must be performed.
In addition, as described above, the conventional geological CAD system selects and removes unnecessary portions by the user. However, since it is three-dimensional image data, it exists in the shadow of other plane images. However, it is difficult for the user to detect unnecessary portions, and there is a problem that a portion that cannot be removed occurs, causing inconvenience in later processing.

  The present invention has been made in view of such circumstances, and a geological structure model creation / display system, method and program capable of removing unnecessary portions without designation and selection by a user. The purpose is to provide.

  The geological structure model creating / displaying system according to the present invention is configured to input a depth distribution of a stratum boundary surface in a predetermined region from information (coordinate values) regarding the position of a specific discrete stratum boundary in an input evaluation target region. When the formation boundary surface generation calculation unit obtained by the calculation of the above and the formation boundary surfaces intersect each other, the intersection point / intersection line calculation for obtaining the intersection line of these two surfaces, and each of the intersecting formation boundary surfaces, An area calculation division unit that divides as different division boundary surfaces by intersection lines, and a correlation determination calculation unit that extracts the division boundary surfaces that cannot exist for related information indicating the vertical relationship of intersecting formations, It is characterized by having.

  In the geological structure model creation / display system according to the present invention, the related information is a relational expression “symmetric plane> constraint plane” indicating that a constraint plane that is a formation cannot exist above the symmetry plane, or a constraint plane Is represented by a relational expression “symmetric plane <constraint plane” indicating that the symbol cannot exist below the symmetry plane.

  The geological structure model creation / display system of the present invention is characterized in that it has an unnecessary part erasure calculation unit that deletes the division boundary surface and reconstructs the displayed image.

  The geological structure model creation / display system according to the present invention includes the formation boundary surface generation calculation unit, the calculated formation boundary surface for each region in a predetermined range, and the geological data for each predetermined depth distance at each coordinate value. From the above, the related information is calculated.

According to the formation modeling method of the present invention, the depth distribution of the formation boundary surface in the region is obtained by a predetermined calculation from the information (coordinate values) regarding the position of the discrete specific formation boundary in the evaluation target region. When the formation boundary surface generation calculation process and the formation boundary surfaces intersect with each other, the intersection point / intersection line calculation process for obtaining the intersection line of these two surfaces, and the intersecting formation boundary surface respectively by the intersection line A region calculation division process to divide as different division boundary surfaces;
Related information indicating the constraint relationship regarding the existence of intersecting strata (Relational expression “target plane> constraint plane” indicating that it cannot exist at the top, or that the constraint plane cannot exist below the target plane. A relational determination operation process for extracting the divided boundary surface that cannot exist, for a relational expression “target plane <constraint plane”).

  The geological structure model creating / displaying system according to the present invention determines the depth distribution of the formation boundary surface in the region from the information (coordinate values) regarding the position of the discrete specific formation boundary in the evaluation target region that is input. When the formation boundary surface generation calculation processing obtained by the calculation of the above and the formation boundary surfaces intersect each other, the intersection point / intersection line calculation processing for obtaining the intersection line of these two surfaces, and each of the intersecting formation boundary surfaces, Correlation determination that extracts the division boundary surface that cannot exist for the region calculation division processing that divides as a different division boundary surface by the intersection line and related information that indicates the constraint relationship regarding the existence of intersecting formations The program is executable by a computer having arithmetic processing.

  With the configuration described above, according to the geological structure model creation / display system of the present invention, the intersecting stratum boundary surfaces are each divided into divided boundary surfaces by intersecting lines, and the strata as constraint conditions set in advance are set. Based on the constraint relationship regarding the existence of the boundary surface, an unnecessary divided boundary surface is selected and this unnecessary portion is deleted. An effect is obtained that the divided boundary surface can be definitely and completely removed.

Hereinafter, a geological structure model creation / display system according to an embodiment of the present invention will be described with reference to the drawings. FIG. 1 is a block diagram showing a configuration example of a geological structure model creation / display system according to the embodiment.
In this figure, an input unit 1 stores geological data obtained from a geological survey by excavation of an evaluation target region by boring from an external device in a database 4 via a control unit 2. (Here, external devices are keyboards, external storage media, actual measuring devices, etc.)
Here, the geological data is used to create the formation boundary surface, and is information (coordinate values) on the position of the discrete specific formation boundary in the evaluation target area obtained from the geological survey such as boring excavation. It consists of and.
This boring excavation is performed at a predetermined density in the region to be evaluated, and the higher the density, the higher the accuracy of the obtained modeling result.

The output unit 3 inputs the formation structure modeled by the calculation of the calculation unit 5 via the control unit 2 and displays it on the display screen or the like.
The database 4 includes information (coordinate values) related to the position of the discrete specific formation boundary in the evaluation target area, the position information in the three-dimensional space of the formation boundary surface obtained by the calculation unit 5, and the formation The identifier and display color of the boundary surface, the related information indicating the vertical relationship between the boundary surfaces of each layer for each predetermined area, and information on the start point, end point, and angle of the tomographic surface are stored.

The calculation unit 5 includes a formation boundary surface generation calculation unit 6, an intersection / intersection line calculation unit 7, a region division calculation unit 8, an interrelation determination calculation unit, and an unnecessary part erasure calculation unit.
The formation boundary surface generation calculation unit 6 inputs a plurality of measurement coordinate values in the evaluation target area (or stored in the database 4), and geological data for each predetermined depth distance at each of these coordinate values. Thus, the depth distribution of the formation boundary surface in the region is obtained by a predetermined calculation.

Here, the formation boundary surface generation calculation unit 6 applies, for example, a region to be evaluated that is limited from information (coordinate values) regarding the position of a discrete specific formation boundary in the region to be evaluated obtained by a geological survey or the like. On the other hand, the depth distribution of the formation boundary surface is estimated by a general mathematical method, and the formation boundary surface is generated as an aggregate of triangular or quadrilateral polygons (mesh).
The above general mathematical enemy method uses the depth of the formation boundary point obtained from the geological information for each distance (depth) in the predetermined depth direction in each measurement coordinate value, and the measured formation boundary point The depth distribution of the formation boundary surface is estimated by complementing the depth of the formation boundary point of the coordinate value that has not been measured by the depth, and the curved surface or the like is expressed as a collection of polygons (mesh).

For example, an aggregate of triangular meshes is generated as information on the formation boundary surface from the inverse distance weighted method shown below.
^{Z = Σ (Zi / Ri n} ) / Σ (1 / Ri n)
Ri ⁿ = ((X−Xi) ² + (Y−Yi) ² ) ^0.5
Here, Xi, Yi, Zi are positional information (actually obtained information) about a specific stratum boundary obtained by geological surveys, etc., and X, Y, Z are coordinate values of the contacts constituting the desired mesh. Yes, n is a weighting factor based on the reciprocal of the distance. Here, the measurement coordinate location of latitude and longitude is indicated by X and Y, Z indicates the depth, and the (X, Y) depth Z of the coordinate value of each vertex of the triangular mesh is estimated (complemented).

When the intersection / intersection line calculation unit 7 detects that the formation boundary surface generated by the formation boundary surface generation calculation unit 6 intersects, the intersection / intersection line calculation unit 7 calculates the coordinates of the intersection line of the two formation boundary surfaces that intersect.
That is, when calculating the intersection line between two formation boundary surfaces, the intersection / intersection line calculation unit 7 uses the calculation formula shown in FIG. 2 for each triangular mesh unit constituting the formation boundary surface to calculate the coordinates of the intersection point. Perform the operation.
The intersection / intersection line calculation unit 7 calculates an intersection that intersects each triangular mesh (FIG. 3A) by using the respective arithmetic expressions shown in FIG. 2 (FIG. 3B).

The area division calculation unit 8 divides the triangular mesh in the vicinity of the intersection by connecting the intersection and the vertex of the nearest triangle mesh to form a sub-triangular mesh (FIG. 3C). The mesh structure representing the formation boundary is reconstructed.
In addition, the area division calculation unit 8 divides each intersecting sub-triangular mesh into an area 1 and an area 2 by an intersection line obtained by combining the intersection points, and each intersecting formation boundary surface is divided into 2 by an intersection line. Divide into two different boundary surfaces, sub-stratum boundary surfaces.

The correlation determination calculation unit 9 stores each sub-stratification by using the correlation information (related information) indicating constraints on the existence of the layer boundary surface for each predetermined region input from the input unit and stored in the internal storage unit. For each boundary surface, the vertical relationship is determined, and the sub-stratum boundary surface is extracted as a surface that cannot exist.
The unnecessary portion erasure calculation unit 10 deletes the sub-stratum boundary surface extracted as a surface that cannot exist, reconstructs a modeled image of the stratum from the remaining sub-stratum boundary surface, and displays it on the display screen. .

In addition, the interrelation information includes not only information on the constraint relationship based on the generation order between intersecting stratum boundary surfaces (information indicating whether the target surface should be above or below the constraint surface), Contains information about which one of the strata boundaries cuts the other.
Here, taking FIG. 11D as an example, the correlation information includes combination data indicating combinations of intersections of intersecting formations such as plane B and plane C, and a constraint relationship regarding the existence of these combinations. It is composed of a relational expression (target plane> constraint plane or target plane <constraint plane) indicating data (indicating a constraint relation regarding existence between two intersecting stratum boundary surfaces).

In the formula of B (target surface)> C (restricted surface), the formation boundary surface that is the target surface is shown on the left side, and the formation boundary surface that is the restriction surface is shown on the right side.
The target plane is a formation boundary surface that is a target whose presence or absence is detected, and the constraint surface is a formation boundary surface that serves as a reference for the symmetry plane.
Further, although the target surface is a target to be deleted when the existence is denied, the constraint surface is a reference for detecting the presence or absence of the target surface, and thus is not a target to be deleted by the unnecessary partial erasure calculation unit 10.

For example, C <B indicates that the formation boundary surface C that is the target surface cannot exist above the formation boundary surface B that is the restriction surface. Accordingly, it is detected that the formation boundary surface C in the dotted line portion above the formation boundary surface B is a surface that does not exist.
There are four types of combinations shown in FIG. 4. As shown in FIG. 4A, the region boundary calculation unit 8 divides the formation boundary surface B into sub-stratum boundary surfaces B1 and B2, and the formation boundary surface C Suppose that it is divided into strata boundary surfaces C1 and C2.

In the case of B> C in FIG. 4 (b), the correlation determination calculation unit 9 cannot make the formation boundary surface B, which is the target surface, below the formation boundary surface C of the restriction surface, as described above. , It is detected that the surface B does not exist below the surface C, and the unnecessary partial erasure calculation unit 10 deletes the surface B2 below the surface C2 as an unnecessary surface.
In the case of B <C in FIG. 4C, the interrelationship determining calculation unit 9 is a constraint relationship indicating that the formation boundary surface B that is the target surface cannot exist above the formation boundary surface C of the constraint surface. Based on the data, it is detected that the surface B does not exist above the surface C, and the unnecessary partial erasure calculation unit 10 deletes the surface B1 above the surface C1.

In the case of C> B in FIG. 4D, the interrelationship determining calculation unit 9 is a constraint relationship indicating that the formation boundary surface C that is the target surface cannot exist below the formation boundary surface B of the restriction surface. Based on the data, it is detected that the surface C does not exist below the surface B, and the unnecessary partial erasure calculation unit 10 deletes the surface C1 under the surface B1 as an unnecessary surface.
In the case of C <B in FIG. 4E, the interrelationship determining calculation unit 9 is a constraint relationship that indicates that the target boundary surface C cannot exist above the target layer boundary surface B. Based on the data, it is detected that the surface C does not exist above the surface B, and the unnecessary partial erasure calculation unit 10 deletes the surface C2 above the surface B2.

Next, the operation of the geological structure model creation / display system of this embodiment will be described with reference to FIGS. FIG. 5 is a flowchart showing an operation example of the geological structure model creation / display system according to the present embodiment.
In step S1, the control unit 2 determines the position coordinates (measurement coordinate values, for example, latitude and longitude information that are excavation positions by boring) of a specific (discrete) formation boundary obtained by a geological survey or the like; (X, Y) coordinates) and measurement data composed of geological information for each predetermined depth (distance in the depth direction of excavation by drilling; depth Z) at the position coordinates are input via the input unit 1, and each position coordinate Every time it is stored in the database 4.
In addition, the control unit 2 inputs the correlation information via the input unit 1 and stores it in the storage unit in the correlation determination calculation unit 9 for each position coordinate.

At this time, the control unit 2 detects the formation boundary surface at the boundary of each formation from the geological information for each depth Z of the input measurement data, and the various formations at this position coordinate by the detected formation boundary surface. The correlation information may be obtained and stored in the storage unit in the correlation determination calculation unit 9 for each position coordinate.
In the measurement data input process described above, the input unit 1 is a keyboard or an external storage medium, and the user inputs the measurement data and the correlation information by performing key input or reading from the external storage medium.

Next, in step S2, when the user performs modeling of a predetermined range of geological structure on the display screen of the geological structure model creation / display system, the control unit 2 selects the predetermined range so that the control unit 2 falls within this range. The measurement data corresponding to the measurement coordinate values existing in is read from the database 4.
In step S3, the formation boundary surface generation calculation unit 6 estimates the boundary depth distribution of each layer from the read data, and generates a formation boundary surface.
Further, the output unit 3 displays the obtained image of the formation structure including the formation boundary surfaces A, B, C, and D on the display screen as shown in FIG.

Next, in step S4, the intersection / intersection line calculation unit 7 intersects the stratum boundary surfaces intersecting each other, that is, the intersection point of the stratum boundary surfaces A and D, the intersection point of the stratum boundary surfaces B and D, and the stratum boundary The intersection of surfaces C and D is obtained, and the intersection line between each surface is obtained from this intersection.
Then, in step S5, the region division calculation unit 8 uses the intersecting lines to change the formation boundary surfaces A, B, C, and D into the sub-stratification boundary surfaces A1 and A2 and the sub-stratification boundary surface as shown in FIG. Divide into B1 and B2, sub-stratum boundary surfaces C1 and C2, and sub-stratum boundary surfaces D1 and D2.

Next, in step S6, the correlation determination calculation unit 9 determines that the generated geological structure image is based on the correlation information for each measurement coordinate value set in the storage unit in advance with the actual formation boundary surface. It is detected whether or not there is consistency.
For example, it is assumed that the correlation information stored in the storage unit is a, b, and c below.
a, C−D C <D The formation boundary surface C does not exist above the formation boundary surface D.
b, B−D B <D The formation boundary surface B does not exist above the formation boundary surface D.
c, A−D A <D The stratum boundary A does not exist above the stratum boundary D.
d, A−D D <A The formation boundary surface D does not exist above the formation boundary surface A.

First, the correlation determination calculation unit 9 uses the sub-stratum based on the constraint relationship data regarding the presence of the correlation information a (that is, the constraint condition that “the stratum boundary surface C does not exist above the stratum boundary surface D”). It is detected that the sub-stratum boundary surface C1 existing above the boundary surface D1 is unnecessary.
In step S7, the unnecessary portion erasure calculation unit 10 erases the image data of the sub-stratum boundary surface C1 based on the detection result, and removes the sub-stratum boundary surface C1 as shown in FIG. Reconfigure.

Next, in step S8, the correlation determination calculation unit 9 detects whether or not the correlation of all the correlation information corresponding to the measurement data in the selected range has been determined, and the correlation information b, c is still detected. , D is not completed, and the process returns to step S6.
In step S6, the correlation determination calculation unit 9 uses the constraint relationship data regarding the presence of the correlation information b (that is, the constraint condition that “the formation boundary surface B does not exist above the formation boundary surface D”). The sub-stratum boundary surface B1 existing above the sub-stratum boundary surface D1 is detected to be unnecessary because it is not consistent with the constraint conditions.

Next, in step S7, the unnecessary portion erasure calculation unit 10 erases the image data of the sub-stratum boundary surface B1 based on the detection result, and removes the sub-stratum boundary surface B1 as shown in FIG. Reconfigure.
Then, in step S8, the correlation determination calculation unit 9 detects whether or not the correlation of all the correlation information corresponding to the selected range of measurement data has been determined. It is detected that the process has not ended, and the process returns to step S6.

Next, in step S 6, the correlation determination calculation unit 9 performs constraint relationship data regarding the presence of the correlation information c (that is, a constraint condition that “the formation boundary surface A does not exist above the formation boundary surface D”). Thus, it is detected that the sub-stratum boundary surface A1 existing above the sub-stratum boundary surface D1 is unnecessary.
Then, in step S7, the unnecessary part erasure calculation unit 10 erases the image data of the sub-stratum boundary surface A1 based on the detection result, and removes the sub-stratum boundary surface A1 as shown in FIG. Reconfigure.

Next, in step S8, the correlation determination calculation unit 9 detects whether or not the correlation of all the correlation information corresponding to the measurement data in the selected range has been determined, and the correlation information d is still finished. If it is not detected, the process returns to step S6.
In step S6, the correlation determination calculation unit 9 uses the constraint relationship data related to the presence of the correlation information d (that is, the constraint condition that “the formation boundary surface D does not exist above the formation boundary surface A”). The sub-stratum boundary surface D2 existing above the sub-stratum boundary surface A2 is detected to be unnecessary.

Next, in step S7, the unnecessary part erasure calculation unit 10 erases the image data of the sub-stratum boundary surface D2 based on the detection result, and removes the sub-stratum boundary surface D2 as shown in FIG. Reconfigure.
In step S8, the correlation determination calculation unit 9 detects whether or not the correlation of all the correlation information corresponding to the measurement data in the selected range has been determined, and the correlation determination for all the correlation information is performed. If it is detected that the process has been completed, the process proceeds to step S9.

Next, in step S9, as shown in FIG. 10, the output unit 3 outputs image data of the finally modeled stratum structure to the display screen.
In step S <b> 10, the control unit 2 attaches a modeling identification number to the final modeled image data shown in FIG. 10 and stores it in the database 4.

In addition, the above-mentioned stratum boundary surface can include a ground surface or a fault surface, and the above-described correlation information is also defined for the ground surface or the fault surface, and is used as a target of the processing for determining the correlation.
For example, if the fault plane does not reach the ground due to inconsistency, the fault plane above the mismatch plane may be deleted in the modeled structure by defining the mismatch plane as a constraint. it can.

  Note that the program for realizing the functions of the formation CAD system in FIG. 1 is recorded on a computer-readable recording medium, and the program recorded on the recording medium is read into the computer system and executed, thereby executing the formation 3 of the formation. Dimensional image generation and display processing may be performed. Here, the “computer system” includes an OS and hardware such as peripheral devices. The “computer system” includes a WWW system provided with a homepage providing environment (or display environment). The “computer-readable recording medium” refers to a portable medium such as a flexible disk, a magneto-optical disk, a ROM, and a CD-ROM, and a storage device such as a hard disk built in the computer system. Further, the “computer-readable recording medium” refers to a volatile memory (RAM) in a computer system that becomes a server or a client when a program is transmitted via a network such as the Internet or a communication line such as a telephone line. In addition, those holding programs for a certain period of time are also included.

  The program may be transmitted from a computer system storing the program in a storage device or the like to another computer system via a transmission medium or by a transmission wave in the transmission medium. Here, the “transmission medium” for transmitting the program refers to a medium having a function of transmitting information, such as a network (communication network) such as the Internet or a communication line (communication line) such as a telephone line. The program may be for realizing a part of the functions described above. Furthermore, what can implement | achieve the function mentioned above in combination with the program already recorded on the computer system, and what is called a difference file (difference program) may be sufficient.

It is a block diagram which shows the structural example of the geological structure model creation and display system by one Embodiment of this invention. It is a conceptual diagram explaining the process which calculates | requires the intersection which two stratum boundary surfaces cross | intersect. It is a conceptual diagram explaining the process which calculates | requires the intersection line which two stratum boundary surfaces cross | intersect from the intersection calculated | required in description of FIG. It is a conceptual diagram explaining the structure and kind of correlation information of two stratum boundary surfaces. It is a flowchart which shows the operation example of the geological structure model creation and display system by one Embodiment of this invention. It is a conceptual diagram explaining the determination operation | movement of the mutual relationship (mainly restrictions relation regarding existence) by the adjacent formation boundary surface by the correlation determination calculating part 9 of FIG. It is a conceptual diagram explaining the determination operation | movement of the mutual relationship between the adjacent geological boundary surfaces by the correlation determination calculating part 9 of FIG. It is a conceptual diagram explaining the determination operation | movement of the mutual relationship between the adjacent geological boundary surfaces by the correlation determination calculating part 9 of FIG. It is a conceptual diagram explaining the determination operation | movement of the mutual relationship between the adjacent geological boundary surfaces by the correlation determination calculating part 9 of FIG. It is a conceptual diagram explaining operation | movement of a geological structure model creation and display system. It is a conceptual diagram explaining operation | movement of a geological structure model creation and display system.

Explanation of symbols

DESCRIPTION OF SYMBOLS 1 ... Input part 2 ... Control part 3 ... Output part 4 ... Database 5 ... Calculation part 6 ... Formation boundary surface generation calculation part 7 ... Intersection / intersection line calculation part 8 ... Area division | segmentation calculation part 9 ... Correlation determination calculation part 10 ... Unnecessary partial erasure calculation section

Claims (6)

From the information regarding the position of the discrete specific stratum boundary in the evaluation target region that is input, the stratum boundary surface generation calculation unit that obtains the depth distribution of the stratum boundary surface in the region by a predetermined calculation;
When the stratum boundary surfaces intersect each other, an intersection / intersection line calculation unit for obtaining an intersection line of these two surfaces;
An area division calculation unit that divides each intersecting formation boundary surface as a different division boundary surface by the intersection line;
A geological structure model creation / display system, comprising: an interrelation determination calculation unit that extracts the division boundary surface that cannot exist with respect to related information indicating a vertical relationship between intersecting strata.   The related information is a relational expression “target plane> constraint plane” indicating that a constraint plane that is a division boundary plane cannot exist above the target plane, or a constraint plane cannot exist below the target plane. The geological structure model creation / display system according to claim 1, wherein the system is represented by a relational expression “target plane <constraint plane”.   A geological structure model creating / displaying system comprising an unnecessary part erasing operation unit that deletes the divided boundary surface that cannot exist and reconstructs a displayed image.   The said formation boundary surface production | generation calculation part calculates the said relevant information from the information regarding the position of a formation boundary surface and a discrete specific formation boundary for every area | region of the predetermined range. The geological structure model creation / display system according to claim 3. Formation boundary surface generation calculation for obtaining a depth distribution of the formation boundary surface in the region by a predetermined operation from a plurality of coordinate values in the input evaluation target region and information on the position of the discrete specific formation boundary. Process,
When the stratum boundary surfaces intersect each other, an intersection / intersection calculation process for obtaining an intersection line between these two surfaces,
An area calculation division process of dividing each intersecting formation boundary surface as a different division boundary surface by the intersection line,
A formation modeling method, comprising: a correlation determination calculation process for extracting the division boundary surface that cannot exist with respect to related information indicating a constraint relationship regarding the existence of intersecting formations. The formation boundary that obtains the depth distribution of the formation boundary surface in the region by a predetermined calculation from a plurality of coordinate values in the input evaluation target region and the geological data for each predetermined depth distance at each coordinate value Surface generation calculation processing,
When the stratum boundary surfaces intersect each other, an intersection / intersection calculation process for obtaining an intersection line of these two surfaces;
Region calculation division processing for dividing each intersecting formation boundary surface as a different division boundary surface by the intersection line;
A computer-executable geological structure model creation / display program comprising: correlation information calculation processing for extracting the division boundary surface that cannot exist with respect to related information indicating a constraint relationship regarding the existence of intersecting strata.

Images