CN109584365B - 一种用于海洋吹填造陆工程土方量造价核算的Voronoi图地质分层模型构建方法 - Google Patents
一种用于海洋吹填造陆工程土方量造价核算的Voronoi图地质分层模型构建方法 Download PDFInfo
- Publication number
- CN109584365B CN109584365B CN201811440239.5A CN201811440239A CN109584365B CN 109584365 B CN109584365 B CN 109584365B CN 201811440239 A CN201811440239 A CN 201811440239A CN 109584365 B CN109584365 B CN 109584365B
- Authority
- CN
- China
- Prior art keywords
- construction
- geological
- layer
- voronoi diagram
- volume
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种用于海洋吹填造陆工程土方量造价核算的Voronoi图地质分层模型构建方法,其针对海洋吹填造陆工程,通过构建钻孔采样点所在局部分块的Voronoi图地质分层模型,并联合双层三角网模型,实现了海洋吹填造陆工程中不同地质层的土方量造价准确、快速计算。本发明步骤是首先构建施工前与施工后海底表面地形的双层三角网模型,然后基于Voronoi图的最邻近特性构建中间地质层的Voronoi图分层模型,得到各地质层的土方量,最后根据不同地质层的单位土方计价得到海洋吹填造陆工程土方量的总造价。
Description
技术领域
本发明属于地理信息系统(GIS)中空间分析研究领域,涉及到一种用于海洋吹填造陆工程土方量造价核算的Voronoi图地质分层模型构建方法。
背景技术
土方量造价核算是工程建设方和施工方结算工程量、避免经济纠纷的主要依据,直接关系到建设方和施工方的经济利益。当前,土方量造价核算方法都是根据施工前与施工后区域的表面地形数据,利用等高线、方格网、断面或三角网等测算方法计算出该区域的土方量,并根据单位土方的计价得出该施工区域的土方量造价。然而,现有土方量造价核算方法应用于海洋吹填造陆工程极易产生经济纠纷,这是因为海底沉积物从上至下容易形成砂砾、碎块、岩石等颗粒大小不同的地质层,不同地质层的开挖方法和使用的船只工具不尽相同,导致不同地质层的施工难度和计价成本均不同。
为了准确核算出海洋吹填造陆工程土方量造价,除了利用多波束测深系统获取施工前与施工后区域的表面地形数据外,还需利用打桩钻孔测量技术获取施工前采样点处不同地质层的深度值。由于每打桩钻孔测量一个采样点所花费的时间、财力、物力等成本非常昂贵,所以采样的钻孔点数量就相当稀少。利用小样本的钻孔点全局插值计算施工后区域表面地形的大样本采样点处的不同地质层深度值,再利用现有土方量测算方法虽然能够核算出海洋吹填造陆工程土方量造价,但存在分层计算过程复杂、计算量大等问题,同时,利用小样本采样点全局插值计算大样本采样点已不具有实际应用意义。因此,本发明提出了一种用于海洋吹填造陆工程土方量造价核算的Voronoi图地质分层模型构建方法,通过构建钻孔采样点所在局部分块的Voronoi图地质分层模型,并基于Voronoi图的最邻近特性,即位于Voronoi图中同一多边形分块内的点具有相同的特性,利用钻孔点的地质分层属性值表示其所在局部分块的地质分层属性值,实现了海洋吹填造陆工程中不同地质层的土方量造价准确、快速计算。
发明内容
本发明要解决的技术问题是:克服现有技术的上述不足,提出了一种用于海洋吹填造陆工程土方量造价核算的Voronoi图地质分层模型构建方法,其能够以较快的速度准确计算不同地质层的土方量造价,为海洋吹填造陆工程建设方和施工方结算工程量提供了可靠的依据。
为实现上述发明目的,本发明提供了一种用于海洋吹填造陆工程土方量造价核算的Voronoi图地质分层模型构建方法,包括以下步骤:
步骤1,数据输入,读取施工前的海底表面水深点集数据、施工后的海底表面水深点集数据、钻孔点集数据和施工边界数据,所述点数据具有各自的序号、坐标数据和水深值数据,所述施工边界数据为简单多边形数据,由多个点数据经首尾顺次相连而成;
步骤2,空间内插,利用施工前的海底表面水深点集数据、施工后的海底表面水深点集数据分别插值计算施工前与施工后的施工边界上点坐标的水深值;
步骤3,剔除判定,分别将施工前的海底表面水深点集数据、施工后的海底表面水深点集数据中位于施工边界范围外的点数据剔除;
步骤4,构建双层三角网模型,方法为:先利用施工边界及其内的施工前的海底表面水深点集数据,构建第零层三角网模型;再利用施工边界及其内的施工后的海底表面水深点集数据,构建最底层三角网模型,所述双层三角网模型构建完毕后,采用三角网法分别计算第零层三角网、最底层三角网与基准海平面构成的体积;
步骤5,构建基于泰森多边形Voronoi图的中间地质层分层模型,方法为:利用施工边界及其内的钻孔点集数据,构建中间地质层的Voronoi图分层模型,得到中间各地质层与基准海平面构成的体积。该方法主要包括构建钻孔点的Voronoi图、构建地质分层模型两个处理环节;
步骤6,土方量造价计算,在计算得到第零层、中间各地质层、最底层与基准海平面构成的体积后,利用两期体积相减法,即下一层体积减去上一层体积,得到不同地质层的真实土方量,并根据不同地质层的单位土方计价,计算各地质层的土方量造价,累加得到海洋吹填造陆工程土方量的总造价。
进一步地所述步骤5中,构建基于泰森多边形Voronoi图的中间地质层分层模型的方法如下:
步骤a,利用钻孔点集数据,构建钻孔点的Voronoi图,并与施工边界进行空间叠置裁剪,得到施工边界范围内钻孔点的Voronoi图;
步骤b,基于Voronoi图的最邻近特性,利用钻孔点的地质分层属性值表示其所在局部分块的地质分层属性值,从而建立Voronoi图中多边形的地质分层模型,得到Voronoi图中每个多边形内的不同地质层的水深属性值。根据多边形的面积和某一地质层的水深属性值可得到该层的体积,并按照分层次序进行体积累加,得到各地质层与基准海平面构成的体积。
本发明与现有技术相比的优点在于:现有土方量造价核算方法通常都忽视了地质分层的影响,即使在海洋吹填造陆工程中考虑了地质分层的影响,但利用小样本的钻孔点全局插值计算施工后区域表面地形的大样本采样点,再利用现有土方量测算方法核算出整体施工区域内不同地质层的土方量造价,存在分层计算过程复杂、计算量大等问题,同时,利用小样本采样点全局插值计算大样本采样点已不具有实际应用意义。本发明构建钻孔采样点的Voronoi图地质分层模型,采用局部分块、分层模型替代了全局整体模型,并联合双层三角网模型,能够准确计算不同地质层的土方量造价,减小了计算量,为海洋吹填造陆工程建设方和施工方结算工程量提供了可靠的计算结果。
附图说明
图1为本发明的方框示意图;
图2为本发明中钻孔采样点的地质层结构及其相应水深值;
图3为本发明中施工边界及其上的坐标点;
图4为本发明基于施工边界的约束三角网;
图5为本发明中基于钻孔点集的Voronoi图区域划分;
具体实施方式
为了更好地实现上述发明目的,下面将结合图1的方框图对本发明中的各个步骤做进一步描述:
步骤1,数据输入,读取施工前的海底表面水深点集数据、施工后的海底表面水深点集数据、钻孔点集数据和施工边界数据,所述数据具有各自的存储结构:施工前与施工后的海底表面水深点集数据都具有相同的存储结构,每个点数据都由序号、x与y坐标、单一水深值组成;钻孔点数据由序号、x与y坐标、多个水深值组成。图2显示的是某钻孔点处的三层地质层结构,其存储结构表示为(id,x,y,h0,h1,h2,h3),对应于海底沉积物从海底表面向下更深地质层依次是砂砾层(第1层)、碎块层(第2层)和岩石层(第3层),第0层为基准海平面到海底表面这一立体空间,h0表示钻孔点处基准海平面到海底表面的水深值,h1、h2、h3分别表示砂砾层(第1层)、碎块层(第2层)和岩石层(第3层)的水深值。施工边界数据为简单多边形数据,由多个点数据经首尾顺次相连而成,每个点的存储结构与海底表面水深点存储结构相同。图3显示的是施工边界多边形及其上的坐标点。
步骤2,空间内插,利用施工前的海底表面水深点集数据,通过反距离权重插值方法计算施工边界上点数据的水深值。同理,利用施工后的海底表面水深点集数据,通过反距离权重插值方法计算施工边界上点数据的水深值。
步骤3,剔除判定,叠加施工前的海底表面水深点集数据和施工边界,通过点-面包含判定方法剔除位于施工边界范围外的水深点数据,得到施工边界范围内的施工前的海底表面水深点集数据。同理,剔除施工后的海底表面水深点集数据中位于施工边界范围外的点数据,得到施工边界范围内的施工后的海底表面水深点集数据。
步骤4,构建双层三角网模型,首先利用施工边界及其内的施工前的海底表面水深点集数据,构建施工边界约束范围内的三角网模型,即第零层三角网模型,如图4所示;再采用公式(1)计算每个三角形与基准海平面构成的三棱柱体积;最后累加三棱柱体积得到第零层三角网与基准海平面构成的体积。同理,利用施工边界及其内的施工后的海底表面水深点集数据,构建最底层三角网模型,得到最底层三角网与基准海平面构成的体积。
上述公式(1)中,vi表示第i个三角形Δabc与基准海平面构成的三棱柱体积,si表示第i个三角形Δabc的面积,za、zb、zc分别表示三角形Δabc的三个顶点a、b、c处的水深值。
步骤5,构建Voronoi图地质分层模型,首先利用钻孔点集数据,构建钻孔点的Voronoi图,并与施工边界进行空间叠置裁剪,得到施工边界范围内钻孔点的Voronoi图,如图5所示;然后,基于Voronoi图的最邻近特性,即位于Voronoi图中同一多边形分块内的点具有相同的特性,利用钻孔点的地质分层的水深属性值表示其所在多边形分块的地质分层属性值,从而建立Voronoi图中多边形的地质分层模型,得到Voronoi图中每个多边形内的不同地质层的水深属性值,并计算多边形内的不同地质层与基准海平面构成的体积;最后采用公式(2)累加各层次体积,得到中间各地质层与基准海平面构成的体积。
上述公式(2)中,Vk表示施工边界范围内的第k地质层与基准海平面构成的体积,Si表示Voronoi图中第i个多边形的面积,hk表示钻孔点的第k地质层的水深值,M表示钻孔点或Voronoi图中多边形的数量,K表示地质层数量。
步骤6,土方量造价计算,利用两期体积相减法,如公式(3)所示,得到不同地质层的土方量,并根据不同地质层的单位土方计价,计算各地质层的土方量造价,累加得到海洋吹填造陆工程土方量的总造价。
V'k=Vk-Vk-1(1≤k≤K) (3)
上述公式(3)中,V'k表示第k地质层的土方量。
本发明未详细阐述部分属于本领域技术人员的公知技术。
以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可理解想到的变换或替换,都应涵盖在本发明的包含范围之内,因此,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (1)
1.一种用于海洋吹填造陆工程土方量造价核算的Voronoi图地质分层模型构建方法,其特征在于:先构建双层三角网模型,精确计算海洋吹填造陆工程总的土方量,再构建Voronoi图地质分层模型,准确计算各地质层的土方量,最后根据不同地质层的单位土方计价得到海洋吹填造陆工程的土方量总造价,具体步骤如下:
步骤1,数据输入,读取施工前的海底表面水深点集数据、施工后的海底表面水深点集数据、钻孔点集数据和施工边界数据,所述点数据具有各自的序号、坐标数据和水深值数据,所述施工边界数据为简单多边形数据,由多个点数据经首尾顺次相连而成;
步骤2,空间内插,利用施工前的海底表面水深点集数据、施工后的海底表面水深点集数据,分别插值计算施工前与施工后的施工边界上点坐标的水深值;
步骤3,剔除判定,分别将施工前的海底表面水深点集数据、施工后的海底表面水深点集数据中位于施工边界范围外的点数据剔除;
步骤4,构建双层三角网模型,方法为:先利用施工边界及其内的施工前的海底表面水深点集数据,构建第零层三角网模型;再利用施工边界及其内的施工后的海底表面水深点集数据,构建最底层三角网模型,所述双层三角网模型构建完毕后,采用三角网法分别计算第零层三角网、最底层三角网与基准海平面构成的体积;
步骤5,构建基于泰森多边形Voronoi图的中间地质层分层模型,包括构建钻孔点的Voronoi图、构建地质分层模型两个处理环节,具体实现为:利用施工边界及其内的钻孔点集数据,构建中间地质层的Voronoi图分层模型,得到中间各地质层与基准海平面构成的体积;
步骤6,土方量造价计算,在计算得到第零层、中间各地质层、最底层与基准海平面构成的体积后,利用两期体积相减法,即下一层体积减去上一层体积,得到不同地质层的真实土方量,并根据不同地质层的单位土方计价,计算各地质层的土方量造价,累加得到海洋吹填造陆工程土方量的总造价;
所述步骤5中,构建基于泰森多边形Voronoi图的中间地质层分层模型的方法如下:
步骤a,利用钻孔点集数据,构建钻孔点的Voronoi图,并与施工边界进行空间叠置裁剪,得到施工边界范围内钻孔点的Voronoi图;
步骤b,基于Voronoi图的最邻近特性,利用钻孔点的地质分层属性值表示其所在局部分块的地质分层属性值,从而建立Voronoi图中多边形的地质分层模型,得到Voronoi图中每个多边形内的不同地质层的水深属性值;根据Voronoi图中多边形的面积和某一地质层的水深属性值得到该层的体积,并按照分层次序进行体积累加,得到各地质层与基准海平面构成的体积。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811440239.5A CN109584365B (zh) | 2018-11-29 | 2018-11-29 | 一种用于海洋吹填造陆工程土方量造价核算的Voronoi图地质分层模型构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811440239.5A CN109584365B (zh) | 2018-11-29 | 2018-11-29 | 一种用于海洋吹填造陆工程土方量造价核算的Voronoi图地质分层模型构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109584365A CN109584365A (zh) | 2019-04-05 |
CN109584365B true CN109584365B (zh) | 2022-05-10 |
Family
ID=65925060
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811440239.5A Active CN109584365B (zh) | 2018-11-29 | 2018-11-29 | 一种用于海洋吹填造陆工程土方量造价核算的Voronoi图地质分层模型构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109584365B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111311746B (zh) * | 2019-12-09 | 2023-09-08 | 中交广州航道局有限公司 | 一种基于钻孔数据的智能三维地质建模方法 |
CN113235655B (zh) * | 2021-05-10 | 2023-03-28 | 中交广州航道局有限公司 | 海底沉管隧道基槽分土质工程土方量计算方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103093447A (zh) * | 2013-01-18 | 2013-05-08 | 南京大学 | 一种计算机图片集中图片剪切拼接方法 |
CN103927783A (zh) * | 2014-04-17 | 2014-07-16 | 武汉大学 | 一种三维三角网构建填挖空间的图割方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7307629B2 (en) * | 2003-03-31 | 2007-12-11 | Council Of Scientific And Industrial Research | Generation of three dimensional fractal subsurface structure by Voronoi Tessellation and computation of gravity response of such fractal structure |
-
2018
- 2018-11-29 CN CN201811440239.5A patent/CN109584365B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103093447A (zh) * | 2013-01-18 | 2013-05-08 | 南京大学 | 一种计算机图片集中图片剪切拼接方法 |
CN103927783A (zh) * | 2014-04-17 | 2014-07-16 | 武汉大学 | 一种三维三角网构建填挖空间的图割方法 |
Non-Patent Citations (4)
Title |
---|
"基于ArcGIS Engine的地上下集成建模关键技术研究";张军强;《中国优秀硕士学位论文全文数据库》;20130115;全文 * |
Application of Voronoi and Delaunay Diagrams in Gravimetric Geoid Determination;Newton Pereira dos Santos ,Iris Pereira Escobar;《 3rd International Symposium on Voronoi Diagrams in Science and Engineering (ISVD 2006)》;20060101;全文 * |
Delaunay-三角网构建及其在土地平整中的应用;张渭军;《矿业研究与开发》;20110228(第01期);全文 * |
基于B/S的工程地质三维建模系统与应用;黄牧等;《工程勘察》;20171101(第11期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109584365A (zh) | 2019-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Virbulis et al. | Hydrogeological model of the baltic artesian basin | |
CA2823989C (en) | Method and system of updating a geocellular model | |
CN108647473A (zh) | 一种盾构法修建海底隧道合理埋深计算方法 | |
CN108842821B (zh) | 一种钻爆法修建海底隧道合理埋深的计算方法 | |
CN108986213B (zh) | 一种基于叠置技术的三维地层建模方法 | |
CN109584365B (zh) | 一种用于海洋吹填造陆工程土方量造价核算的Voronoi图地质分层模型构建方法 | |
CN112465972B (zh) | 基于bim模型及桩基建模进行地质研判方法 | |
Sundell et al. | A probabilistic approach to soil layer and bedrock-level modeling for risk assessment of groundwater drawdown induced land subsidence | |
CN104634295A (zh) | 碳酸盐岩洞穴型储层有效体积估算方法 | |
CN114332391A (zh) | 一种三维地质体概率模型建模方法 | |
CN104899454A (zh) | 一种适用于土岩组合地层的地铁隧道裂隙水渗水量预测方法 | |
CN113433589A (zh) | 一种基于数理统计的风化壳储层底界面识别方法 | |
CN115019004A (zh) | 一种基于储层地质建模的油气储量计算方法 | |
CN110687591B (zh) | 基于先验数据的波形匹配确定煤层及围岩物性参数的方法 | |
Serpescu et al. | 3D geological model of Bucharest city Quaternary deposits | |
Zhang et al. | A novel three-dimensional mine area hydrogeological model based on groundwater modeling systems | |
Cowdery et al. | Hydrologic change in the St. Louis River Basin from iron mining on the Mesabi Iron Range, northeastern Minnesota | |
Keay et al. | 14 Challenges of Port Landscapes. Integrating Geophysics, Open Area Excavation and Computer Graphic Visualisation at Portus and the Isola Sacra | |
Matildes et al. | The contribution of GeoSIS_Lx database and geoscientific information system to the reconstruction of Lisbon’s coastline | |
Stafleu et al. | AN INTEGRATED MODELLING APPROACH AT TNO–GEOLOGICAL SURVEY OF THE NETHERLANDS | |
Bersezio et al. | Aquifer analogues to assist modeling of groundwater flow: the Pleistocene aquifer complex of the Agri Valley (Basilicata) | |
CN117930327A (zh) | 一种激发井深提取方法、装置和电子设备 | |
Chapple | 3D GEOLOGICAL MAPPING OF FRACTURE NETWORKS IN THE HAWKESBURY SANDSTONE AND THEIR IMPLICATIONS TO GROUNDWATER FLOW | |
Hou et al. | Refined Integration of Multilayer DEM-GTP-TEN for Three-Dimensional Quaternary Geological Modeling: A Case Study in Nantong City, China | |
Liu et al. | 3D Geological Modelling and Management System |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |