发明名称 一种表征地面沉降的方法
摘要 本发明涉及一种表征地面沉降的方法,所述指标体系包含地面沉降表征指标体系和综合沉降指数。地面沉降表征体系包括:点指标体系、线指标体系、面指标体系和体指标体系。点指标体系包括:沉降量、最大沉降量、沉降速率、最大沉降速率、累计沉降量、平均沉降量和矢量方向,用来表征地面沉降的量值和速率。线指标体系包括:沉降差、倾斜和不均匀沉降系数,用来表征不同点之间的差异沉降。面指标体系包括:沉降面积、沉降面积变化率和漏斗密度,用来表征地面沉降在平面上的展布规律和漏斗发育的密度。体指标体系包括:沉降漏斗体积、沉降漏斗体积变化率、极限沉降时间、单位沉降体积损失和漏斗类型,用来表征地面沉降的三维特征。综合沉降指数用来评价地面沉降的危险性。本发明可用于地面沉降信息发布、地面沉降控制、地面沉降危险性分析以及风险区划。
申请公布号 CN102494667B 申请公布日期 2013.04.10
申请号 CN201110388387.9 申请日期 2011.11.30
申请人 同济大学 发明人 王建秀;谷雪影;蒋玉坤;司鹏飞;吴远斌;吴治厚
分类号 G01C5/00(2006.01)I 主分类号 G01C5/00(2006.01)I
代理机构 上海正旦专利代理有限公司 31200 代理人 张磊
主权项 1.一种表征地面沉降的方法,其特征在于具体步骤为:(1)、搜集沉降观测数据城市一般设有多个地面沉降监测点,需要测量各监测点的高程(z)和平面坐标(x,y);搜集各观测点的高程进行初步分析得到沉降量、最大沉降量、累计沉降量、平均沉降量、沉降速率、最大沉降速率;分析各点的沉降量和平面坐标的变化,即可得其矢量方向;沉降量为各监测点相对上一周期的沉降量;最大沉降量为各监测点沉降量中的最大值;累计沉降量为各测点第一次沉降监测到本期沉降监测的沉降量的总和;平均沉降量为各监测点沉降量的平均值;沉降速率为各监测点的沉降量与沉降时间的比值;最大沉降速率为监测点沉降速率中的最大值;矢量方向则以该监测点该监测周期的初始位置为起点,监测周期结束时的位置为终点作一个向量,该向量为该观测点运动的矢量方向;(2)、根据点指标结合沉降等值线图得到线指标沉降差<img file="705577DEST_PATH_IMAGE001.GIF" wi="29" he="18" />:两个观测点1、2 之间的沉降差的计算公式:<img file="988659DEST_PATH_IMAGE002.GIF" wi="109" he="28" />(5)式中:<i>H</i><sub>1</sub>——观测点1的沉降量;<i>H</i><sub>2</sub>——观测点2的沉降量;倾斜<i>T</i>:指两点的沉降差与这两点水平距离L的比值,无量纲,计算公式:<img file="637946DEST_PATH_IMAGE003.GIF" wi="145" he="48" />(6)式中:<i>T</i><sub>2-1</sub>——观测点1和点2之间的倾斜;<i>H</i><sub>2</sub>——观测点2的沉降量;<i>H</i><sub>1</sub>——观测点1的沉降量;<img file="690085DEST_PATH_IMAGE004.GIF" wi="29" he="25" />——观测点1和点2的水平距离;不均匀沉降系数<img file="953576DEST_PATH_IMAGE005.GIF" wi="16" he="20" />:计算公式如下:<img file="474687DEST_PATH_IMAGE006.GIF" wi="69" he="28" />(7)<img file="911354DEST_PATH_IMAGE007.GIF" wi="113" he="49" /><img file="834310DEST_PATH_IMAGE008.GIF" wi="128" he="49" />式中:<i>n——</i>统计个数;<img file="336879DEST_PATH_IMAGE009.GIF" wi="24" he="20" />——相邻等值线间的沉降速率变化量,mm/a;<img file="345286DEST_PATH_IMAGE010.GIF" wi="25" he="18" />——相邻等值线间距,mm;(3)、根据地面沉降等值线图划分沉降漏斗得到面指标从地面沉降监测中心发布的地面沉降等值线图中划分地面沉降漏斗;圈出地面沉降漏斗,并以1,2,3……表示;将地面沉降漏斗导入CAD软件,并进行合适的缩放,描出每个地面沉降漏斗的轮廓,统计每个漏斗的面积、沉降区域的总面积;沉降漏斗面积变化率的计算可以通过漏斗的面积变化量与产生变化的时间之比;漏斗密度是漏斗的数目与研究区域的面积之比;(4)、将地面沉降等值线图导入GIS 软件,采用三维分析工具进行分析得到体指标沉降漏斗体积<i>V</i>:利用三维分析可以统计某一个平面以上或以下的体积;沉降漏斗体积即地表以下的体积;数学表达式:<img file="585644DEST_PATH_IMAGE011.GIF" wi="353" he="36" />(1)沉降漏斗体积变化率<img file="363107DEST_PATH_IMAGE012.GIF" wi="29" he="42" />:区域沉降漏斗经多期观测会得到不同的漏斗体积,漏斗体积的变化量与时间的比值即可得到沉降漏斗体积变化率;沉降漏斗体积变化率<img file="30717DEST_PATH_IMAGE012.GIF" wi="29" he="42" />:表示单位时间内一般为1年沉降漏斗体积的增长量,单位:m<sup>3</sup>/a,在空间上反应沉降漏斗的发展速度;极限沉降时间:海平面以上地层的体积与沉降漏斗体积变化率的比值;海平面以上地层体积可以由GIS三维分析得到;单位:a;其计算公式如下:<img file="260841DEST_PATH_IMAGE013.GIF" wi="65" he="65" />(2)式中:<img file="304890DEST_PATH_IMAGE014.GIF" wi="21" he="25" />——极限沉降时间;<img file="936859DEST_PATH_IMAGE015.GIF" wi="18" he="25" />——海平面以上地层的体积;<img file="713054DEST_PATH_IMAGE012.GIF" wi="29" he="42" />——沉降漏斗体积变化率;单位沉降体积损失<i>m</i>:统计每个沉降漏斗的体积和相应的最大沉降量,沉降漏斗体积与最大沉降量的比值即为单位沉降体积损失;单位:m<sup>3</sup>/mm;计算公式如下:<img file="945321DEST_PATH_IMAGE016.GIF" wi="66" he="46" />(3)式中:<img file="278214DEST_PATH_IMAGE017.GIF" wi="18" he="16" />——单位沉降体积损失;<i>V</i>——沉降漏斗体积;<i>W</i><sub>max</sub>——最大沉降量;沉降漏斗类型:若<i>m</i>≥<img file="13957DEST_PATH_IMAGE018.GIF" wi="18" he="25" />m<sup>3</sup>/mm为平缓型漏斗;若<img file="39682DEST_PATH_IMAGE019.GIF" wi="17" he="25" />≤<i>m&lt;</i><img file="493666DEST_PATH_IMAGE018.GIF" wi="18" he="25" />m<sup>3</sup>/mm为过渡型漏斗;若 <i>m</i>&lt;<img file="630250DEST_PATH_IMAGE019.GIF" wi="17" he="25" />m<sup>3</sup>/mm定义为急剧型漏斗;(5)、计算综合沉降指数并判断沉降危险性综合沉降指数计算公式如下:<img file="175761DEST_PATH_IMAGE020.GIF" wi="69" he="62" />(4)式中 <i>I</i>=综合沉降指数;<img file="621655DEST_PATH_IMAGE021.GIF" wi="17" he="25" />=第<i>i</i>项指标的取值;<img file="500618DEST_PATH_IMAGE022.GIF" wi="20" he="25" />=第<i>i</i>项指标的权重;根据公式(4)计算地面沉降综合指数,并且给出地面沉降危险性等级;在综合沉降指数的计算模型中,从点、线、面、体四个指标体系中各选出一个指标参与计算,选择累计沉降量、不均匀沉降系数、沉降面积和沉降漏斗体积变化率;第1项指标累计沉降量<i>H</i><sub>cum</sub>分级<img file="175313DEST_PATH_IMAGE023.GIF" wi="17" he="25" />取值方法如下:<i> H</i><sub><i>cum</i></sub><200mm,<img file="885649DEST_PATH_IMAGE023.GIF" wi="17" he="25" />=1;200mm≤<i>H</i><sub>cum</sub><400mm,<img file="253176DEST_PATH_IMAGE023.GIF" wi="17" he="25" />=3;400mm≤<i>H</i><sub>cum</sub><600mm,<img file="353856DEST_PATH_IMAGE023.GIF" wi="17" he="25" />=5;600mm≤<i>H</i><sub>cum</sub><800mm,<img file="81510DEST_PATH_IMAGE023.GIF" wi="17" he="25" />=7;<i>H</i><sub>cum</sub>≥800mm,<img file="849614DEST_PATH_IMAGE023.GIF" wi="17" he="25" />=9;第2项指标不均匀沉降系数<i>m</i>分级<img file="122464DEST_PATH_IMAGE024.GIF" wi="18" he="25" />取值方法如下:<i>m</i><2 mm /a<sup>.</sup>km,<img file="38336DEST_PATH_IMAGE024.GIF" wi="18" he="25" />=1;2 mm /a<sup>.</sup>km ≤ <i>m</i><3 mm /a<sup>.</sup>km,<img file="54834DEST_PATH_IMAGE024.GIF" wi="18" he="25" />=3;3 mm /a<sup>.</sup>km ≤ <i>m</i><4 mm /a<sup>.</sup>km,<img file="745621DEST_PATH_IMAGE024.GIF" wi="18" he="25" />=5;4 mm /a<sup>.</sup>km ≤ <i>m</i><5 mm /a<sup>.</sup>km,<img file="189372DEST_PATH_IMAGE024.GIF" wi="18" he="25" />=7;<i>m</i>≥5 mm /a<sup>.</sup>km,<img file="592540DEST_PATH_IMAGE024.GIF" wi="18" he="25" />=9;第3项指标沉降面积<i>S</i>分级<img file="412729DEST_PATH_IMAGE025.GIF" wi="18" he="25" />取值方法如下:<i>S</i><10 km<sup>2</sup>,<img file="952163DEST_PATH_IMAGE025.GIF" wi="18" he="25" />=1;10 km<sup>2</sup>≤ <i>S</i><50 km<sup>2</sup>,<img file="19345DEST_PATH_IMAGE025.GIF" wi="18" he="25" />=3;50 km<sup>2</sup>≤ <i>S</i><100 km<sup>2</sup>,<img file="394963DEST_PATH_IMAGE025.GIF" wi="18" he="25" />=5;100 km<sup>2</sup>≤ <i>S</i><500 km<sup>2</sup>,<img file="2531DEST_PATH_IMAGE025.GIF" wi="18" he="25" />=7;<i>S</i>≥500 km<sup>2</sup>,<img file="412784DEST_PATH_IMAGE025.GIF" wi="18" he="25" />=9;第4项指标沉降漏斗体积变化率<img file="447605DEST_PATH_IMAGE012.GIF" wi="29" he="42" />分级<img file="310518DEST_PATH_IMAGE026.GIF" wi="18" he="25" />取值方法如下:<img file="721777DEST_PATH_IMAGE012.GIF" wi="29" he="42" /><1×10<sup>6</sup>m<sup>3</sup>/a,<img file="986536DEST_PATH_IMAGE026.GIF" wi="18" he="25" />=1;1×10<sup>6</sup>m<sup>3</sup>/a≤<img file="457838DEST_PATH_IMAGE012.GIF" wi="29" he="42" /><4×10<sup>6</sup>m<sup>3</sup>/a,<img file="808048DEST_PATH_IMAGE026.GIF" wi="18" he="25" />=3;4×10<sup>6</sup>m<sup>3</sup>/a≤<img file="751559DEST_PATH_IMAGE012.GIF" wi="29" he="42" /><7×10<sup>6</sup>m<sup>3</sup>/a,<img file="870824DEST_PATH_IMAGE026.GIF" wi="18" he="25" />=5;7×10<sup>6</sup>m<sup>3</sup>/a ≤<img file="513027DEST_PATH_IMAGE012.GIF" wi="29" he="42" /><10×10<sup>6</sup>m<sup>3</sup>/a,<img file="537484DEST_PATH_IMAGE026.GIF" wi="18" he="25" />=7;<img file="41277DEST_PATH_IMAGE012.GIF" wi="29" he="42" />≥10×10<sup>6</sup>m<sup>3</sup>/a,<img file="264317DEST_PATH_IMAGE026.GIF" wi="18" he="25" />=9;四个指标的权重<img file="828154DEST_PATH_IMAGE022.GIF" wi="20" he="25" />由层次分析法得到,分别为:0.2513、0.2346、0.2393、0.2748;地面沉降危险性分级包括:非常危险、危险、中等、低、非常低五个等级;根据计算出的<i>I</i>值,判断地面沉降危险性分级,1.0≤<i>I</i>≤2.6,危险性非常低;2.6&lt;<i> I </i>≤4.2,危险性低;4.2&lt;<i> I </i>≤5.8,危险性中等;5.8&lt; <i>I </i>≤ 7.4 危险性高;7.4&lt;<i> I</i>≤9.0,危险性非常高。
地址 200092 上海市杨浦区四平路1239号