CN103344297B - The reservoir storage survey method of underground water seal cave depot - Google Patents

The reservoir storage survey method of underground water seal cave depot Download PDF

Info

Publication number
CN103344297B
CN103344297B CN201310271053.2A CN201310271053A CN103344297B CN 103344297 B CN103344297 B CN 103344297B CN 201310271053 A CN201310271053 A CN 201310271053A CN 103344297 B CN103344297 B CN 103344297B
Authority
CN
China
Prior art keywords
section
cave depot
water seal
underground water
point
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
Application number
CN201310271053.2A
Other languages
Chinese (zh)
Other versions
CN103344297A (en
Inventor
吕宏权
刘斌
何文浩
付志伦
徐学斌
张世飞
牛春霞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Railway Tunnel Group Yichu Co Ltd
Original Assignee
China Railway Tunnel Group Yichu Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Railway Tunnel Group Yichu Co Ltd filed Critical China Railway Tunnel Group Yichu Co Ltd
Priority to CN201310271053.2A priority Critical patent/CN103344297B/en
Publication of CN103344297A publication Critical patent/CN103344297A/en
Application granted granted Critical
Publication of CN103344297B publication Critical patent/CN103344297B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Length Measuring Devices With Unspecified Measuring Means (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The present invention proposes a kind of reservoir storage survey method of underground water seal cave depot, comprise the steps: step 1, set up underground water seal cave depot SURVEYING CONTROL NETWORK; Step 2, after having set up SURVEYING CONTROL NETWORK, definition underground water seal cave depot section survey parameter, scans underground water seal cave depot; Data after scanning are arranged, finally obtain the three-dimensional data of underground water seal cave depot by step 3; Step 4, carries out three-dimensional modeling to the three-dimensional data of underground water seal cave depot shaping, namely sets up some adjacent triangulation networks, forms solid model, then carries out according to the axle vertical direction of definition the section that cutting forms corresponding height; Step 5, calculate the area of section, with the height value of corresponding equalization, obtain the volume of described section, according to the overall height of the underground water seal cave depot measured, obtain underground water seal cave depot aggregate storage capacity, the positional accuracy of adjustment underground water seal cave depot section survey and some bit density obtain the accurate storage capacity volume of underground water seal cave depot.

Description

The reservoir storage survey method of underground water seal cave depot
Technical field
The present invention relates to cave depot storage capacity cubic content measurement field, particularly relate to a kind of reservoir storage survey method of underground water seal cave depot.
Background technology
Society, volume for underground karst cavity or cave depot is measured accurately and is belonged to technical matters urgently to be resolved hurrily, although current 3-D scanning technology has larger achievement in the collection to contour of object three-dimensional coordinate, but because the span of cave depot is large, high side wall, irregular rock wall, all bring very large technical barrier to measurement volumes, the collection of the larger cave depot volume of very difficult realization, even if measure corresponding volume due to its condition of scanning harshness to be also inaccurate, and error is very large, and instrument and equipment price costly, be widely used and also have some limitations.Measurement aspect for cave depot capacity does not use high-tech product to realize corresponding cave depot measurement means at present yet.
Summary of the invention
The present invention is intended at least solve the technical matters existed in prior art, especially innovatively proposes a kind of reservoir storage survey method of underground water seal cave depot, comprises the steps:
Step 1, sets up underground water seal cave depot SURVEYING CONTROL NETWORK;
Step 2, after having set up SURVEYING CONTROL NETWORK, definition underground water seal cave depot section survey parameter, scans underground water seal cave depot;
Data after scanning are arranged, finally obtain the three-dimensional data of underground water seal cave depot by step 3;
Step 4, carries out three-dimensional modeling to the three-dimensional data of underground water seal cave depot shaping, namely sets up some adjacent triangulation networks, forms solid model, then carries out according to the axle vertical direction of definition the section that cutting forms corresponding height;
Step 5, calculate the area of section, with the height value of corresponding equalization, obtain the volume of described section, according to the overall height of the underground water seal cave depot measured, obtain underground water seal cave depot aggregate storage capacity, the positional accuracy of adjustment underground water seal cave depot section survey and some bit density obtain the accurate storage capacity volume of underground water seal cave depot.
The reservoir storage survey method of described underground water seal cave depot, preferably, described step 1 comprises:
Step 1-1, ground control point or in-cavity construction reference mark are that benchmark sets up unified SURVEYING CONTROL NETWORK to whole cave depot, to the encryption needing the position measured to carry out reference mark on the basis of Controling network, adopt accurate fourth class tracerse survey method, after compensating computation, obtain the three-dimensional coordinate at each reference mark.
The reservoir storage survey method of described underground water seal cave depot, preferably, described step 2 comprises:
Step 2-1, after having set up underground water seal cave depot SURVEYING CONTROL NETWORK, according to the design transversal section parameter of design data definition vertical curve, horizontal curve and each measuring point, obtains the design reference section closed;
Step 2-2, definition vertical curve, wherein the data of vertical curve are mileage, elevation, the radius of knick point and tangent line, Origin And Destination;
Step 2-3, definition horizontal curve, wherein horizontal curve is the data for describing, determining center line of road accurate location, and horizontal curve data are mileage, line style, turning radius, the coordinate of line element starting point;
Step 2-4, definition transversal section parameter, wherein transversal section supplemental characteristic is that right-handed coordinate system is tried to achieve, coordinate system with circuit at section mileage place designed elevation point for initial point, with normal on the right side of line midline working direction for X-axis, with initial point zenith direction for Y-axis, during In-put design section, according to inputting every section of straight line and circular arc successively clockwise.
The reservoir storage survey method of described underground water seal cave depot, preferably, described step 3 comprises:
Step 3-1, total powerstation is set up near center line and centering leveling bottom underground water seal cave depot, adopt the method for Free Station, obtain the three-dimensional coordinate of survey station point, select, horizontal curve, vertical curve, transversal section type, to arrange according to the dot spacing of 1 meter, mileage direction and vertical direction 0.5 meter, then carry out measurements and scan;
Step 3-2, record transversal section mileage, the period of measurement point, partially middle distance, the discrepancy in elevation, three-dimensional coordinate.
The reservoir storage survey method of described underground water seal cave depot, preferably, described step 4 comprises:
Step 4-1, adopt measured database to carry out manipulate measurement data, cross-sectioning data are converted to line file form, each the section point first section of being combined in line file, all sections by a cave depot form line string;
Step 4-2, by some adjacent triangular facets, is rolled into the entity that inside and outside is ventilative; Entity is by a series of point on line, and be linked to be the triangulation network that inside and outside is ventilative, must there be adjacent triangular facet on the limit of any triangular facet, 3 summits of any triangular facet, must be attached on effective point;
Step 4-3, if each triangular facet of composition entity exists self intersection, without adjacent limit, repeats limit, invalid edges, then entity is not effective entity, and invalid entity can not calculate volume, space constraint, logical operation;
Step 4-4, by the validation verification of entity, between entities, the computing of the common factor in space, union and difference set is carried out between entity and face, definition is using Z axis as cross-sectional axis, cutting is carried out with the vertical direction of axle, definition storage capacity measuring point elevation, input section absolute altitude, automatically will calculate and draw out the section of this absolute altitude.
The reservoir storage survey method of described underground water seal cave depot, preferably, described step 4-2 also comprises the step creating the triangulation network:
Step 4-5, between two sections: connect the triangulation network between different sections;
Step 4-6, in a section: the connection automatically completing the triangulation network in a section chosen;
Step 4-7, section is to a point: a section is to a some connection triangulation network;
Step 4-8, between multiple sections: automatically connect the triangulation network between a series of line and section;
Step 4-9, uses control line: define section with control line, and then produces the triangulation network;
Step 4-10, according to manual reconnaissance: the manual definition triangulation network retrains, and connects the triangulation network between section;
Step 4-11, single triangle: define each triangle by selected triangular apex;
Step 4-12, from a section to two sections: a closed female section and two sub-sections or point, use associating (Union) concept, control branch line, and then produce the triangulation network;
Step 4-13, from a section to multiple sections: connect the triangulation network a female section and multiple sub-section;
Step 4-14, uses center line and section: utilize section to connect the triangulation network along center line, generally generates underworkings entity by this function.
The reservoir storage survey method of described underground water seal cave depot, preferably, described step 5 comprises:
Step 5-1, calculates all horizontal or vertical cross-sectional areas,
wherein V 1for initial computing volume, S 0be the first cross sectional area, S 1be the area of the second transversal section, L 1be the first computing length,
Aggregate storage capacity is tried to achieve again with the length integration of each measuring unit,
V = Σ i = 1 n ( Si × Li ) , V is aggregate storage capacity volume, and i, n are positive integer;
Step 5-2, after calculating underground water seal cave depot aggregate storage capacity, adjustment positional accuracy and some bit density;
Step 5-3, carries out the positional accuracy of section survey and puts bit density to the impact analysis of aggregate storage capacity; During section survey, each data point that section is upsampled to is the most elementary cell forming storage capacity, and the set of point determines underground water seal cave depot aggregate storage capacity, section is put the area measurement precision that the precision of position and Effects of Density section; Collect each point coordinate by sampled point is carried out angle measurement, survey limit calculate;
Step 5-4, underground water seal cave depot section is with city gate type and three-core round section composition, and wherein type section in city gate is mainly distributed in and connects tunnel, if main cave depot take radius as the circle cross section of R, then has:
Area S=π R2,
If area measurement medial error is m s, radius measurement medial error is m r, according to law of propagation of errors, can obtain:
wherein subscript S is for measuring area, subscript R for measuring radius;
If cross-sectional area measuring accuracy is 1/10000, then the cave depot of 20m radius, it surveys edge error can not be greater than ± 1mm; Equally, if the direction mean square error of angle observation of instrument is m β, then have:
wherein subscript S measures the direction angle measurement that area, subscript β are instrument;
When direction, mean square error of angle observation is less, also less on the impact of cross-sectional area measuring accuracy;
Step 5-5, some bit density is the spacing distance between points of sampling, to the angle step value of sampled point in section survey; Because hole wall surface is coarse uneven, sampled point is exactly measure on coarse surface, and concavo-convex difference is not of uniform size; There is representive error problem for section sampling, concerning whole cave depot, concavo-convex difference meets normal distribution law, just can be solved the representive error of cave depot appearance like this by the density of constraint fixed sample point position; The maximum radius of underground cavern section is 26 meters, minimum be 4 meters, actual when gathering section point by radius size set point between arc length be 0.2--0.5m, corresponding angle step value is 30-10; Under the condition of same measuring accuracy and sample density, cross sectional radii is larger, and the angle step value of sampled point can suitably amplify; The density of point increases, and curve between points just can shorten as far as possible;
Step 5-6, hold manually to calculate cave depot, volume folded between its two parallel cutting surfaces in conventional section is pressed trapezoid cylinder method and is calculated.By differentiate and can obtain medial error m v, wherein subscript V is for measuring aggregate storage capacity volume:
m V = | ΔL 2 ( S 0 2 + S 1 2 ) 4 · ( m s S ) 2 + ( S 0 + S 1 ) 2 4 · ΔL 2 · ( m ΔL ΔL ) 2 | , Wherein subscript S is for measuring area,
Above formula pushes away the phase error of centralization of cave depot volume is:
m V V = | | ( m ΔL L ) 2 + S 0 2 + S 1 2 ( S 0 + S 1 ) 2 ( m S S ) 2 | | , Wherein subscript S is for measuring area,
The phase error of centralization of cave depot volume, mainly by the phase error of centralization of cave depot length L direction increment and the impact of the cross-sectional area phase error of centralization, for ensureing the positional accuracy of cave depot aggregate storage capacity, needs Measure section point bit density.
In sum, owing to have employed technique scheme, the invention has the beneficial effects as follows:
At present, the method solving underground water seal cave depot reservoir storage survey problem, except above-mentioned section survey method, also has one to be exactly utilize three-dimensional laser scanner to carry out three-dimensional measurement, then resolves through set.But this instrument belongs to new things, expensive, being widely used also has some limitations.Underground water seal cave depot belongs to construction frontier, and section survey method is relatively convenient, easy to operate, and by analyzing, the phase error of centralization of cross-sectional area also easily ensures.
Additional aspect of the present invention and advantage will part provide in the following description, and part will become obvious from the following description, or be recognized by practice of the present invention.
Accompanying drawing explanation
Above-mentioned and/or additional aspect of the present invention and advantage will become obvious and easy understand from accompanying drawing below combining to the description of embodiment, wherein:
Fig. 1 is the reservoir storage survey method flow diagram of underground water seal cave depot of the present invention;
Fig. 2 is the reservoir storage survey method Design of Vertical Curve figure of underground water seal cave depot of the present invention;
Fig. 3 is the reservoir storage survey method Section Design figure of underground water seal cave depot of the present invention;
Fig. 4 is the free instrumentation plan of reservoir storage survey method of underground water seal cave depot of the present invention;
Fig. 5 is the reservoir storage survey method autoscan schematic diagram of underground water seal cave depot of the present invention;
Fig. 6 is the reservoir storage survey method cave depot sectional schematic diagram of underground water seal cave depot of the present invention;
Fig. 7 is that the reservoir storage survey method of underground water seal cave depot of the present invention creates triangulation network schematic diagram;
Fig. 8 is the reservoir storage survey method three-dimensional modeling schematic diagram of underground water seal cave depot of the present invention;
Fig. 9 is that the reservoir storage survey method entity of underground water seal cave depot of the present invention calculates schematic diagram;
Figure 10 is that the reservoir storage survey method of underground water seal cave depot of the present invention calculates aggregate storage capacity simulation drawing.
Embodiment
Be described below in detail embodiments of the invention, the example of described embodiment is shown in the drawings, and wherein same or similar label represents same or similar element or has element that is identical or similar functions from start to finish.Being exemplary below by the embodiment be described with reference to the drawings, only for explaining the present invention, and can not limitation of the present invention being interpreted as.
In describing the invention, it will be appreciated that, term " longitudinal direction ", " transverse direction ", " on ", D score, "front", "rear", "left", "right", " vertically ", " level ", " top ", " end " " interior ", the orientation of the instruction such as " outward " or position relationship be based on orientation shown in the drawings or position relationship, only the present invention for convenience of description and simplified characterization, instead of indicate or imply that the device of indication or element must have specific orientation, with specific azimuth configuration and operation, therefore can not be interpreted as limitation of the present invention.
In describing the invention, unless otherwise prescribed and limit, it should be noted that, term " installation ", " being connected ", " connection " should be interpreted broadly, such as, can be mechanical connection or electrical connection, also can be the connection of two element internals, can be directly be connected, also indirectly can be connected by intermediary, for the ordinary skill in the art, the concrete meaning of above-mentioned term can be understood as the case may be.
The hole reservoir storage survey of underground storage is exactly carry out volumetric measurement to the shape of hole tank, and cave depot volumetric measurement position comprises main cavern, pump hole, water collecting basin, and enters connection tunnel in the tank of hole, construction tunnel, vertical shaft.Because the construction of underground water seal cave depot does not also have corresponding job specifications, hole reservoir storage survey does not have ripe measuring method and correlation technique data yet, the present invention has inquired into the Combination application of high-precision automatic measuring total powerstation, tunnel cross-section measurement system and d solid modeling, according to the section point set data gathered, create cave depot three-dimensional model, carry out the cutting of horizontal direction according to the cave depot of modeling faced by selected reference level, obtain certain depth volumometer scale, cutting slice map and the cave depot three-dimensional model section reserve report that regulation requires.
As shown in Figure 1, the present invention discloses a kind of reservoir storage survey method of underground water seal cave depot, comprises the steps:
Step 1, sets up underground water seal cave depot SURVEYING CONTROL NETWORK;
Step 2, after having set up SURVEYING CONTROL NETWORK, definition underground water seal cave depot section survey parameter, scans underground water seal cave depot;
Data after scanning are arranged, finally obtain the three-dimensional data of underground water seal cave depot by step 3;
Step 4, carries out three-dimensional modeling to the three-dimensional data of underground water seal cave depot shaping, namely sets up some adjacent triangulation networks, forms solid model, then carries out according to the axle vertical direction of definition the section that cutting forms corresponding height;
Step 5, calculate the area of section, with the height value of corresponding equalization, obtain the volume of described section, according to the overall height of the underground water seal cave depot measured, obtain underground water seal cave depot aggregate storage capacity, the positional accuracy of adjustment underground water seal cave depot section survey and some bit density obtain the accurate storage capacity volume of underground water seal cave depot.
1 hole reservoir storage survey
1.1 reservoir storage survey principles
Cave depot volumetric measurement adopts non-prism intelligent total powerstation+airborne section survey, gathers horizontal or vertical section along cave depot direction according to certain spacing.The section line of surveying according to each section calculates the volume that cross-sectional area and section spacing calculate each cave depot, then to the 3 d measurement data gathered, carry out the foundation of cave depot Spatial three-dimensional dispersion model, make the mould shapes produced can be used for three-dimensional visualization, thus the cave depot volume calculated under a certain liquid level elevation and arbitrary orientation cutting profile, obtain physical capacity report.
1.2 main measuring equipments
The collection of three-dimensional coordinate data, the process of measurement data, the foundation etc. of model is completed in cave depot measurement.
Table 1 measuring equipment complete list
Table1measuringequipmentsoftwarelist
Device name Model Precision Remarks
Total powerstation TM/TS/TCRA or other I level total powerstation Non-prism, intelligent
Airborne fractograph analysis Airborne tunnel cross-section Meet measuring requirement Supporting
D solid modeling Three-dimensional artificial modeling Meet measuring requirement
Prism ----- ----- Supporting
Tripod Total station foot rest ----- Supporting
Hygrothermograph ----- ±1.5℃、7%RH
Strong-light torch ----- ----- Low-voltage
Platform truck 25m height semi-girder ----- 25m-ascends a height equipment
Securing band ----- Meet the demands-----
Reflection plaster ----- Meet the demands-----
1.3 cave depots hold measuring method step
Employing two precision in reservoir storage survey is 0.5 ", the come card TM30 of 1mm+1ppm and TS30 precision measurement robot, have employed supporting total powerstation positioning system (TPS) Section measuring system.Before carrying out reservoir storage survey, the impurity high pressure water flushing such as dust, chip ballast in hole tank and all measuring points should be cleaned down and keeps Guan Jige measuring point, hole base plate without ponding.
1.3.1 set up cave depot SURVEYING CONTROL NETWORK
Unified SURVEYING CONTROL NETWORK is set up to whole cave depot for benchmark in the ground control point handed over design institute or in-cavity construction reference mark, just having started to set up measurement net is in order to the unified coordinate system when image data, make image data in the same coordinate system, to the encryption needing the position measured to carry out reference mark on the basis of Controling network, measuring method adopts accurate fourth class tracerse survey, obtains the three-dimensional coordinate at each reference mark after compensating computation.Wherein precise traversing precision is not more than 2.5 according to mean square error of angle observation ", relative closing error of extreme length is not more than 1/35000, the relative mean square error of a point of consecutive point is not more than 8mm and carries out; The measurement of higher degree is according to the fourth class, and precision is not more than every kilometer of full medial error MW=± 10mm of the discrepancy in elevation by the full medial error of the discrepancy in elevation and carries out testing.
1.3.2 definition of T PS tunnel cross-section measurement systematic parameter
After the design parameter of cave depot is checked, according to the design section parameter of design data definition horizontal curve, vertical curve and each measuring point in the Section measuring system of total powerstation, obtain the design reference section closed.
(1) vertical curve typing
The typing work of vertical curve is relatively easy, comparatively fast can complete typing and work is checked in comparison.Wherein note clear and definite: large mileage direction, knick point; Input will have: mileage, elevation, radius and tangent line, and wherein radius divides positive and negative; Not radius and tangent line, the Origin And Destination of knick point.Be illustrated in figure 2 the reservoir storage survey method Design of Vertical Curve figure of underground water seal cave depot.
(2) horizontal curve typing
Plane alignment refers to one group of data for describing, determining center line of road accurate location, comprises the data such as the mileage of line element starting point, line style, turning radius, coordinate.The typing of plane alignment elements mainly utilizes line element method to carry out circuit and checks, can embody the calculation review advantage of line element method preferably.
The typing of transversal section also may be used for checking design drawing parameter, it should be noted that the typing form corresponding to different metering systems is slightly distinguished, need conscientiously to check drawing, guarantee errorless after logging data is imported airborne measurement device, carry out checking and be convenient to further field data collecting work.
(3) design transversal section parameter typing
Section coordinate system: transversal section curve element data are tried to achieve according to right-handed coordinate system, coordinate system with circuit at section mileage place designed elevation point for initial point, with normal on the right side of line midline working direction for X-axis, with initial point zenith direction for Y-axis, during In-put design section, according to inputting every section of straight line and circular arc successively clockwise.Wherein note: the determination of coordinate system and initial point, large mileage direction, clockwise direction; Input key element: the coordinate of each line style change place, line style, radius, central angle; Section line style is divided: straight line and circle, section must be closed.Be illustrated in figure 3 the reservoir storage survey method Section Design figure of underground water seal cave depot.
1.3.3 gather section three-dimensional data
Total powerstation is set up and centering leveling near cave depot (tunnel) base plate center line, adopt the method for Free Station, according to the mode of Fig. 4, instrument hypothesis is at a0 point, enter Free Station with Total Station electronic Thacheometer method, respectively A point and B point being measured according to prompting, obtain the three-dimensional coordinate of survey station point, during through setting precision of standing meet the demands and use the 3rd point to check errorless, total powerstation can be put as without prism pattern, select, horizontal curve, vertical curve, the design parameters such as transversal section type, dot spacing and section spacing.Arrange according to the dot spacing of 1 meter, mileage direction with vertical direction 0.5 meter, then automatically start to measure scanning, see Fig. 5, simultaneously and record section mileage, the period of measurement point, partially in distance, the discrepancy in elevation, three-dimensional coordinate.
The data of cave depot end socket and crossing section gather according to the technique of painting pointwise encryption of intersection
As shown in Figure 6, each section with " >> " boundary, after immediately following this section mileage and designed elevation.Next behavior measurement point information again, respectively with period, partially middle distance, the discrepancy in elevation, X-coordinate, Y-coordinate, the arrangement of H elevation.Bring great convenience to the process of follow-up data territory.
1.3.4 the modeling of cave depot and cutting
(1) the pre-service of measurement data
Employing measured database carrys out manipulate measurement data, first cross-sectioning data are converted to line file form, understand the data layout of lower line file first: each the section point first section of being combined in line file, all sections by a cave depot form line string, process is difficult to for our this measuring method, understood the form of line file, change raw measurement data form now, source code is as follows:
Sub Data Format Transform ()
DimnAsInteger
DimText()asdouble
DimFileNameAsString
DimHZAsVariant
FileName=Application.GetOpenFilename (" storage capacity profileometry data (* .OBS), * .OBS ")
OpenFileNameForInputAs#1
DoWhileNotEOF(1)
n=n+1
ReDimPreserveText(1Ton)
LineInput#1,Text(n)
Loop
Close#1
OpenApp.Path+ " storage capacity .str " ForOutputAs#1
Print#1,"aaa,20-Jun-12,,ssi_styles:autocad.ssi"
Fori=0ToUBound(Text)-1
HZ=Split(Text(n),",")
IfMid(HZ(0),1,2)=">>"Then
Print#1,0,0#,0#,0#,0#,0#,0#
Else
HZ=Split(Text(n),",")
Print#1,HZ(3)&","&HZ(4)&","&HZ(5)
EndIf
Nexti
MsgBox " changes complete ", " storage capacity section "
EndSub
Reservoir storage survey data can be realized so easily and be transformed into line file data layout .str file needed for three-dimensional model.
(2) the modeling of cave depot:
Solid model is by a series of adjacent triangular facet, is rolled into the entity that inside and outside is ventilative.As shown in Figure 7, entity is by a series of point on line, and be linked to be the triangulation network that inside and outside is ventilative, must there be adjacent triangular facet on the limit of any triangular facet, 3 summits of any triangular facet, must be attached on effective point, otherwise entity is open or invalid.Solid model is used to describe three-dimensional object, solid model is also the principle based on digital surface model (DTM), the a series of Points And lines formed by the point that line string comprises, is linked to be face between three points, and the profile of object is reacted in the set of all.Adopt which type of algorithm, make the set of triangular facet can react the profile of object, these triangular facets may be overlapping on plane angle, but think not overlapping or crossing in three dimensional network, triangle in solid model is a complete totally enclosed structure, it is defined as the model of solid or empty body, the shape of generation can be used for three-dimensional visible, calculate the cutting profile etc. of volume, arbitrary orientation
Table triangulation network constructive process complete list is as table 2:
Table 2triangularmeshcreationprocesslist
According to our measuring method, select profile line modeling, for cave depot, by section line modeling, create the triangulation network, form solid model.At the two ends of entity, be closed, just define the entity of cave depot, as shown in Figure 8.
(3) the validation verification of entity:
If each triangular facet of composition entity exists self intersection, without adjacent limit, repeats limit, invalid edges, then entity is not effective entity, and invalid entity can not calculate volume, space constraint, logical operation etc.
(4) cave depot physical manipulations
By the validation verification of entity, between entities, between entity and face, carry out the computing of the common factor in space, union and difference set.Spatial operation can be carried out between entity, between entity and face.As shown in Figure 9, introduce operation relation between them below: between entity and entity, if there is public part, the entity of public part can be drawn, one can be merged into, can from certain entity deduction another one.
Entity occurs simultaneously: A ∩ B
Entity union: A ∪ B
Entity difference set: A-B
Between entity and face:
Entity and face: A is on B
A is under B
Face and entity:
B is outside A
B is in A
By physical manipulations, first form the split figure of each cavern, then form the merging body figure of cave depot, finally form the solid model figure of cave depot.
Validation verification to be carried out to entity after having merged, be verified and .3dm preservation is carried out to model, so that later to model manipulation.
(5) the cutting of cave depot solid model
Definition, using Z axis as cross-sectional axis, carries out cutting with the vertical direction of axle, definition storage capacity measuring point elevation, input section absolute altitude, automatically will calculate and draw out the section of this absolute altitude.
The calculating of 2 cave depot aggregate storage capacities
2.1 modelings calculate automatically
According to the three-dimensional example illustraton of model of the cave depot built up, by the validation verification of entity, just can automatically generate entity report file, the volume complete list including cave depot cumulative volume in file and calculate by elevation spacing value.
Namely elevation spacing value: 0.01m(cuts into slices according to 1/cm section)
Body=: 1
The triangulation network=: 1
Verify=pass through
State=solid
Triangulation network coordinate range:
X is minimum: 219759.146X is maximum: 219903.713
Y is minimum: 101356.327Y is maximum: 101691.681
Z is minimum :-116.485Z is maximum :-87.259
Grid DEM volume complete list pressed by table 3
Table3byelevationtocalculatevolumelist
2.2 computation process
As shown in Figure 10, calculate all horizontal or vertical cross-sectional areas, then try to achieve aggregate storage capacity with the length integration of each measuring unit.
Computing formula is:
V 1 = ( S 0 + S 1 ) 2 × L 1 V = Σ i = 1 n ( Si × Li ) - - - ( 1 )
The aggregate storage capacity of artificial calculating and the aggregate storage capacity automatically calculated can comparatively validates mutually.
3. section survey positional accuracy and some bit density to the impact analysis of aggregate storage capacity
3.1 positional accuracy
During section survey, each data point that section is upsampled to is the most elementary cell forming storage capacity, and the set of point determines the aggregate storage capacity of cave depot, section is put the area measurement precision that the precision of position and Effects of Density section.Collect in Section measuring system each point coordinate by sampled point is carried out angle measurement, survey limit calculate, visible, the measuring accuracy of cross-sectional area is decided by the angle measurement of total powerstation, the precision of side, also has close relationship with the shape of section.
Underground water seal cave depot section is mainly with city gate type and the round two kinds of sections composition of three-core, and wherein type section in city gate is mainly distributed in and connects tunnel, has the section of 5 kinds of Different structural parameters altogether; Three-core is round is mainly distributed in main cave depot, and section structure parameter only has a kind, and its maximum height is 26 meters, breadth extreme 18 meters.Easy for analyzing, supposing that main cave depot is take radius as the circle cross section of R, then have: area S=π R2, if area measurement medial error be mS, radius measurement medial error is mR, according to law of propagation of errors, can obtain:
m s S = 2 m R R - - - ( 2 )
If cross-sectional area measuring accuracy is 1,/10,000, then the cave depot of 20m radius, it surveys edge error can not be greater than ± 1mm; Equally, if the direction mean square error of angle observation of instrument is m β, then have:
m s S = m β πρ 2 - - - ( 3 )
Known, the direction mean square error of angle observation of instrument is less, also less on the impact of cross-sectional area measuring accuracy.Therefore suggestion adopts high-precision I level intelligence automatic total instruments to measure.
3.2 bit densitys
Point bit density be exactly sampling spacing distance between points, in Section measuring system namely total powerstation to the angle step value of sampled point.After cave depot construction, do not carry out secondary lining, hole wall has only carried out first spray, and hole wall surface is coarse uneven, and sampled point is exactly measure on coarse surface, and concavo-convex difference is not of uniform size.There is representive error problem for section sampling like this, concerning whole cave depot, concavo-convex difference meets normal distribution law, just can be solved the representive error of cave depot appearance like this by the density of constraint fixed sample point position.The maximum radius of underground cavern section is 26 meters, minimum be 4 meters, actual when gathering section point by radius size set point between arc length be 0.2--0.5m, corresponding angle step value is 30--10.Under the condition of same measuring accuracy and sample density, cross sectional radii is larger, and the angle step value of sampled point can suitably amplify.The density of point increases, and curve between points just can shorten as far as possible, so just more can embody actual profile situation.
4. the density of Measure section and precision are to the impact analysis of aggregate storage capacity
Hold manually to calculate cave depot, volume folded between its two parallel cutting surfaces in conventional section is pressed trapezoid cylinder method and is calculated.Differentiated by above-mentioned computing formula (1) and medial error mV can be obtained:
mV = | ΔL 2 ( S 0 2 + S 1 2 ) 4 · ( m s S ) 2 + ( S 0 + S 1 ) 2 4 · ΔL 2 · ( m ΔL ΔL ) 2 | - - - ( 4 )
The phase error of centralization that be can be derived from cave depot volume by (4) formula is:
m V V = | | ( m ΔL L ) 2 + S 0 2 + S 1 2 ( S 0 + S 1 ) 2 ( m S S ) 2 | | - - - ( 5 )
As can be seen from formula, the phase error of centralization of cave depot volume is mainly by the phase error of centralization of cave depot length L direction increment and the impact of the cross-sectional area phase error of centralization.So, for ensureing the relative accuracy of cave depot aggregate storage capacity, certain Measure section density will be considered.
Its beneficial effect of comprehensive technique scheme is: at present, and the method solving underground water seal cave depot reservoir storage survey problem, except above-mentioned section survey method, also has one to be exactly utilize three-dimensional laser scanner to carry out three-dimensional measurement, then resolves through set.But this instrument belongs to new things, expensive, being widely used also has some limitations.Underground water seal cave depot belongs to construction frontier, and section survey method is relatively convenient, easy to operate, and by analyzing, the phase error of centralization of cross-sectional area also easily ensures.
In the description of this instructions, specific features, structure, material or feature that the description of reference term " embodiment ", " some embodiments ", " example ", " concrete example " or " some examples " etc. means to describe in conjunction with this embodiment or example are contained at least one embodiment of the present invention or example.In this manual, identical embodiment or example are not necessarily referred to the schematic representation of above-mentioned term.And the specific features of description, structure, material or feature can combine in an appropriate manner in any one or more embodiment or example.
Although illustrate and describe embodiments of the invention, those having ordinary skill in the art will appreciate that: can carry out multiple change, amendment, replacement and modification to these embodiments when not departing from principle of the present invention and aim, scope of the present invention is by claim and equivalents thereof.

Claims (5)

1. a reservoir storage survey method for underground water seal cave depot, is characterized in that, comprise the steps:
Step 1, sets up underground water seal cave depot SURVEYING CONTROL NETWORK;
Step 2, after having set up SURVEYING CONTROL NETWORK, definition underground water seal cave depot section survey parameter, scans underground water seal cave depot;
Data after scanning are arranged, finally obtain the three-dimensional data of underground water seal cave depot by step 3;
Step 4, carries out three-dimensional modeling to the three-dimensional data of underground water seal cave depot shaping, namely sets up some adjacent triangulation networks, forms solid model, then carries out according to the axle vertical direction of definition the section that cutting forms corresponding height;
Step 4-1, adopt measured database to carry out manipulate measurement data, cross-sectioning data are converted to line file form, each the section point first section of being combined in line file, all sections by a cave depot form line string;
Step 4-2, by some adjacent triangular facets, is rolled into the entity that inside and outside is ventilative; Entity is by a series of point on line, and be linked to be the triangulation network that inside and outside is ventilative, must there be adjacent triangular facet on the limit of any triangular facet, 3 summits of any triangular facet, must be attached on effective point;
Described step 4-2 also comprises the step creating the triangulation network:
Step 4-3, between two sections: connect the triangulation network between different sections;
Step 4-4, in a section: the connection automatically completing the triangulation network in a section chosen;
Step 4-5, section is to a point: a section is to a some connection triangulation network;
Step 4-6, between multiple sections: automatically connect the triangulation network between a series of line and section;
Step 4-7, uses control line: define section with control line, and then produces the triangulation network;
Step 4-8, according to manual reconnaissance: the manual definition triangulation network retrains, and connects the triangulation network between section;
Step 4-9, single triangle: define each triangle by selected triangular apex;
Step 4-10, from a section to two sections: a closed female section and two sub-sections or point, use associating concept, control branch line, and then the generation triangulation network;
Step 4-11, from a section to multiple sections: connect the triangulation network a female section and multiple sub-section;
Step 4-12, uses center line and section: utilize section to connect the triangulation network along center line, generally generates underworkings entity by this function;
Step 4-13, if each triangular facet of composition entity exists self intersection, without adjacent limit, repeats limit, invalid edges, then entity is not effective entity, and invalid entity can not calculate volume, space constraint, logical operation;
Step 4-14, by the validation verification of entity, between entities, the computing of the common factor in space, union and difference set is carried out between entity and face, definition is using Z axis as cross-sectional axis, cutting is carried out with the vertical direction of axle, definition storage capacity measuring point elevation, input section absolute altitude, automatically will calculate and draw out the section of this absolute altitude;
Step 5, calculate the area of section, with the height value of corresponding equalization, obtain the volume of described section, according to the overall height of the underground water seal cave depot measured, obtain underground water seal cave depot aggregate storage capacity, the positional accuracy of adjustment underground water seal cave depot section survey and some bit density obtain the accurate storage capacity volume of underground water seal cave depot.
2. the reservoir storage survey method of underground water seal cave depot according to claim 1, is characterized in that, described step 1 comprises:
Step 1-1, ground control point or in-cavity construction reference mark are that benchmark sets up unified SURVEYING CONTROL NETWORK to whole cave depot, to the encryption needing the position measured to carry out reference mark on the basis of Controling network, adopt accurate fourth class tracerse survey method, after compensating computation, obtain the three-dimensional coordinate at each reference mark.
3. the reservoir storage survey method of underground water seal cave depot according to claim 1, is characterized in that, described step 2 comprises:
Step 2-1, after having set up underground water seal cave depot SURVEYING CONTROL NETWORK, according to the design transversal section parameter of design data definition vertical curve, horizontal curve and each measuring point, obtains the design reference section closed;
Step 2-2, definition vertical curve, wherein the data of vertical curve are mileage, elevation, the radius of knick point and tangent line, Origin And Destination;
Step 2-3, definition horizontal curve, wherein horizontal curve is the data for describing, determining center line of road accurate location, and horizontal curve data are mileage, line style, turning radius, the coordinate of line element starting point;
Step 2-4, definition transversal section parameter, wherein transversal section supplemental characteristic is that right-handed coordinate system is tried to achieve, coordinate system with circuit at section mileage place designed elevation point for initial point, with normal on the right side of line midline working direction for X-axis, with initial point zenith direction for Y-axis, during In-put design section, according to inputting every section of straight line and circular arc successively clockwise.
4. the reservoir storage survey method of underground water seal cave depot according to claim 1, is characterized in that, described step 3 comprises:
Step 3-1, total powerstation is set up near center line and centering leveling bottom underground water seal cave depot, adopt the method for Free Station, obtain the three-dimensional coordinate of survey station point, select horizontal curve, vertical curve, transversal section type, arrange according to the dot spacing of 1 meter, mileage direction with vertical direction 0.5 meter, then carry out measurement scanning;
Step 3-2, record transversal section mileage, the period of measurement point, partially middle distance, the discrepancy in elevation, three-dimensional coordinate.
5. the reservoir storage survey method of underground water seal cave depot according to claim 1, is characterized in that, described step 5 comprises:
Step 5-1, calculates all horizontal or vertical cross-sectional areas,
wherein V 1for initial computing volume, S 0be the first cross sectional area, S 1be the area of the second transversal section, L 1be the first computing length,
Aggregate storage capacity is tried to achieve again with the length integration of each measuring unit,
v is aggregate storage capacity volume, and i, n are positive integer;
Step 5-2, after calculating underground water seal cave depot aggregate storage capacity, adjustment positional accuracy and some bit density;
Step 5-3, carries out the positional accuracy of section survey and puts bit density to the impact analysis of aggregate storage capacity; During section survey, each data point that section is upsampled to is the most elementary cell forming storage capacity, and the set of point determines underground water seal cave depot aggregate storage capacity, section is put the area measurement precision that the precision of position and Effects of Density section; Collect each point coordinate by sampled point is carried out angle measurement, survey limit calculate;
Step 5-4, underground water seal cave depot section is with city gate type and three-core round section composition, and wherein type section in city gate is mainly distributed in and connects tunnel, if main cave depot take radius as the circle cross section of R, then has:
Area S=π R 2,
If area measurement medial error is m s, radius measurement medial error is m r, according to law of propagation of errors, can obtain:
wherein subscript S is for measuring area, subscript R for measuring radius;
If cross-sectional area measuring accuracy is 1/10000, then the cave depot of 20m radius, it surveys edge error can not be greater than ± 1mm; Equally, if the direction mean square error of angle observation of instrument is m β, then have:
wherein subscript S measures the direction angle measurement that area, subscript β are instrument;
When direction, mean square error of angle observation is less, also less on the impact of cross-sectional area measuring accuracy;
Step 5-5, some bit density is the spacing distance between points of sampling, to the angle step value of sampled point in section survey; Because hole wall surface is coarse uneven, sampled point is exactly measure on coarse surface, and concavo-convex difference is not of uniform size; There is representive error problem for section sampling, concerning whole cave depot, concavo-convex difference meets normal distribution law, just can be solved the representive error of cave depot appearance like this by the density of constraint fixed sample point position; The maximum radius of underground cavern section is 26 meters, minimum be 4 meters, actual when gathering section point by radius size set point between arc length be 0.2--0.5m, corresponding angle step value is 30-10; Under the condition of same measuring accuracy and sample density, cross sectional radii is larger, and the angle step value of sampled point can suitably amplify; The density of point increases, and curve between points just can shorten as far as possible;
Step 5-6, hold manually to calculate cave depot, volume folded between its two parallel cutting surfaces in conventional section is pressed trapezoid cylinder method and is calculated, by differentiate and can obtain medial error m v, wherein subscript V is for measuring aggregate storage capacity volume:
wherein subscript S is for measuring area,
Above formula pushes away the phase error of centralization of cave depot volume is:
wherein subscript S is for measuring area,
The phase error of centralization of cave depot volume, mainly by the phase error of centralization of cave depot length L direction increment and the impact of the cross-sectional area phase error of centralization, for ensureing the positional accuracy of cave depot aggregate storage capacity, needs Measure section point bit density.
CN201310271053.2A 2013-06-28 2013-06-28 The reservoir storage survey method of underground water seal cave depot Active CN103344297B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310271053.2A CN103344297B (en) 2013-06-28 2013-06-28 The reservoir storage survey method of underground water seal cave depot

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310271053.2A CN103344297B (en) 2013-06-28 2013-06-28 The reservoir storage survey method of underground water seal cave depot

Publications (2)

Publication Number Publication Date
CN103344297A CN103344297A (en) 2013-10-09
CN103344297B true CN103344297B (en) 2015-12-02

Family

ID=49279113

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310271053.2A Active CN103344297B (en) 2013-06-28 2013-06-28 The reservoir storage survey method of underground water seal cave depot

Country Status (1)

Country Link
CN (1) CN103344297B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103643935B (en) * 2013-12-20 2016-04-06 山东大学 Underground water seal oil storage cave depot identification of formation and water curtain hole optimization device and method
CN104268344B (en) * 2014-09-28 2018-04-17 中国水利水电科学研究院 A kind of Groundwater Reservoir storage capacity computational methods
CN104989428A (en) * 2015-07-20 2015-10-21 顾达容 Tunnel construction lining quantity determining device and tunnel construction lining quantity determining construction method
CN107525557B (en) * 2017-10-19 2018-08-17 中国矿业大学(北京) A kind of accurate method for measuring coal mine underground reservoir capacity
CN108168429B (en) * 2017-12-20 2020-06-09 中国计量大学 Large cave depot capacity measuring method based on three-dimensional laser scanner
CN113481799B (en) * 2021-07-21 2022-06-10 中国路桥工程有限责任公司 Concrete road surface maintenance device

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009192332A (en) * 2008-02-13 2009-08-27 Konica Minolta Sensing Inc Three-dimensional processor and method for controlling display of three-dimensional data in the three-dimensional processor
CN101882171A (en) * 2010-05-21 2010-11-10 中交第二公路勘察设计研究院有限公司 Method for fast establishing interactive tunnel and wall rock body three-dimensional models
CN102174893A (en) * 2011-03-16 2011-09-07 中铁第一勘察设计院集团有限公司 Testing method for laying ballastless tracks in early stage before overall completion of super long tunnel
CN102622526A (en) * 2012-04-06 2012-08-01 中国矿业大学(北京) Digital mine tunneling search prediction method
CN102831647A (en) * 2012-08-06 2012-12-19 电子科技大学 Delaunay triangulation network subdividing method based on space curved surface constraint

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4219347B2 (en) * 2005-08-17 2009-02-04 富士通株式会社 Measurement processing method and measurement apparatus

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009192332A (en) * 2008-02-13 2009-08-27 Konica Minolta Sensing Inc Three-dimensional processor and method for controlling display of three-dimensional data in the three-dimensional processor
CN101882171A (en) * 2010-05-21 2010-11-10 中交第二公路勘察设计研究院有限公司 Method for fast establishing interactive tunnel and wall rock body three-dimensional models
CN102174893A (en) * 2011-03-16 2011-09-07 中铁第一勘察设计院集团有限公司 Testing method for laying ballastless tracks in early stage before overall completion of super long tunnel
CN102622526A (en) * 2012-04-06 2012-08-01 中国矿业大学(北京) Digital mine tunneling search prediction method
CN102831647A (en) * 2012-08-06 2012-12-19 电子科技大学 Delaunay triangulation network subdividing method based on space curved surface constraint

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
三维建模技术在容积计量中的应用研究;余成道;《中国优秀博硕士学位论文全文数据库(硕士) 基础科学辑》;20060515(第5期);第12-48页 *
洞库容量测量若干问题的研究;卫建东等;《解放军测绘学院学报》;19991231;第16卷(第4期);第252-256页 *

Also Published As

Publication number Publication date
CN103344297A (en) 2013-10-09

Similar Documents

Publication Publication Date Title
CN103344297B (en) The reservoir storage survey method of underground water seal cave depot
CN102251451B (en) Track geometric state measurement system and method based on multisource information fusion technology
CN114329749B (en) Method for horizontal-vertical integrated auxiliary optimization design of existing railway
CN102661736B (en) Highway reorganization and expansion surveying method
CN103217688B (en) Airborne laser radar point cloud adjustment computing method based on triangular irregular network
CN106767686B (en) Structural settlement monitoring method and system based on inclination angle and structural joint deformation measurement
CN105043342A (en) Unidirectional precise distance measuring triangulated height measuring method
CN111968233B (en) BIM and point cloud positioning and track interval communication equipment installation method
CN107815935A (en) High speed railway track geometry status real time monitor method and system
CN107700280B (en) Existing double railway lines line position reconstructing method
CN106909740B (en) Method and device for calculating total map earthwork amount of oil and gas storage and transportation station
CN103727920B (en) The method of height differences by leveling is measured based on geoid's model
CN107330140A (en) The method that transformer station is quickly vertically arranged is realized based on BIM technology
CN112461205A (en) Method for manufacturing cross section of existing railway line based on unmanned aerial vehicle oblique photogrammetry
CN112100715A (en) Three-dimensional oblique photography technology-based earthwork optimization method and system
CN107945264A (en) Roadbed three-dimensional modeling method
CN110598239B (en) Application method based on track area point cloud big data
CN102620721A (en) Fine digital terrain model based road surveying method
CN107144265B (en) A kind of measuring method based on the graphic three-dimensional property right body space boundary mark of projection
CN102750413A (en) Data processing and mapping method of topographic surveying of electric transmission line tower positions
CN103196425A (en) Estimation method of extra-long tunnel horizontal through error
CN101539410A (en) Existing railway survey method based on total station
CN111854692A (en) Method for measuring unmanned aerial vehicle image matching point cloud in road survey
CN111426303A (en) Karst slope parameter measuring method
CN108917706A (en) A kind of GNSS leveling measuring method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant