CN101038677B - Interface data conversion method and device - Google Patents
Interface data conversion method and device Download PDFInfo
- Publication number
- CN101038677B CN101038677B CN2006100114915A CN200610011491A CN101038677B CN 101038677 B CN101038677 B CN 101038677B CN 2006100114915 A CN2006100114915 A CN 2006100114915A CN 200610011491 A CN200610011491 A CN 200610011491A CN 101038677 B CN101038677 B CN 101038677B
- Authority
- CN
- China
- Prior art keywords
- point
- processing module
- data
- data file
- grid format
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Abstract
The invention discloses an interface data conversion method and an apparatus thereof, the method is suitable for conversion between data of a three-dimensional modeling processing module and an earthquake inversion processing module, the method comprises: 1, generating and outputting an enhanced raster-format data file of a three-dimensional geological model by the three-dimensional modeling processing module; 2, reading the enhanced raster-format data file by the earthquake inversion processing module, and performing an earthquake inversion operation based on the enhanced raster-format data file; wherein the enhanced raster-format data file has a discrete data structure with a data in X and Y directions coincident to earthquake data and a same number of rasters with the earthquake data, and contains structure information in which key points and point pairs of the three-dimensional geological model are continuously recorded in a Z direction. The enhanced raster-format data file provided in the invention realizes the conversion between interface data of the three-dimensional modeling processing module and the earthquake inversion processing module; describes the complex three-dimensional geological model by a simple structure, and is convenient to be used by the earthquake inversion processing module.
Description
Technical field
The present invention relates to a kind of interface data conversion method and device thereof, particularly relate to a kind of method and conversion equipment that interface data is changed between three-dimensional modeling processing module/software and the seismic inversion processing module/software that be used for.
Background technology
Three-dimensional modeling processing module/software adopts the geometric model of triangle gridding to describe complex geological tectonic.Because architectonic complicacy has determined geometric model to adopt complex topology structure record stratum shape and tomography shape.Simultaneously, geometric model also will write down the relation between stratum and stratum, stratum and tomography and tomography and the tomography.
General earthquake inversion processing module/software adopts discrete aspect Raster Data Model record tectonic structure information, and as the model of P701 form, this model all disperses on X, Y, three directions of Z.The aspect Raster Data Model is fairly simple, but can't write down complex geological tectonic information comprehensively.In the geological structure model that tomography exists, can't in the stratum, add a new layer according to the information among the P701, so the aspect Raster Data Model can not be described near the geologic structure of tomography exactly according to certain ratio.
If geometric model is converted to existing aspect Raster Data Model, as the model of P701 form, so a lot of important geology tectonic informations can be lost; If allow the seismic inversion processing module directly geometric model as input, need so the earthquake inversion processing module is done bigger modification.
Summary of the invention
Technical matters to be solved by this invention is to provide a kind of interface data conversion method and device thereof, is used for simply, conveniently carries out interface data conversion between three-dimensional modeling processing module/software and the seismic inversion processing module/software.
To achieve these goals, the invention provides a kind of interface data conversion method, be applicable to the conversion of data between three-dimensional modeling processing module and the seismic inversion processing module, it is characterized in that, comprise the steps:
Described enhancing grid format data file has on X, Y direction and geological data unanimity, discrete data structure that Gate lattice number is identical with geological data, and be included in key point and the right structural information of point that continuous recording on the Z direction has three-dimensional geological model, this key point is the incident line and the stratal surface of three-dimensional geological model and the intersection point of fault surface on the imaginary Z direction in grid node place, and it is right that per two the neighbouring key points in same grid node place are combined into a point; Described step 2 further comprises name, key point that type is identical is summed up as same ground prime element, and the identical have a few in affiliated layer position is to being summed up as same ground prime element.
Described interface data conversion method, wherein, described step 1 specifically comprises again: described three-dimensional modeling processing module is set up a three-dimensional geological model, mark off the grid consistent on the XY face in this geologic model, penetrate all stratum and tomography to the incident line at each imagination Z of grid node place with geological data; The intersection point of described every incident line and stratal surface and fault surface is a key point; And with in this three-dimensional geological model at same grid node place per two neighbouring key points be combined into a little right, note all key points, point to and be stored in the file that strengthens grid format, generate and strengthen the grid format data file.
Described interface data conversion method, wherein, described step 2 is specially: described seismic inversion processing module is read in the data and the right data of point of the key point of three-dimensional geological model from described enhancing grid format data file, and the data of described key point are carried out the seismic inversion computing with the right data of point as input.
Described interface data conversion method, wherein, described step 2 further comprises: the absolute value that the degree of depth of described some centering end point deducts the degree of depth of starting point is the thickness on each stratum, grid node place.
Described interface data conversion method, wherein, described step 2 further comprises: the step of asking for the property value of arbitrfary point in the three-dimensional geological model by described enhancing grid format data file in conjunction with the logging trace data.
Described interface data conversion method, wherein, the described step of asking for the property value of arbitrfary point in the three-dimensional geological model specifically comprises again:
Step 701 determines that according to the XYZ coordinate of this arbitrfary point its grid node of living in, point of living in are right, and further determines stratum of living in, this arbitrfary point and the ratio in stratum of living in;
Step 702 calculates a new aspect according to described ratio and stratum of living in; And
Step 703 is asked for the intersection point of logging trace track and this new aspect, obtains the property value of this arbitrfary point by interpolation according to the property value of described intersection point.
Described interface data conversion method wherein, in the described step 2, also comprises by described seismic inversion processing module and asking in the stratum and end face and the step of bottom surface apart from the calculating aspect that becomes fixed proportion according to described enhancing grid format data file.
Described interface data conversion method wherein, in the described step 2, also comprises by described seismic inversion processing module and asks in the stratum step with the equidistant calculating aspect of end face according to described enhancing grid format data file.
Described interface data conversion method wherein, in the described step 2, also comprises by described seismic inversion processing module and asks in the stratum step with the equidistant calculating aspect in bottom surface according to described enhancing grid format data file.
To achieve these goals, the present invention also provides a kind of interface data conversion equipment, it is characterized in that, comprise a three-dimensional modeling processing module, a seismic inversion processing module that is connected with described three-dimensional modeling processing module, described three-dimensional modeling processing module comprises an interface data modular converter, an interface data output module again;
Described interface data modular converter is used to generate the enhancing grid format data file of a three-dimensional geological model; Described interface data output module is used to export described enhancing grid format data file; Described seismic inversion processing module is used to read in described enhancing grid format data file, and carries out the seismic inversion computing according to described enhancing grid format data file;
Described enhancing grid format data file has on X, Y direction and geological data unanimity, discrete data structure that Gate lattice number is identical with geological data, and be included in key point and the right structural information of point that continuous recording on the Z direction has three-dimensional geological model, this key point is the incident line and the stratal surface of three-dimensional geological model and the intersection point of fault surface on the imaginary Z direction in grid node place, and it is right that per two the neighbouring key points in same grid node place are combined into a point; Described seismic inversion processing module is summed up as same ground prime element with name, key point that type is identical, and the identical have a few in affiliated layer position is to being summed up as same ground prime element.
Described interface data conversion equipment, wherein, described seismic inversion processing module is used to read the interface that strengthens Gate sound of laughing formula data file by one and reads in described enhancing grid format file.
Have the structural information that the data file that strengthens grid format has comprised necessity of three-dimensional geological model among the present invention, this enhancing grid format is a kind of simple, the discrete data layout that can express the three-dimensional geological model topological structure again exactly; This strengthens grid format can greatly make things convenient for the use of seismic inversion processing module with the three-dimensional geological model of simple structrual description complexity.
Describe the present invention below in conjunction with the drawings and specific embodiments, but not as a limitation of the invention.
Description of drawings
Fig. 1 strengthens grid format X, Y direction grid synoptic diagram for the present invention;
Fig. 2 strengthens the diagrammatic cross-section of the vertical Z direction of grid format for the present invention;
Fig. 3 strengthens the application synoptic diagram of grid format for the present invention;
Fig. 4 is the method flow diagram of interface data conversion between three-dimensional modeling processing module of the present invention and the seismic inversion processing module;
Fig. 5 is a specific embodiment of the present invention;
Fig. 6 asks for the process flow diagram of any some property values in the three-dimensional geological model for the present invention;
Fig. 7 is an interface data conversion equipment synoptic diagram of the present invention.
Embodiment
See also shown in Figure 1, for the present invention strengthens grid format X, Y direction grid synoptic diagram.Strengthening grid format is the architectonic data layout of a kind of description.As shown in Figure 1, on the XY direction, employing of the present invention and the consistent grid format of geological data (as SEGY and P701 form), the number of grid is identical with geological data, at the node place of each grid, suppose that Z is to there being a vertical incident line to penetrate all stratum and tomography.
See also shown in Figure 2ly, strengthen the diagrammatic cross-section of the vertical Z direction of grid format for the present invention.And in conjunction with Fig. 1, as shown in Figure 2, the intersection point of every incident line and stratal surface and fault surface is called key point P (Point), as P1, P2, P3, P4, P5, P6, P7, P8, P9, P10.The information of key point P comprises: coordinate X, Y, Z, intersection point type Type, the name Name of the aspect of handing over.Intersection point type Type is used to write down the intersection point of this intersection point and stratal surface or the intersection point of this intersection point and fault surface; The name Name of institute's friendship aspect is used to write down the name of stratum or tomography.Intersect as P1 and stratum L1, the name Name of P1 is L1.
Per two neighbouring key points at same grid node place are formed a point to PC (PointCouple), and point comprises the information of PC record: starting point Pup, end point Pdn, starting point ratio Rup, end point ratio Rdn, affiliated layer position Lname.Layer position Lname represents this aspect with the name of the end face of a layer position usually wherein.Key point, the right file layout of point are as follows:
Points
{
//No. x y z type Name
0， 12530，?3580，?-900.1，?L， layer1；
1， 12530， 3580， -990.5， F，?fault1；
}
Point?Couples
{
//Pup?/Pdn?/Rup /Rdn/LayerName
0， 1， 0， 1， layer1；
1， 2， 0， 1， layer2；
20030，20031，0，0.73，layer5；
30016，30017，0.335，1，layer6；
}
Strengthening grid format is a kind of discrete data structure that combines discrete and continuous characteristics: in the XY direction, data structure disperses, and is the grid format consistent with geological data, and the number of grid is identical with geological data; Make progress at Z, data structure be can be on the Z direction structural information of continuous recording three-dimensional geological model.The data that strengthen grid format comprise two kinds of main information: it is right that the key point on stratal surface and the fault surface reaches the point of being made up of these key points.These key points and the point to having comprised three-dimensional geologic structure structure of models information.Can combine closely three-dimensional modeling processing module/software and seismic inversion processing module/software (as Jason software) by the data that strengthen grid format.
See also shown in 3, strengthen the application synoptic diagram of grid format for the present invention.Shown in Fig. 1,2, this is used by using and strengthens the raster data form, can make things convenient for generation and end face and bottom surface apart from the calculating aspect that becomes fixed proportion, and equidistant calculating aspect of end face and the calculating aspect equidistant with the bottom surface.This application shown in Figure 3 30 comprises following several aspect:
31, ask for the calculating aspect that becomes fixed proportion with the distance of end face and bottom surface; Obtain the stratum title at the aspect place that will calculate and be in end face and the bottom surface between ratio, just can obtain the depth value of calculating aspect according to above-mentioned enhancing grid format data structure at each node place.
For example, ask among Fig. 2 in the L1 of stratum distance with end face and bottom surface than the calculating aspect that is 0.3.With some pt1 on the calculating aspect and pt3 is example, and the calculating aspect is arranged in stratum L1, and the ratio at top, the end of being in is 0.3 position:
Point pt1 on incident line S1, and be on the incident line S1 among the L1 of stratum point to for (P1, P2), the ratio 0.3 according to the degree of depth, ratio and the calculating aspect of P1 and P2 can obtain the depth value of calculating aspect at this point place:
Zpt1＝ZP1+(ZP2-ZP1)(Rpt1-RP1)/(RP2-RP1) (1)
If the value of Rpt1 is less than RP1, perhaps the value of Rpt1 illustrates the grid node place at this incident line S1 place greater than RP2, does not meet the point of ratio among the L1 of stratum.
Point pt3 on incident line S3, and be on the incident line S3 in the L1 stratum point to for (P6, P7), the ratio 0.3 according to the degree of depth, ratio and the calculating aspect of P6 and P7 can obtain the depth value of calculating aspect at this point place:
Zpt3＝ZP6+(ZP7-ZP6)(Rpt3-RP6)/(RP7-RP6) (2)
Calculate all grid node places successively with top end distance than the degree of depth that is 0.3 point, just obtained the aspect of asking.
32, ask for a calculating aspect equidistant with end face; Obtain the stratum title at the aspect place that will calculate and, just can obtain the depth value of calculating aspect according to above-mentioned enhancing grid format data structure at each grid node place from the distance of end face.
For example, the distance from end face is the calculating plane of d rice among the calculating stratum Layer1.
At some grid nodes place, the point that intersects of incident line S and stratum Layer1 to be (Pup, Pdn), calculate the depth value computing method of the some p on plane and be at this node place so:
|Zp|＝|Zup|+d (3)
If | Zp| shows at this grid node place greater than | Zdn|, among the Layer1 of stratum not the distance from end face be the point of d.
Calculate the depth value on the calculating plane at all node places successively.
33, ask for a calculating aspect equidistant with the bottom surface.
Similar with the method for asking for the equidistant calculating plane of end face, the formula that calculates with the equidistant aspect in bottom surface is:
|Zp|＝|Zdn|-d (4)
If | Zp| shows at this grid node place less than | Zup|, is not the point of d from the bottom surface distance among the Layer1 of stratum.
Use enhancing grid format that the present invention proposes and logging trace data can the interpolation three-dimensional geological model in the property value of any point p.Judge at first which bar incident line the p point is positioned on, then according to the depth value of p obtain p which stratum, which point between and ratio in the stratum.According to ratio and layer of living in position, just can calculate a new aspect.Try to achieve the intersection point of well track and this new aspect,, can interpolation obtain the property value that p is ordered according to the property value at these intersection point places.If p or well point not just in time at the grid node place, can obtain the property value of these points by the interpolation to grid node.
See also shown in Figure 4ly, be the method flow diagram of interface data conversion between three-dimensional modeling processing module of the present invention and the seismic inversion processing module.And in conjunction with shown in Figure 2, this flow process comprises the steps:
The step that generates the enhancing grid format data file of three-dimensional geological model specifically comprises:
(11), set up three-dimensional geological model (geometric model), at surface level XY upwards, the grid that employing of the present invention is consistent with geological data at the node place of each grid, supposes that Z is to there being a vertical incident line to penetrate all stratum and tomography.The intersection point of every incident line and stratal surface and fault surface goes on record, and these points are called key point P (Point), computation grid node promptly of the present invention on stratal surface and fault surface vertical projection key point P and write down the information of these key points;
(12), will be combined into a point at per two neighbouring key point P at same grid node place to PC (Point Couple), and write down this point right;
(13), write down all key points, the institute that writes down all grid node places is somewhat right, and is stored in the file or data structure that strengthens grid format, has so just generated the enhancing raster data from three-dimensional model.These data can be saved to strengthening the data file of grid format.
The point of incident line S1 to be (P1, P2), starting point Pup=P1, end point Pdn=P2, starting point ratio Rup=0, end point ratio Rdn=1, affiliated layer position Lname=L1 (representing this aspect with the name of the end face of a layer position usually).
The point of incident line S3 to be (P6, P7), starting point Pup=P6, end point Pdn=P7, starting point ratio Rup=0, end point ratio Rdn=0.4 (supposition be 0.4 in the tomography right side ratio of P7 between L1, L2), affiliated layer position Lname=L1.
Another point of incident line S3 to be (P7, P8), starting point Pup=P7, end point Pdn=P8, starting point ratio Rup=0.8 (supposition be 0.8 in the tomography right side ratio of P7 between L1, L2), end point ratio Rdn=1, affiliated layer position Lname=L1.
For the grid of the capable N row of M, there is the incident line Smn position of each node, and the intersection point P of incident line Smn and stratum and tomography goes on record, and the point that per two adjacent intersection point P form also goes on record to PC.The data that strengthen grid format have so just generated.
In step 401,, that is, has the data file that strengthens grid format by the enhancing grid format data file of three-dimensional geological modeling processing module output three-dimensional geological model.
The step of the enhancing grid format data file of output three-dimensional geological model specifically comprises again:
(21), at first the three-dimensional modeling processing module obtains the XY coordinate of grid node according to grid node of choice of location of geological data grid, promptly determines the XY coordinate of an imaginary incident line according to the position of grid node in the geological data; The coordinate of the stratum of calculating the imaginary incident line and its then and being intersected and the intersection point of tomography comprises Z coordinate and XY coordinate;
(22), the information such as title of the coordinate of intersection point, intersection point type and institute's friendship aspect are recorded in certain data structure successively, as chained list or array or other known data structure (by the data structure decision of three-dimensional modeling processing module); Next right per two the neighbouring points that are on this incident line in the research range as point, information such as the right starting point of point, end point, starting point ratio, end point ratio, affiliated layer position are also recorded in the known data structure;
(23), calculate the point at all grid node places successively and put rightly, and record in the data structure; At last the point in the data structure with put right data and be written in the file that strengthens grid format.
In this step, the seismic inversion processing module need be set up one can read the interface that strengthens the grid format data file, strengthens the grid format data file to read; This interface can realize reading in all key points from strengthen the grid format data file and point is right, then these data is carried out the computing of seismic inversion as an input.The step that the interface that strengthens the grid format data file is read in foundation is specially:
(31), name, point that type is identical are summed up as same ground prime element, form same stratal surface or fault surface as all names, point that type is identical; And the information of these points is respectively put in the data structure of these ground prime elements, as the point that belongs to same stratal surface is put in the data structure of this stratal surface; An affiliated layer identical have a few to being summed up as same ground prime element,, is put into the right information of point in the corresponding formation data structure forming same stratum as the identical point of layer under all;
(32), from point with put information such as the position that can obtain the stratum the right information, form, thickness, can obtain the information such as position, form of tomography, can also obtain the information of the mutual alignment relation of stratum and tomography, geologic model just is read into the seismic inversion module and has suffered like this.
Also can from the file of EGF form, read in a little with put right in just point and point to the data structure that is respectively put into the different ground prime element under it in.
In the step 402, name, point that type is identical are summed up as same ground prime element, as being Layer1 with all names, type is that the point of L (stratum) is formed same stratal surface Layer1; All names are F1, and type is that the point of F (tomography) is formed same fault surface F1.Stratal surface can be represented stratum end face or bottom surface, stratum.Stratal surface and fault surface do not have thickness.
In the step 402, the identical have a few in affiliated layer position is summed up as same ground prime element to PC, as be that the point of Layer1 is to composition stratum Layer1 with layer position under all, there is thickness on this stratum, the thickness on each stratum, grid node place is for putting the absolute value that the degree of depth of end point among the PC is deducted the degree of depth of starting point, and geologic model has just generated like this.Any some position and attribute in the space can obtain by the neighbor point interpolation to this point in the geologic model.
See also shown in Figure 5ly, be a specific embodiment of the present invention.This embodiment has described the embodiment of data-switching.Three-dimensional geological model as shown in Figure 5 is a continuous three-dimensional geologic model (geometric model), comprises a tomography F1 and three stratum d2-1, d2-1 and d2-3 in this model, and intersect respectively on tomography F1 and three stratum.The present invention is a research object with stratum d2-1 and d2-2, and on X, Y direction, the present invention gets 7 * 7 grid, and the distance of grid is 250 meters, and on the Z direction, the present invention studies d2-1 and two stratum of d2-2.At first, use the three-dimensional modeling processing module to set up this three-dimensional geological model; Secondly, re-use data file that strengthens grid format of three-dimensional modeling processing module output, comprised the key point information of three-dimensional geological model and point in this file information.The seismic inversion processing module is read in this data file that strengthens grid format, just can obtain three-dimensional three-dimensional geological model, just can carry out the seismic inversion computing in conjunction with other data.
The enhancing grid format data file that generates is as follows:
// file begins
point
{
//No.x y z type Name
0 0 0 800.1 L d2-1
1 0 0 880.6 L d2-2
2 0 0 960.8 L d2-3
3 250 0 815.6 L d2-1
4 250 0 896.2 L d2-2
5 250 0 970.1 L d2-3
6 500 0 826.3 L d2-1
7 500 0 905.6 L d2-2
8 500 0 980.7 L d2-3
9 750 0 830.6 L d2-1
10 750 0 909.2 L d2-2
11 750 0 985.5 L d2-3
12 1000 0 829.5 L d2-1
13 1000 0 908.6 L d2-2
14 1000 0 985.0 L d2-3
15 1250 0 827.7 L d2-1
16 1250 0 907.2 L d2-2
17 1250 0 986.5 L d2-3
18 1500 0 831.1 L d2-1
19 1500 0 913.4 L d2-2
20 1500 0 990.3 L d2-3
21 0 250 810.5 L d2-1
22 0 250 889.3 L d2-2
23 0 250 968.7 L d2-3
24 250 250 824.9 L d2-1
25 250 250 903.7 L d2-2
26 250 250 977.1 L d2-3
27 500 250 833.6 L d2-1
28 500 250 912.6 L d2-2
29 500 250 987.5 L d2-3
30 750 250 836.6 L d2-1
31 750 250 915.8 L d2-2
32 750 250 992.5 L d2-3
33 1000 250 838.1 L d2-1
34 1000 250 915.8 L d2-2
35 1000 250 953.6 F F1
36 1000 250 993.5 L d2-3
37 1250 250 833.6 L d2-1
38 1250 250 912.2 L d2-2
39 1250 250 993.9 L d2-3
40 1500 250 837.6 L d2-1
41 1500 250 919.3 L d2-2
42 1500 250 996.4 L d2-3
43 0 500 813.6 L d2-1
44 0 500 892.4 L d2-2
45 0 500 970.3 L d2-3
46 250 500 827.6 L d2-1
47 250 500 907.2 L d2-2
48 250 500 980.6 L d2-3
49 500 500 836.4 L d2-1
50 500 500 915.8 L d2-2
51 500 500 991.7 L d2-3
52 750 500 839.9 L d2-1
53 750 500 878.9 F F1
54 750 500 918.8 L d2-2
55 750 500 995.9 L d2-3
56 1000 500 843.2 L d2-1
57 1000 500 876.5 F F1
58 1000 500 919.2 L d2-2
59 1000 500 996.8 L d2-3
60 1250 500 835.7 L d2-1
61 1250 500 917.8 L d2-2
62 1250 500 957.2 F F1
63 1250 500 996.5 L d2-3
64 1500 500 843.1 L d2-1
65 1500 500 922.8 L d2-2
66 1500 500 998.1 L d2-3
67 0 750 815.5 L d2-1
68 0 750 891.6 L d2-2
69 0 750 972.4 L d2-3
70 250 750 829.6 L d2-1
71 250 750 906.9 L d2-2
72 250 750 979.9 L d2-3
73 500 750 837.1 L d2-1
74 500 750 916.5 L d2-2
75 500 750 989.8 L d2-3
76 750 750 840.1 L d2-1
77 750 750 919.2 L d2-2
78 750 750 996.6 L d2-3
79 1000 750 843.3 L d2-1
80 1000 750 877.9 F F1
81 1000 750 920.2 L d2-2
82 1000 750 995.3 L d2-3
83 1250 750 836.7 L d2-1
74 1250 750 856.6 F F1
85 1250 750 917.5 L d2-2
86 1250 750 995.5 L d2-3
87 1500 750 848.3 L d2-1
88 1500 750 921.8 L d2-2
89 1500 750 997.6 L d2-3
90 0 1000 813.9 L d2-1
91 0 1000 890.6 L d2-2
92 0 1000 971.5 L d2-3
93 250 1000 828.6 L d2-1
94 250 1000 904.3 L d2-2
95 250 1000 977.6 L d2-3
96 500 1000 835.5 L d2-1
97 500 1000 914.3 L d2-2
98 500 1000 987.2 L d2-3
99 750 1000 838.1 L d2-1
100 750 1000 916.5 L d2-2
101 750 1000 994.3 L d2-3
102 1000 1000 841.9 L d2-1
103 1000 1000 876.4 F F1
104 1000 1000 918.6 L d2-2
105 1000 1000 994.3 L d2-3
106 1250 1000 835.6 L d2-1
107 1250 1000 916.2 L d2-2
108 1250 1000 993.7 L d2-3
109 1500 1000 846.5 L d2-1
110 1500 1000 920.2 L d2-2
111 1500 1000 995.3 L d2-3
112 0 1250 813.3 L d2-1
113 0 1250 889.6 L d2-2
114 0 1250 969.2 L d2-3
115 250 1250 826.1 L d2-1
116 250 1250 902.3 L d2-2
117 250 1250 975.3 L d2-3
118 500 1250 833.6 L d2-1
119 500 1250 912.6 L d2-2
120 500 1250 985.3 L d2-3
121 750 1250 836.2 L d2-1
122 750 1250 913.3 L d2-2
123 750 1250 942.2 L d2-3
124 1000 1250 841.6 L d2-1
125 1000 1250 918.3 L d2-2
126 1000 1250 995.2 L d2-3
127 1250 1250 837.2 L d2-1
128 1250 1250 917.3 L d2-2
129 1250 1250 994.3 L d2-3
130 1500 1250 848.6 L d2-1
131 1500 1250 922.5 L d2-2
132 1500 1250 998.6 L d2-3
133 0 1500 812.6 L d2-1
134 0 1500 887.3 L d2-2
135 0 1500 967.5 L d2-3
136 250 1500 825.1 L d2-1
137 250 1500 901.1 L d2-2
138 250 1500 973.3 L d2-3
139 500 1500 831.2 L d2-1
140 500 1500 912.6 L d2-2
141 500 1500 983.5 L d2-3
142 750 1500 835.9 L d2-1
143 750 1500 913.3 L d2-2
144 750 1500 941.9 L d2-3
145 1000 1500 843.3 L d2-1
146 1000 1500 919.8 L d2-2
147 1000 1500 997.5 L d2-3
148 1250 1500 839.1 L d2-1
149 1250 1500 919.5 L d2-2
150 1250 1500 996.4 L d2-3
151 1500 1500 850.7 L d2-1
152 1500 1500 925.6 L d2-2
153 1500 1500 999.1 L d2-3
}
point?couples
{
//Pup/Pdn/Rup?/Rdn/LayerName
0 1 0 1 d2-1
1 2 0 1 d2-2
3 4 0 1 d2-1
4 5 0 1 d2-2
6 7 0 1 d2-1
7 8 0 1 d2-2
9 10?0 1 d2-1
10 11?0 1 d2-2
12 13?0 1 d2-1
13 14?0 1 d2-2
15 16?0 1 d2-1
16 17?0 1 d2-2
18 19?0 1 d2-1
19 20?0 1 d2-2
21 22?0 1 d2-1
22 23?0 1 d2-2
24 25?0 1 d2-1
25 26?0 1 d2-2
27 28?0 1 d2-1
28 29?0 1 d2-2
30 31 0 1 d2-1
31 32 0 1 d2-2
33 34 0 1 d2-1
34 35 0 0.88 d2-2
35 36 0.88?1 d2-2
37 38 0 1 d2-1
38 39 0 1 d2-2
40 41 0 1 d2-1
41 42 0 1 d2-2
43 44 0 1 d2-1
44 45 0 1 d2-2
46 47 0 1 d2-1
47 48 0 1 d2-2
49 50 0 1 d2-1
50 51 0 1 d2-2
52 53 0 0.53 d2-1
53 54 0.53?1 d2-1
54 55 0 1 d2-2
56 57 0 0.58 d2-1
57 58 0.58?1 d2-1
58 59 0 1 d2-2
60 61 0 1 d2-1
61 62 0 0.89 d2-2
62 63 0.89?1 d2-2
64 65 0 1 d2-1
65 66 0 1 d2-2
67 68 0 1 d2-1
68 69 0 1 d2-2
70 71 0 1 d2-1
71 72 0 1 d2-2
73 74 0 1 d2-1
74 75 0 1 d2-2
76 77 0 1 d2-1
77 78 0 1 d2-2
79 80 0 0.58 d2-1
80 81 0.58?1 d2-1
81 82 0 1 d2-2
83 84 0 0.57 d2-1
84 85 0.57?1 d2-1
85 86 0 1 d2-2
87 88 0 1 d2-1
88 89 0 1 d2-2
90 91 0 1 d2-1
91 92 0 1 d2-2
93 94 0 1 d2-1
94 95 0 1 d2-2
96 97 0 1 d2-1
97 98 0 1 d2-2
99 100?0 1 d2-1
100 101?0 1 d2-2
102 103?0 0.58 d2-1
103 104?0.58?1 d2-1
104 105?0 1 d2-2
106 107?0 1 d2-1
107 108?0 1 d2-2
109 110?0 1 d2-1
110 111?0 1 d2-2
112 113?0 1 d2-1
113 114?0 1 d2-2
115 116?0 1 d2-1
116 117 0 1?d2-2
118 119 0 1?d2-1
119 120 0 1?d2-2
121 122 0 1?d2-1
122 123 0 1?d2-2
124 125 0 1?d2-1
125 126 0 1?d2-2
127 128 0 1?d2-1
128 129 0 1?d2-2
130 131 0 1?d2-1
131 132 0 1?d2-2
133 134 0 1?d2-1
134 135 0 1?d2-2
136 137 0 1?d2-1
137 138 0 1?d2-2
139 140 0 1?d2-1
140 141 0 1?d2-2
142 143 0 1?d2-1
143 144 0 1?d2-2
145 146 0 1?d2-1
146 147 0 1?d2-2
148 149 0 1?d2-1
149 150 0 1?d2-2
151 152 0 1?d2-1
152 153 0 1?d2-2
} // the end of file
See also shown in Figure 6, for the present invention asks in the three-dimensional geological model process flow diagram of some property values arbitrarily; This flow process comprises the steps:
Step 601 determines according to the XYZ coordinate of arbitrfary point which grid node the arbitrfary point is in, which point between, and determine that further the arbitrfary point is in which stratum and the ratio in the stratum;
Step 602 calculates a new aspect according to described ratio and layer of living in position; And
Step 603 is asked for the intersection point of logging trace track and this new aspect, obtains the property value of described arbitrfary point by interpolation according to the property value of described intersection point.
See also shown in Figure 7ly, be interface data conversion equipment synoptic diagram of the present invention.This interface data conversion equipment 70 comprises three-dimensional modeling processing module 71, seismic inversion processing module 72, and three-dimensional modeling processing module 71 comprises interface data modular converter 711, interface data output module 712 again; Wherein interface data modular converter 711 is used to generate the enhancing grid format data file of three-dimensional geological model, and interface data output module 712 is used to export the enhancing grid format data file of three-dimensional geological model; Seismic inversion processing module 72 is used to read in enhancing grid format data file, and carries out the seismic inversion computing according to strengthening the grid format data file.
Enhancing grid format of the present invention very accurate description contains the three-dimensional geological model of trap-up, and adds a layer in containing the three-dimensional geological model of trap-up arbitrarily.
The present invention proposes a kind of interface data conversion method, realized the conversion of data between simple, a convenient three-dimensional modeling processing module (software) that exchanges and the seismic inversion processing module (software).Related enhancing grid format has comprised the structural information of necessity of three-dimensional geological model, discrete data structure adapts with the data layout of seismic inversion processing module again, this form can be a kind of simple, the discrete data layout that can express the three-dimensional geological model topological structure again exactly with the three-dimensional geological model of simple structrual description complexity; By the function served as bridge of this enhancing grid format, can make geology three-dimensional modeling processing module and the seismic inversion processing module is extraordinary combines, make things convenient for the earthquake inversion processing module to use.
Certainly; the present invention also can have other various embodiments; under the situation that does not deviate from spirit of the present invention and essence thereof; those of ordinary skill in the art work as can make various corresponding changes and distortion according to the present invention, but these corresponding changes and distortion all should belong to the protection domain of the appended claim of the present invention.
Claims (11)
1. an interface data conversion method is applicable to the conversion of data between three-dimensional modeling processing module and the seismic inversion processing module, it is characterized in that, comprises the steps:
Step 1 is by the enhancing grid format data file of described three-dimensional modeling processing module generation and output three-dimensional geological model;
Step 2 is read in described enhancing grid format data file by described seismic inversion processing module, and carries out the seismic inversion computing according to described enhancing grid format data file;
Described enhancing grid format data file has on X, Y direction and geological data unanimity, discrete data structure that Gate lattice number is identical with geological data, and be included in key point and the right structural information of point that continuous recording on the Z direction has three-dimensional geological model, this key point is the incident line and the stratal surface of three-dimensional geological model and the intersection point of fault surface on the imaginary Z direction in grid node place, and it is right that per two the neighbouring key points in same grid node place are combined into a point;
Wherein, step 2 further comprises: name, key point that type is identical are summed up as same ground prime element, and the identical have a few in affiliated layer position is to being summed up as same ground prime element.
2. interface data conversion method according to claim 1, it is characterized in that, described step 1 specifically comprises again: described three-dimensional modeling processing module is set up a three-dimensional geological model, mark off the grid consistent on the XY face in this geologic model, penetrate all stratum and tomography to the incident line at each imagination Z of grid node place with geological data; The intersection point of described every incident line and stratal surface and fault surface is a key point; And with in this three-dimensional geological model at same grid node place per two neighbouring key points be combined into a little right, note all key points, point to and be stored in the file that strengthens grid format, generate and strengthen the grid format data file.
3. interface data conversion method according to claim 1, it is characterized in that, described step 2 is specially: described seismic inversion processing module is read in the data and the right data of point of the key point of three-dimensional geological model from described enhancing grid format data file, and the data of described key point are carried out the seismic inversion computing with the right data of point as input.
4. interface data conversion method according to claim 1 is characterized in that, described step 2 further comprises: the absolute value that the degree of depth of described some centering end point deducts the degree of depth of starting point is the thickness on each stratum, grid node place.
5. interface data conversion method according to claim 3 is characterized in that, described step 2 further comprises: the step of asking for the property value of arbitrfary point in the three-dimensional geological model by described enhancing grid format data file in conjunction with the logging trace data.
6. interface data conversion method according to claim 5 is characterized in that, the described step of asking for the property value of arbitrfary point in the three-dimensional geological model specifically comprises again:
Step 701 determines that according to the XYZ coordinate of this arbitrfary point its grid node of living in, point of living in are right, and further determines stratum of living in, this arbitrfary point and the ratio in stratum of living in;
Step 702 calculates a new aspect according to described ratio and stratum of living in; And
Step 703 is asked for the intersection point of logging trace track and this new aspect, obtains the property value of this arbitrfary point by interpolation according to the property value of described intersection point.
7. interface data conversion method according to claim 3, it is characterized in that, in the described step 2, also comprise by described seismic inversion processing module and asking in the stratum and end face and the step of bottom surface apart from the calculating aspect that becomes fixed proportion according to described enhancing grid format data file.
8. interface data conversion method according to claim 3 is characterized in that, in the described step 2, also comprises by described seismic inversion processing module and asks in the stratum step with the equidistant calculating aspect of end face according to described enhancing grid format data file.
9. interface data conversion method according to claim 3 is characterized in that, in the described step 2, also comprises by described seismic inversion processing module and asks in the stratum step with the equidistant calculating aspect in bottom surface according to described enhancing grid format data file.
10. interface data conversion equipment, it is characterized in that, comprise a three-dimensional modeling processing module, a seismic inversion processing module that is connected with described three-dimensional modeling processing module, described three-dimensional modeling processing module comprises an interface data modular converter, an interface data output module again;
Described interface data modular converter is used to generate the enhancing grid format data file of a three-dimensional geological model; Described interface data output module is used to export described enhancing grid format data file; Described seismic inversion processing module is used to read in described enhancing grid format data file, and carries out the seismic inversion computing according to described enhancing grid format data file;
Described enhancing grid format data file has on X, Y direction and geological data unanimity, discrete data structure that Gate lattice number is identical with geological data, and be included in key point and the right structural information of point that continuous recording on the Z direction has three-dimensional geological model, this key point is the incident line and the stratal surface of three-dimensional geological model and the intersection point of fault surface on the imaginary Z direction in grid node place, and it is right that per two the neighbouring key points in same grid node place are combined into a point;
Wherein, described seismic inversion processing module is summed up as same ground prime element with name, key point that type is identical, and the identical have a few in affiliated layer position is to being summed up as same ground prime element.
11. interface data conversion equipment according to claim 10 is characterized in that, described seismic inversion processing module is used to read the interface that strengthens Gate sound of laughing formula data file by one and reads in described enhancing grid format file.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2006100114915A CN101038677B (en) | 2006-03-14 | 2006-03-14 | Interface data conversion method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2006100114915A CN101038677B (en) | 2006-03-14 | 2006-03-14 | Interface data conversion method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101038677A CN101038677A (en) | 2007-09-19 |
| CN101038677B true CN101038677B (en) | 2011-09-21 |
Family
ID=38889559
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2006100114915A Active CN101038677B (en) | 2006-03-14 | 2006-03-14 | Interface data conversion method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101038677B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101303414B (en) * | 2008-05-22 | 2011-03-30 | 北京航空航天大学 | Method for generating ground layer surface and geologic body based on level set |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1316722A (en) * | 2000-04-03 | 2001-10-10 | 上海大唐移动通信设备有限公司 | Method for topographic contour line conversion |
| CN1321249A (en) * | 1998-09-02 | 2001-11-07 | 菲利浦石油公司 | Automated seismic fault detection and picking |
| US20020013687A1 (en) * | 2000-03-27 | 2002-01-31 | Ortoleva Peter J. | Methods and systems for simulation-enhanced fracture detections in sedimentary basins |
| US20020120429A1 (en) * | 2000-12-08 | 2002-08-29 | Peter Ortoleva | Methods for modeling multi-dimensional domains using information theory to resolve gaps in data and in theories |
| WO2005076226A1 (en) * | 2004-02-03 | 2005-08-18 | British Telecommunications Public Limited Company | Animation techniques using quaternions |
-
2006
- 2006-03-14 CN CN2006100114915A patent/CN101038677B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1321249A (en) * | 1998-09-02 | 2001-11-07 | 菲利浦石油公司 | Automated seismic fault detection and picking |
| US20020013687A1 (en) * | 2000-03-27 | 2002-01-31 | Ortoleva Peter J. | Methods and systems for simulation-enhanced fracture detections in sedimentary basins |
| CN1316722A (en) * | 2000-04-03 | 2001-10-10 | 上海大唐移动通信设备有限公司 | Method for topographic contour line conversion |
| US20020120429A1 (en) * | 2000-12-08 | 2002-08-29 | Peter Ortoleva | Methods for modeling multi-dimensional domains using information theory to resolve gaps in data and in theories |
| WO2005076226A1 (en) * | 2004-02-03 | 2005-08-18 | British Telecommunications Public Limited Company | Animation techniques using quaternions |
Non-Patent Citations (4)
| Title |
|---|
| 陈立潮.基于栅格的GIS三维空间数据模型.计算机工程30 8.2004,30(8),4-6. |
| 陈立潮.基于栅格的GIS三维空间数据模型.计算机工程30 8.2004,30(8),4-6. * |
| 龚健雅.矢量与栅格集成的三维数据模型.武汉测绘科技大学学报22 1.1997,22(1),7-15. |
| 龚健雅.矢量与栅格集成的三维数据模型.武汉测绘科技大学学报22 1.1997,22(1),7-15. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101038677A (en) | 2007-09-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Dell’Unto et al. | Experiencing ancient buildings from a 3D GIS perspective: a case drawn from the Swedish Pompeii Project | |
| CN101882171B (en) | Method for fast establishing interactive tunnel and wall rock body three-dimensional models | |
| Saygi et al. | Evaluation of GIS and BIM roles for the information management of historical buildings | |
| CN102612682B (en) | For the method and apparatus of reservoir modeling and simulation | |
| Jusuf et al. | Integrated modeling of CityGML and IFC for city/neighborhood development for urban microclimates analysis | |
| CN110210120B (en) | Automatic modeling method for railway engineering based on parameterization | |
| CN102129715B (en) | Generation method of quadrilateral grid of geometric model with any internal feature constraints | |
| Sprague et al. | Interpretive tools for 3-D structural geological modelling part II: Surface design from sparse spatial data | |
| CN103236088B (en) | Three-dimensional geological structure modeling method and system based on composite grid | |
| CN108256218B (en) | Underground communication roadway fine modeling method based on actually measured waist line data | |
| CN109101732B (en) | Classification-free river channel two-dimensional structure grid subdivision method based on topographic feature boundary line | |
| CN111754626B (en) | Topology-based borehole geological profile modeling method | |
| CN116152461B (en) | Geological modeling method, device, computer equipment and computer readable storage medium | |
| CN107315847B (en) | A kind of generation method and device in river and underground water coupled simulation parameter | |
| CN110263434A (en) | A kind of flow unit method for numerical simulation based on multiple dimensioned mixed finite element | |
| CN108153917A (en) | A kind of real estate investigation banking process for visualizing with defining people room relationship | |
| CN101038677B (en) | Interface data conversion method and device | |
| CN113190910B (en) | Slice analysis system and method based on BIM three-dimensional geological tunnel model | |
| CN106649776A (en) | Method of semi-automating comprehensive vector polygon | |
| CN107886573B (en) | Slope three-dimensional finite element grid generation method under complex geological conditions | |
| CN113779663A (en) | BIM-based subway station enclosure structure three-dimensional modeling method, system and medium | |
| CN105719342A (en) | Three-dimensional modeling visualization method and device of ground fracture geologic body | |
| CN109858160A (en) | A kind of modeling method of the rail traffic geological information model based on BIM technology | |
| CN103714494B (en) | Power circuit method for building up and device | |
| JP2007264952A (en) | Ground-analyzing mesh generation method and ground-analyzing mesh generation program |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C53 | Correction of patent for invention or patent application | ||
| CB03 | Change of inventor or designer information |
Inventor after: Meng Xianhai Inventor after: Yang Qin Inventor after: Li Jigang Inventor before: Yang Qin |
|
| COR | Change of bibliographic data |
Free format text: CORRECT: INVENTOR; FROM: YANG QIN TO: MENG XIANHAI YANG QIN LI JIGANG |
|
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20170602 Address after: 100088 Beijing City, Haidian District Zhichun Road No. 6 (Jinqiu International Building) A block, room 1501 Patentee after: Beijing grid world software technology Limited by Share Ltd Address before: 100083 Haidian District, Xueyuan Road, No. 37, Patentee before: Beihang University |