CN110263365A - 一种基于离散光滑插值技术的变形监测数据处理及分析方法 - Google Patents

一种基于离散光滑插值技术的变形监测数据处理及分析方法 Download PDF

Info

Publication number
CN110263365A
CN110263365A CN201910371971.XA CN201910371971A CN110263365A CN 110263365 A CN110263365 A CN 110263365A CN 201910371971 A CN201910371971 A CN 201910371971A CN 110263365 A CN110263365 A CN 110263365A
Authority
CN
China
Prior art keywords
deformation
analysis
monitoring
data processing
model
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.)
Pending
Application number
CN201910371971.XA
Other languages
English (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.)
PowerChina Huadong Engineering Corp Ltd
Original Assignee
PowerChina Huadong Engineering Corp 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 PowerChina Huadong Engineering Corp Ltd filed Critical PowerChina Huadong Engineering Corp Ltd
Priority to CN201910371971.XA priority Critical patent/CN110263365A/zh
Publication of CN110263365A publication Critical patent/CN110263365A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Architecture (AREA)
  • Mathematical Analysis (AREA)
  • Structural Engineering (AREA)
  • Pure & Applied Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Civil Engineering (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

本发明涉及一种基于离散光滑插值技术的变形监测数据处理及分析方法,适用于交通、水利、水电、矿山等行业复杂地质条件下地下洞室和边坡的变形监测数据处理和分析。本发明提供一种基于离散光滑插值技术的变形监测数据处理及分析方法,通过有限点的变形监测数据插值获得典型断面和三维全局的变形场分析特征,并通过比对数值模型的反演分析而获得的变形场分布特征,从而实现从局部和全局的维度来分析评价监测数据本身的合理性以及地下洞室围岩、边坡的稳定特征。

Description

一种基于离散光滑插值技术的变形监测数据处理及分析方法
技术领域
本发明涉及一种基于离散光滑插值技术的变形监测数据处理及分析方法,适用于交通、水利、水电、矿山等行业复杂地质条件下地下洞室和边坡的变形监测数据处理和分析。
背景技术
岩土工程中,地下洞室和边坡变形监测数据的分析、处理和基于合理监测数据的反演分析是评价其稳定特征的重要手段。变形监测一般采用在典型剖面选择典型高程位置进行布设,有限测点或者监测断面在分析整体变形模式和稳定特征时通常具有其局限性。通过有限点的监测数据获得典型断面的变形云图甚至获得整个建筑物三维空间的变形分布特征,快速实现实际监测数据的空间展示、以及与数值反演分析的变形云图进行比较分析对判断围岩稳定而言显得尤为重要。
发明内容
本发明的目的在于,针对现有技术中存在的不足,提供一种基于离散光滑插值技术的变形监测数据处理及分析方法。
为此,本发明的上述目的通过以下技术方案来实现:
一种基于离散光滑插值技术的变形监测数据处理方法,所述基于离散光滑插值技术的变形监测数据处理方法包括如下步骤:
(1)在三维插值技术分析软件中,根据所要分析的开挖体型轮廓建立几何模型,该几何模型与要分析的目标对象及其所处的施工阶段的体型轮廓一致,并初步划分网格单元,网格单元大小遵循开挖面附近较小且向几何模型外边界逐渐增大的原则,并根据监测点点位的疏密程度对网格单元大小进行调整,目的是使得划分后的网格单元在开挖面附近和监测点附近具有更好的插值结果,使得云图过渡更加平顺;
(2)将工程现场已经安装并已经获得数据的变形监测仪器的具***置和变形量/变形增量形成一个多维数组,并且变形监测仪器的位置坐标和步骤(1)中网格化的模型具有相同的坐标系,并可以在网格化的模型中标定出变形监测仪器的具***置,也即是监测点的具***置;
(3)将步骤(2)中所形成的带有变形属性特征的监测点的数组集合导入至步骤(1)中所建立好的网格化模型中,根据分析需求将所有监测点的变形量/增量或者其中一个变形分量设置为插值的约束点;
(4)通过离散光滑插值技术,给定网格单元相应的初始值以及将(3)中的约束点的变形量设置为约束插值条件,可以直接插值得到整个网格模型的变形场分布特征。
在采用上述技术方案的同时,本发明还可以采用或者组合采用以下进一步的技术方案:
优选地,只需要建立与实际分析对象轮廓相符的几何模型。
优选地,在建立好的几何轮廓模型基础上可以直接划分单元网格模型,网格大小可以任意指定,并且网格尺寸的大小可以根据分析对象的需要,可以在局部进行加密调整。
优选地,将监测数据包含的空间位置坐标以及监测成果写成多维数组的格式直接导入,可以实现数据的批量处理。
优选地,通过有限空间点的监测数据,可以插值得到空间范围内全局的变形分布特征,可以更加直观地用于与数值分析获得的变形特征对比。
本发明还有一个目的在于,针对现有技术中存在的不足,提供一种基于离散光滑插值技术的变形监测数据分析方法。
为此,本发明的上述目的通过以下技术方案来实现:
一种基于离散光滑插值技术的变形监测数据分析方法,所述基于离散光滑插值技术的变形监测数据分析方法基于前文所述的基于离散光滑插值技术的变形监测数据处理方法,并包括如下步骤:根据工程现场的围岩开挖响应,宏观分析模型的整体变形特征,并对比基于数值模型的反演分析获得变形场特征,判断基于实际监测数据的围岩变形模式和围岩变形规律的合理性。
在采用上述技术方案的同时,本发明还可以采用或者组合采用以下进一步的技术方案:
优选地,若经离散光滑插值技术获得后的变形场分析特征与现场实际的整体开挖响应、经验判断以及技术数值模型的反演分析结果存在矛盾,可以删除不合理监测点的变形,并作为新的约束条件,然后再次插值从而获得合理的变形场分析特征。
本发明提供一种基于离散光滑插值技术的变形监测数据处理及分析方法,通过有限点的变形监测数据插值获得典型断面和三维全局的变形场分析特征,并通过比对数值模型的反演分析而获得的变形场分布特征,从而实现从局部和全局的维度来分析评价监测数据本身的合理性以及地下洞室围岩、边坡的稳定特征。
附图说明
图1a示出了地下洞室的网格模型。
图1b示出了边坡的网格模型。
图2a为地下洞室的变形监测点布置与网格模型的对应关系图。
图2b为边坡的变形监测点布置与网格模型的对应关系图。
图3a示出了地下洞室基于离散光滑插值技术处理后的变形场特征。
图3b示出了边坡基于离散光滑插值技术处理后的变形场特征。
具体实施方式
参照附图和具体实施例对本发明作进一步详细地描述。
一种基于离散光滑插值技术的变形监测数据处理方法,包括如下步骤:
(1)在三维插值技术分析软件中,根据所要分析的开挖体型轮廓建立几何模型,参照图1a和1b,该几何模型与要分析的目标对象及其所处的施工阶段的体型轮廓一致,并初步划分网格单元,网格单元大小遵循开挖面附近较小且向几何模型外边界逐渐增大的原则,并根据监测点点位的疏密程度对网格单元大小进行调整,目的是使得划分后的网格单元在开挖面附近和监测点附近具有更好的插值结果,使得云图过渡更加平顺;
(2)将工程现场已经安装并已经获得数据的变形监测仪器的具***置和变形量/变形增量形成一个多维数组,并且变形监测仪器的位置坐标和步骤(1)中网格化的模型具有相同的坐标系,并可以在网格化的模型中标定出变形监测仪器的具***置,也即是监测点的具***置,参照图2a和2b;
(3)将步骤(2)中所形成的带有变形属性特征的监测点的数组集合导入至步骤(1)中所建立好的网格化模型中,根据分析需求将所有监测点的变形量/增量或者其中一个变形分量设置为插值的约束点;
(4)通过离散光滑插值技术,给定网格单元相应的初始值以及将(3)中的约束点的变形量设置为约束插值条件,可以直接插值得到整个网格模型的变形场分布特征,参照图3a和3b。
以图3a地下洞室的监测成果为例,一般在地下洞室的监测断面会在边顶拱布置多支多点式位移计,每支多点位移计会布置2个以上的测点,每个测点距离开挖面的距离不同,比如四点式多点位移计为例,假定四个测点距离开挖面的距离为1 m,2 m,4 m,6 m,当获得该断面7支多点位移计共28个监测点的数据后,可以将监测点的空间坐标(x,y,z)以及变形量以文本格式写成一个多维数组:
(x1, y1, z1,xdisp1,ydisp1,zdisp1,disp1)
(x2, y2, z2,xdisp2,ydisp2,zdisp2,disp2)
….
(xn, yn, zn,xdispn,ydispn,zdispn,dispn)
然后将这些数据导入划分好网格单元的模型中,根据需要对所有网格单元的变形属性特征赋初始值0或者其他数值,将导入的所有监测点数据中的一个变形值(xdisp,ydisp,zdisp,disp其中的任意一个)作为约束点,即在监测点位置的网格单元其变形值在插值过程中始终与监测点实测值保持一致,然后对整个模型范围内进行插值,最终得到基于有限插值点数据获得的全局变形分布特征,并在插值后得到的变形云图中得到变形等值线d1,d2,d3,…,通过对比基于实测数据和光滑离散插值技术获得的全局变形分布和数值分析成果得到的变形分布,实现从局部和全局的维度来分析评价监测数据本身的合理性以及地下洞室围岩、边坡的稳定特征。
上述具体实施方式用来解释说明本发明,仅为本发明的优选实施例,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改、等同替换、改进等,都落入本发明的保护范围。

Claims (7)

1.一种基于离散光滑插值技术的变形监测数据处理方法,其特征在于,所述基于离散光滑插值技术的变形监测数据处理方法包括如下步骤:
(1)在三维插值技术分析软件中,根据所要分析的开挖体型轮廓建立几何模型,该几何模型与要分析的目标对象及其所处的施工阶段的体型轮廓一致,并初步划分网格单元,网格单元大小遵循开挖面附近较小且向几何模型外边界逐渐增大的原则,并根据监测点点位的疏密程度对网格单元大小进行调整,目的是使得划分后的网格单元在开挖面附近和监测点附近具有更好的插值结果,使得云图过渡更加平顺;
(2)将工程现场已经安装并已经获得数据的变形监测仪器的具***置和变形量/变形增量形成一个多维数组,并且变形监测仪器的位置坐标和步骤(1)中网格化的模型具有相同的坐标系,并可以在网格化的模型中标定出变形监测仪器的具***置,也即是监测点的具***置;
(3)将步骤(2)中所形成的带有变形属性特征的监测点的数组集合导入至步骤(1)中所建立好的网格化模型中,根据分析需求将所有监测点的变形量/增量或者其中一个变形分量设置为插值的约束点;
(4)通过离散光滑插值技术,给定网格单元相应的初始值以及将(3)中的约束点的变形量设置为约束插值条件,可以直接插值得到整个网格模型的变形场分布特征。
2.根据权利要求1所述的基于离散光滑插值技术的变形监测数据处理方法,其特征在于,只需要建立与实际分析对象轮廓相符的几何模型。
3.根据权利要求1所述的基于离散光滑插值技术的变形监测数据处理方法,其特征在于,在建立好的几何轮廓模型基础上可以直接划分单元网格模型,网格大小可以任意指定,并且网格尺寸的大小可以根据分析对象的需要,可以在局部进行加密调整。
4.根据权利要求1所述的基于离散光滑插值技术的变形监测数据处理方法,其特征在于,将监测数据包含的空间位置坐标以及监测成果写成多维数组的格式直接导入,可以实现数据的批量处理。
5.根据权利要求1所述的基于离散光滑插值技术的变形监测数据处理方法,其特征在于,通过有限空间点的监测数据,可以插值得到空间范围内全局的变形分布特征,可以更加直观地用于与数值分析获得的变形特征对比。
6.一种基于离散光滑插值技术的变形监测数据分析方法,其特征在于,所述基于离散光滑插值技术的变形监测数据分析方法基于权利要求1所述的基于离散光滑插值技术的变形监测数据处理方法,并包括如下步骤:根据工程现场的围岩开挖响应,宏观分析模型的整体变形特征,并对比基于数值模型的反演分析获得变形场特征,判断基于实际监测数据的围岩变形模式和围岩变形规律的合理性。
7.根据权利要求6所述的基于离散光滑插值技术的变形监测数据分析方法,其特征在于,若经离散光滑插值技术获得后的变形场分析特征与现场实际的整体开挖响应、经验判断以及技术数值模型的反演分析结果存在矛盾,可以删除不合理监测点的变形,并作为新的约束条件,然后再次插值从而获得合理的变形场分析特征。
CN201910371971.XA 2019-05-06 2019-05-06 一种基于离散光滑插值技术的变形监测数据处理及分析方法 Pending CN110263365A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910371971.XA CN110263365A (zh) 2019-05-06 2019-05-06 一种基于离散光滑插值技术的变形监测数据处理及分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910371971.XA CN110263365A (zh) 2019-05-06 2019-05-06 一种基于离散光滑插值技术的变形监测数据处理及分析方法

Publications (1)

Publication Number Publication Date
CN110263365A true CN110263365A (zh) 2019-09-20

Family

ID=67914298

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910371971.XA Pending CN110263365A (zh) 2019-05-06 2019-05-06 一种基于离散光滑插值技术的变形监测数据处理及分析方法

Country Status (1)

Country Link
CN (1) CN110263365A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110704934A (zh) * 2019-10-12 2020-01-17 上海勘察设计研究院(集团)有限公司 一种基于BIM模型的Web端轻量化展示方法
CN110851984A (zh) * 2019-11-12 2020-02-28 长江存储科技有限责任公司 一种数据处理方法及装置
CN111122402A (zh) * 2019-12-10 2020-05-08 北京蛙鸣华清环保科技有限公司 基于离散监测点数据的污染物路况图生成方法和***
CN117095081A (zh) * 2023-10-11 2023-11-21 武汉大学 谷幅变形监测范围确定方法及***

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103884291A (zh) * 2014-03-27 2014-06-25 河海大学 基于nurbs参数曲面的建筑物表面柔性变形监测方法
CN106682325A (zh) * 2016-12-30 2017-05-17 武汉英思工程科技股份有限公司 基于bim的拱坝施工期变形场动态拟合分析方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103884291A (zh) * 2014-03-27 2014-06-25 河海大学 基于nurbs参数曲面的建筑物表面柔性变形监测方法
CN106682325A (zh) * 2016-12-30 2017-05-17 武汉英思工程科技股份有限公司 基于bim的拱坝施工期变形场动态拟合分析方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王长海等: "基于离散光滑插值的三维地质体构造网格模型", 《深圳大学学报理工版》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110704934A (zh) * 2019-10-12 2020-01-17 上海勘察设计研究院(集团)有限公司 一种基于BIM模型的Web端轻量化展示方法
CN110851984A (zh) * 2019-11-12 2020-02-28 长江存储科技有限责任公司 一种数据处理方法及装置
CN110851984B (zh) * 2019-11-12 2021-05-07 长江存储科技有限责任公司 一种数据处理方法及装置
CN111122402A (zh) * 2019-12-10 2020-05-08 北京蛙鸣华清环保科技有限公司 基于离散监测点数据的污染物路况图生成方法和***
CN111122402B (zh) * 2019-12-10 2022-06-07 北京蛙鸣华清环保科技有限公司 基于离散监测点数据的污染物路况图生成方法和***
CN117095081A (zh) * 2023-10-11 2023-11-21 武汉大学 谷幅变形监测范围确定方法及***
CN117095081B (zh) * 2023-10-11 2023-12-19 武汉大学 谷幅变形监测范围确定方法及***

Similar Documents

Publication Publication Date Title
CN110263365A (zh) 一种基于离散光滑插值技术的变形监测数据处理及分析方法
CN107503797B (zh) 基于3d仿真平台的矿井瓦斯浓度趋势预警***
CN102609986B (zh) 数字矿山中矿体建模***及其建模和检测方法
CN100478996C (zh) 水利水电工程地质信息的三维统一模型构建方法
CN102434210B (zh) 地下工程画像信息与监测信息安全监控的方法与***
CN108710732A (zh) 一种盾构刀具服役期磨损的可视化评估和预测方法
CN205262429U (zh) 一种施工现场土方量测定***
KR102496592B1 (ko) 토공 정보 관리 시스템, 방법, 및 상기 방법을 실행시키기 위한 컴퓨터 판독 가능한 프로그램을 기록한 기록 매체
CN106952003A (zh) 高地应力区层状岩体地下洞室围岩破坏类型预测方法
Thum et al. 2D and 3D GIS-based geological and geomechanical survey during tunnel excavation
CN103205972A (zh) 一种分析基坑变形与坑外地面沉降关系的方法
CN115510527B (zh) 基于安全度指标的隧洞围岩稳定性判别和量化评估方法
CN110705178A (zh) 一种基于机器学习的隧道/地铁施工全过程围岩变形动态预测方法
CN109885866A (zh) 深埋隧道围岩形变荷载计算方法
CN111986321A (zh) 三维地质建模方法、装置、设备及存储介质
CN115238553A (zh) 一种地埋管线渗漏浸蚀的危险区域划分方法和***
Díaz-Fernández et al. Computation of influence functions for automatic mining subsidence prediction
CN100560940C (zh) 渗流-应力耦合现场试验的卸荷方法
CN102877471B (zh) 一种垂直深基坑失稳塌方的预警方法
CN102609614A (zh) 一种针对复杂矿体的开挖分析方法
CN107507179A (zh) 基于gocad的岩土体量化分析方法
CN115880449A (zh) 基于隧道超前钻探数据建立三维地质模型的方法及***
Bandpey et al. Validation of 3D discrete fracture network model focusing on areal sampling methods-a case study on the powerhouse cavern of Rudbar Lorestan pumped storage power plant, Iran
CN104462649A (zh) 一种矿体块段模型储量的自动更新方法
CN105718625A (zh) 协调的智能建模器

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20190920