CN112508399A - 一种地下空间开发适宜性评价方法及系统 - Google Patents

一种地下空间开发适宜性评价方法及系统 Download PDF

Info

Publication number
CN112508399A
CN112508399A CN202011408822.5A CN202011408822A CN112508399A CN 112508399 A CN112508399 A CN 112508399A CN 202011408822 A CN202011408822 A CN 202011408822A CN 112508399 A CN112508399 A CN 112508399A
Authority
CN
China
Prior art keywords
evaluation
index
underground space
suitability
coordinate system
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.)
Granted
Application number
CN202011408822.5A
Other languages
English (en)
Other versions
CN112508399B (zh
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.)
CREEC Chongqing Survey Design and Research Co Ltd
Original Assignee
CREEC Chongqing Survey Design and Research 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 CREEC Chongqing Survey Design and Research Co Ltd filed Critical CREEC Chongqing Survey Design and Research Co Ltd
Priority to CN202011408822.5A priority Critical patent/CN112508399B/zh
Publication of CN112508399A publication Critical patent/CN112508399A/zh
Application granted granted Critical
Publication of CN112508399B publication Critical patent/CN112508399B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Development Economics (AREA)
  • Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Educational Administration (AREA)
  • Operations Research (AREA)
  • Marketing (AREA)
  • Game Theory and Decision Science (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提出了地下空间开发适宜性评价方法及系统，该方法为建立评价区域的地下空间三维体素坐标系∑0，在2000国家大地坐标系∑1中查得坐标系∑0的原点在坐标系∑1的坐标，建立坐标系∑0和坐标系∑1仿射变换关系；建立二维地表栅格坐标系∑s，在评价区域内形成栅格集合G，记栅格g(x,y)∈G，其地下空间体素构成集合Vxy；基于g(x,y)建立数据结构T(x,y,I)用以描述地下空间体素集合Vxy各评价指标数据；对评价指标进行敏感性分类，分为一般指标和敏感指标；对所述一般指标通过定量指标和定性指标进行综合评价,计算一般指标的评价值；提高评价值过低或极高的指标的权重，调整得到该评价单元的综合评价值，并根据该值划分地下空间适宜性评价等级，本方法实现地下空间开发适宜性三维精细化评价。

Description

一种地下空间开发适宜性评价方法及系统
技术领域
本发明涉及地质勘探处理领域，具体涉及一种地下空间开发适宜性评价方法及系统。
背景技术
近年来，伴随着城市社会和经济的快速发展，城市中心区的交通、环境、用地紧张等问题也越来越突出，开发利用城市地下空间已势在必行。由于城市地下空间一经开发便很难逆转，且同一平面范围内，地下不同深度可以安排不同的功能空间，所以开展山城市地下空间开发适宜性三维评价，对于其地下空间三维规划及开发利用，具有十分重要的作用。
目前，国内外学者已针对城市地下空间适宜性评价开展了一系列的研究，然而仍存在局限之处：
(1)目前地下空间适宜性评价多以大跨度分层的方式划分地下空间，其本质仍是二维分析，无法满足城市地下空间三维规划的需求；
(2)针对城市地下空间的三维评估分析多以大尺寸体素(棱长数十米) 作为评价基本单元体，无法满足复杂城市地下空间精细化分析的需求；
(3)评价指标权重及质量等级的确定过于主观，不能很好解决评价指标的模糊性及敏感性，评价指标权重的“木桶效应”等问题。
因此，复杂城市地下空间的适宜性评价存在不准确、不精细的问题。
发明内容
为了克服上述现有技术中存在的缺陷，本发明的目的是提供一种地下空间开发适宜性评价方法及系统。
为了实现本发明的上述目的，本发明提供了一种地下空间开发适宜性评价方法，包括以下步骤：
S1,建立评价区域的地下空间三维体素坐标系∑0，并将坐标系∑0的X、Y、 Z方向分别以Rx、Ry和Rz的单位长度进行划分，在评价空间区域内形成空间体素集合V，记v(x,y,z)为空间体素集合V中的体素，记s为v对应的地表体素，则 s的坐标为(x,y,0)，即s是v在Z方向的原点，s(x,y,0)地表体素构成集合S；
S2,在2000国家大地坐标系∑1中查得坐标系∑0的原点在坐标系∑1的坐标，并记为(x0,y0,z0)，建立坐标系∑0和坐标系∑1仿射变换关系；
S3,建立二维地表栅格坐标系∑s，该坐标系∑s与坐标系∑0的原点、X方向、 Y方向以及各方向划分的单位长度均相同，在评价区域内形成栅格集合G，记栅格g(x,y)∈G，则s(x,y,0)∈S是g(x,y)对应的三维地表体素，其地下空间体素构成集合Vxy；
S4,基于g(x,y)建立数据结构T(x,y,I)用以描述地下空间体素集合Vxy各评价指标数据，其中x和y描述地表位置信息，I为城市地下空间开发适宜性评价指标集合，确定坐标系∑0中体素单元v(x,y,z)的地下空间适宜性评价指标Ii的属性值p,p值即为该体素单元参与评价运算的指标数据值；
S5,建立适宜性评价指标体系：对评价指标进行敏感性分类，分为一般指标和敏感指标；
S6,确定一般指标常权权重；
S7，结合指标数据值，对所述一般指标通过定量指标和定性指标进行综合评价，计算一般指标的评价值；
S8,提高评价值过低或极高的指标的权重，调整得到该评价单元的综合评价值，并根据该值划分地下空间适宜性评价等级。
该地下空间开发适宜性评价方法中评价指标权重及质量等级的确定客观，解决了评价指标的模糊性及敏感性以及评价指标权重的“木桶效应”等问题，最终实现城市地下空间开发适宜性三维精细化评价。
该评价方法的优选方案：所述地下空间三维体素坐标系∑0的X轴指向 BIH1984.0定义的零子午面与协议赤道的交点，Z轴指向历元2000.0的地球参考极的方向，Y轴按右手坐标系确定，以评价区域地表西南边界上的一点作为 X和Y方向原点，取地表体素中心点作为Z方向原点，并将坐标系的X、Y、Z 方向分别以Rx、Ry和Rz的单位长度进行划分。
该评价方法的优选方案：所述坐标系∑0和坐标系∑1仿射变换关系为：
地下空间点P在∑1中的坐标是(xp,yp,zp)，对应的地表点坐标是(xp,yp,zs)，则点P在∑0中的坐标(x,y,z)计算公式为：
Figure RE-GDA0002908622460000031
其中，floor是向下取整函数。
将坐标系∑0和坐标系∑1中的体素对应起来。
该评价方法的优选方案：所述步骤S4采用链式存储结构存储评价指标集合，根据评价需求调整评价指标体系：
对于每一评价指标Ii，T(x,y,Ii)采用链式存储结构存放三元组集合，每个三元组[t,b,p]以游程编码的方式描述地下空间体素集合Vxy中一个深度区间的评价指标数据，其中t为区间顶部埋深、b为区间底部埋深、p为该深度区间的评价指标属性；
确定坐标系∑0中体素单元v(x,y,z)的地下空间适宜性评价指标Ii的属性值p, p＝{pk|z·Rz∈[tk,bk)}，p值即为该体素单元参与评价运算的指标数据值，其中Rz为 Z轴划分的单位长度。
该评价方法的优选方案：步骤S7的具体方法为：
将定量指标中对评价结果影响方向为优的作为正向指标，对评价结果影响方向为劣的作为反向指标，对正向指标和反向指标进行标准化处理，具体如下：
正向指标标准值
Figure RE-GDA0002908622460000041
反向指标标准值
Figure RE-GDA0002908622460000042
其中：x为评价指标数据原始值；amax、amin分别为各原始评价等级区间的最大值及最小值，di为各评价等级的区间标准长度，D 为各区间标准长度之和；
利用定量指标数据标准值计算该评价指标关于四个评价等级的隶属关系，计算函数为：
Figure RE-GDA0002908622460000043
Figure RE-GDA0002908622460000044
Figure RE-GDA0002908622460000045
Figure RE-GDA0002908622460000046
其中，当该定量指标为正向指标时，x0为正向指标标准值x01，当该定量指标为反向指标时，x0为反向指标标准值x02；
对隶属函数计算结果进行归一化，得到定量指标隶属度μⅣ、μⅢ、μⅡ、μⅠ, 如下式：
Figure RE-GDA0002908622460000051
其中，(i＝Ⅰ,Ⅱ,Ⅲ,Ⅳ)；对于定性指标，采用模糊定性方法确定定性指标隶属度；计算单个一般指标的评价值
Figure RE-GDA0002908622460000052
该步骤对评价指标进行了有效的模糊评价。
该评价方法的优选方案：步骤S8的具体方法为：
计算评价单元适宜性综合得分，计算公式为：评价单元适宜性综合得分
Figure RE-GDA0002908622460000053
其中：
Figure RE-GDA0002908622460000054
为变权综合决策模型， Wj(X')＝wj(x1',...,xm')是与指标评价值有关的变权向量；j＝1,2,…,m，X'＝(x1',...,xm')为步骤S6得到的评价指标评价值向量；W＝(w1,...,wm)为评价指标常权向量； S(X')＝(S1(X'),…,Sm(X'))为评价指标状态变权向量；bs为与敏感指标相关的权重参数，若计算单元体内敏感指标的属性值p＝1，则bs＝0；若计算单元体内敏感指标的属性值p＝0，则bs＝1；
状态变权向量通过构造惩罚-激励协调均衡函数来获得，具体如下：
Figure RE-GDA0002908622460000055
其中，0<xj'≤0.1为重度惩罚区，0.1<xj'≤0.3为中度惩罚区，0.3<xj'≤0.4为轻度惩罚区，0.4<xj'≤0.8为合格区，0.8<xj'<1为激励区；
根据评价单元适宜性综合得分M划分地下空间适宜性评价等级。
这里通过建立惩罚-激励复合型变权综合决策模型，通过“惩罚”及“激励”的方式提高评价值过低或极高的因子的权重，从而调整该评价单元的综合评价值，使其更为全面地反应其适宜性程度。
本发明还提供了一种地下空间开发适宜性评价系统，包括存储器、处理器和显示器，所述处理器和所述存储器通过总线完成相互间的通信，所述处理器与显示器通信连接，所述存储器存储有可被所述处理器执行的程序指令，所述处理器调用所述程序指令能够执行如上述的方法对地下空间开发适宜性进行评价，所述处理器将评价结果于显示器上进行可视化展视。
本发明的有益效果是：本发明将多类别评价指标数据压缩至二维地表栅格，并支持评价过程中在体素空间内的三维重现，从而实现城市地下空间评价基本单元体的精细化划分；本发明建立地下空间三维适宜性评价指标体系，将敏感指标模型、模糊综合评价法以及变权思想进行优化复合，建立敏感指标- 模糊-变权复合评价计算模型，解决了评价指标权重及质量等级的确定过于主观，不能很好解决评价指标的模糊性及敏感性以及评价指标权重的“木桶效应”等问题，最终实现城市地下空间开发适宜性三维精细化评价。
本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1是本发明地下空间开发适宜性三维精细化评价流程图；
图2是本发明地下空间开发适宜性评价计算模型步骤流程图；
图3是基于二维地表栅格的数据结构图；
图4是指标要素栅格化实例图；
图5是本发明具体实施方式的运行结果。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。
在本发明的描述中，除非另有规定和限定，需要说明的是，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是机械连接或电连接，也可以是两个元件内部的连通，可以是直接相连，也可以通过中间媒介间接相连，对于本领域的普通技术人员而言，可以根据具体情况理解上述术语的具体含义。
如图1所示，本发明提供了一种地下空间开发适宜性评价方法，尤其适用于城市地下空间三维规划。
该方法具体步骤如下：
S1,建立评价区域的地下空间三维体素坐标系∑0，并将坐标系∑0的X、Y、 Z方向分别以Rx、Ry和Rz的单位长度进行划分，在评价空间区域内形成空间体素集合V，记v(x,y,z)为空间体素集合V中的体素，记s为v对应的地表体素，则 s的坐标为(x,y,0)，即s是v在Z方向的原点，s(x,y,0)地表体素构成集合S。
这里地下空间三维体素坐标系∑0的X轴指向BIH1984.0定义的零子午面与协议赤道的交点，Z轴指向历元2000.0的地球参考极的方向，Y轴按右手坐标系确定，以评价区域地表西南边界上的一点作为X和Y方向原点，取地表体素中心点作为Z方向原点，并将坐标系的X、Y、Z方向分别以Rx、Ry和Rz的单位长度进行划分。
在建立评价区域的地下空间三维体素坐标系∑0之前，可通过地质勘查工具采集、测量得到评价区域的地表尺寸以及地下尺寸等，也可以通过地图或其它地质资料得到评价区域的地表尺寸以及地下尺寸等。
S2,在2000国家大地坐标系∑1中查得坐标系∑0的原点在坐标系∑1的坐标，并记为(x0,y0,z0)，建立坐标系∑0和坐标系∑1仿射变换关系，具体如下：
坐标系∑0和坐标系∑1仿射变换关系为：
地下空间点P在∑1中的坐标是(xp,yp,zp)，对应的地表点坐标是(xp,yp,zs)，则点P在∑0中的坐标(x,y,z)计算公式为：
Figure RE-GDA0002908622460000081
其中，floor是向下取整函数。
S3,建立二维地表栅格坐标系∑s，该坐标系∑s与坐标系∑0的原点、X方向、 Y方向以及各方向划分的单位长度均相同，在评价区域内形成栅格集合G，记栅格g(x,y)∈G，则s(x,y,0)∈S是g(x,y)对应的三维地表体素，其地下空间体素构成集合Vxy。这里即划分好了三维精细化评价基本单元体素。
S4,基于g(x,y)建立数据结构T(x,y,I)用以描述地下空间体素集合Vxy各评价指标数据，其中x和y描述地表位置信息，I为城市地下空间开发适宜性评价指标集合，确定坐标系∑0中体素单元v(x,y,z)的地下空间适宜性评价指标Ii的属性值p,p值即为该体素单元参与评价运算的指标数据值。
具体的，采用链式存储结构存储评价指标集合，如图3所示，根据评价需求调整评价指标体系：
对于每一评价指标Ii，T(x,y,Ii)采用链式存储结构存放三元组集合，每个三元组[t,b,p]以游程编码的方式描述地下空间体素集合Vxy中一个深度区间的评价指标数据，其中t为区间顶部埋深、b为区间底部埋深、p为该深度区间的评价指标属性；
确定坐标系∑0中体素单元v(x,y,z)的地下空间适宜性评价指标Ii的属性值p, p＝{pk|z·Rz∈[tk,bk)}，其中Rz为Z轴划分的单位长度。
通俗的讲，当评价三维体素空间的某体素单元v(x,y,z)时，找到其地表体素 s(x,y,0)对应的二维地表栅格g(x,y)，则T(x,y,I)存储了s(x,y,0)对应地下空间体素集合Vxy的各评价指标数据，对于体素单元v(x,y,z)的某一评价指标Ii，遍历T(x,y,Ii) 存储的三元组集合，记[tk,bk,pk]为第k个三元组，则可以根据p＝{pk|z·Rz∈[tk,bk)}来确定指标属性，即v(x,y,z)的评价指标Ii的属性值为p。
以部分评价指标为例进行说明，将评价区域内建设现状指标的栅格图和工程地质的结构化数据表转换为二维地表栅格数据结构，再按照上述映射关系还原至地下三维体素空间。从图5可以看出，此数据模型实现了多类别评价指标数据的压缩和二维地表栅格数据结构到地下体素空间的三维重现，并将多类别评价指标数据统一集成在同一空间坐标系内。
S5,建立适宜性评价指标体系：对评价指标进行敏感性分类，分为一般指标和敏感指标。
其中一般指标有一级一般指标和二级一般指标，敏感指标有一级敏感指标和二级敏感指标。其中一级一般指标包括区域地质A、环境地质B、岩体性质 C、水文地质D、地形地貌E、建设现状F、区位价值G；二级一般指标是对一级一般指标的细化，包括地震震级A1、洪水A2、地质构造A3、地质灾害B1、绿地B2、岩石饱和抗压强度C1、岩体完整性C2、松散岩类裂隙水D1、基岩裂隙水D2、地面标高E1、地面坡度E2、文物保护区F1、广场F2、道路F3、区域GDPG1、区域人口G2、土地等级G3；一级敏感指标包括地下空间开发现状H；二级敏感指标是对一级敏感指标的细化，包括基础及其影响范围H1、既有地下建筑物H2、城市轨道交通H3、既有隧道H4、人防硐室H5。
本实施例中，可在地下设置与指标相对应的传感器或检测装置来采集一般指标和/或敏感指标，包括一般指标的定性指标和/或定量指标，也可以通过网络在地质局、城市开发规划研究院等机构公开的信息中获取一般指标和/或敏感指标，包括一般指标的定性指标和/或定量指标。
S6，利用层次分析法确定一般指标常权权重。
层次分析法基本流程如下：
(1)建立层次结构模型：将决策的目标、对象及纳入考虑的因素按照其相互关系分为最高层、中间层和最低层。
(2)构造判断矩阵：采用一致矩阵法，即将各个因素两两比较，得出相对权重，以提高多目标条件下的准确度。这样得到的矩阵就称为判断矩阵，判断矩阵中的量化值与对应因素相对重要性的关系如下表：
层次分析法判断矩阵量化值含义
因素i与因素j相比 量化值
同等重要 1
稍微重要 3
明显重要 5
强烈重要 7
极端重要 9
相邻重要性之间的过渡值 2,4,6,8
(3)一致性和随机性检验：先假设有一同阶正则向量A,使得存在 XA＝λmaxA,X为上述的判断矩阵，解此特征方程所得的A经正规化后即为各评价因子的权值。由于客观事物的复杂及对事物认识的片面性，构造的判断矩阵不一定是一致性矩阵(也不强求是一致性矩阵)，但当偏离一致性过大时,会导致一些问题的产生。因此得到λmax后,还需进行一致性和随机性检验。检验公式为：
Figure RE-GDA0002908622460000111
式中，CI为一致性指标；λmax为最大特征根，n为矩阵阶数；RI为平均随机一致性指标，取值见下表：层次分析法判断矩阵随机一致性指标标准值，CR 为随机一致性比率，只有当CR<0.10时，判断矩阵才具有满意的一致性，所获取值才比较合理。
层次分析法判断矩阵随机一致性指标标准值
Figure RE-GDA0002908622460000112
S7，单项指标模糊评价：经过敏感性分类后的一般指标包括定量指标和定性指标，因此结合指标数据值，对一般指标通过定量指标和定性指标进行综合评价,计算一般指标的评价值,其中，定量指标是可以准确数量定义、精确衡量的评价指标，如地震震级，地面标高等。定性指标：指不能直接量化的评价指标，如岩体完整性，土地等级等。
具体方法为：
将定量指标中对评价结果影响方向为优(即指标数据值p越大，评价结果越适宜)的作为正向指标，对评价结果影响方向为劣(即指标数据值p越大，评价结果越不适宜)的作为反向指标，对正向指标和反向指标进行标准化处理，具体如下：
正向指标标准值
Figure RE-GDA0002908622460000113
反向指标标准值
Figure RE-GDA0002908622460000114
其中：x为评价指标数据原始值；amax、amin分别为各原始评价等级区间的最大值及最小值，di为各评价等级的区间标准长度，D 为各区间标准长度之和；
利用定量指标数据标准值计算该评价指标关于四个评价等级的隶属关系，计算函数为：
Figure RE-GDA0002908622460000121
Figure RE-GDA0002908622460000122
Figure RE-GDA0002908622460000123
Figure RE-GDA0002908622460000124
其中，当该定量指标为正向指标时，x0为正向指标标准值x01，当该定量指标为反向指标时，x0为正向指标标准值x02；
对隶属函数计算结果进行归一化，得到定量指标隶属度μⅣ、μⅢ、μⅡ、μⅠ, 如下式：
Figure RE-GDA0002908622460000125
其中，(i＝Ⅰ,Ⅱ,Ⅲ,Ⅳ)，；
对于定性指标，采用模糊定性方法确定定性指标隶属度。本实施例中采用的模糊定性方法为根据定性指标适宜性评价等级，查询隶属度分布表如下，确定定性指标隶属度。
定性指标隶属度分布表
Figure RE-GDA0002908622460000126
计算单个一般指标的评价值
Figure RE-GDA0002908622460000131
如果是定量指标就代入定量指标的隶属度，如果是定性指标就代入定性指标的隶属度计算。
S8,提高评价值过低或极高的指标的权重，调整得到该评价单元的综合评价值，并根据该值划分地下空间适宜性评价等级。
本实施例中采用惩罚-激励型变权实现：采用惩罚-激励复合型变权综合决策模型，通过“惩罚”及“激励”的方式提高评价值过低或极高的指标的权重，调整得到该评价单元的综合评价值，使其更为全面地反应其适宜性程度，并根据该值划分地下空间适宜性评价等级。通常可设置评价值区间，当评价值低于评价值区间最低值时即为评价值过低，评价值高于评价值区间最高值时即为极高。
步骤S8的具体方法为：
计算评价单元适宜性综合得分，计算公式为评价单元适宜性综合得分
Figure RE-GDA0002908622460000132
其中：
Figure RE-GDA0002908622460000133
为变权综合决策模型， Wj(X')＝wj(x1',...,xm')是与指标评价值有关的变权向量；j＝1,2,…,m，X'＝(x1',...,xm')为经过S5模糊化后的评价指标评价值向量；W＝(w1,...,wm)为评价指标常权向量； S(X')＝(S1(X'),…,Sm(X'))为评价指标状态变权向量；bs为与敏感指标相关的权重参数，由于敏感指标会直接影响评价的不利结果，在评价单元内，敏感指标只有存在和不存在两种状态，存在直接导致不利，因此若计算单元体内敏感指标的属性值p为存在(＝1)，则bs＝0；若计算单元体内敏感指标的属性值p为不存在(＝0)，则bs＝1。
状态变权向量通过构造惩罚-激励协调均衡函数来获得，具体如下：
Figure RE-GDA0002908622460000141
其中，0<xj'≤0.1为重度惩罚区，0.1<xj'≤0.3为中度惩罚区，0.3<xj'≤0.4为轻度惩罚区，0.4<xj'≤0.8为合格区，0.8<xj'<1为激励区。
根据评价单元适宜性综合得分M划分地下空间适宜性评价等级。
本实施例中，地下空间开发适宜性评价结果共分为4个评价等级，分别为适宜(Ⅰ级)，较适宜(Ⅱ级)，较不适宜(Ⅲ级)，不适宜(Ⅳ级)，M≤2.5 为不适宜区，2.5<M≤5为较不适宜区，5<M≤7.5为较适宜区，7.5<M≤10 为适宜区。
城市地下空间开发时，可根据该评价结果进行开发。
根据三维评价单元体的计算结果数据在软件中进行可视化展视，如图5所示。
本实施例中，根据三维评价单元体的计算结果数据，在Cesium平台中实现三维可视化，如图5所示，对于任意深度，可在ArcMap中绘制二维评价结果。对比图例，地下空间工程建设应在开发较适宜、适宜的区域进行。根据适宜性三维评价结果，可选取指定深度，在ArcGIS中绘制开发适宜性二维评价结果，分析该深度地下空间开发适宜性平面分布规律。
在使用该地下空间开发适宜性评价方法时，还包括建立地下空间开间开发适宜性评价指标体系的步骤，该体系中包含评价指标要素，将所述评价指标要素标准化为shp文件，所述shp文件属性包含：区间顶部埋深t、区间底部埋深b、区间属性p；对所述评价指标要素进行栅格化，并以评价指标要素的区间顶部埋深t、区间底部埋深b、区间属性p作为特征字段输出栅格数据，以图层的方式添加至ArcMap中，并统一坐标系为2000国家大地坐标系，得到包含所有要素的初始条件。
对评价初始条件中的所有指标要素，使用要素转栅格工具(Feature to Raster)，选择输入要素、特征字段、输出位置及像元大小，对所有指标要素进行栅格化操作，以要素的区间顶部埋深t、区间底部埋深b、区间属性p作为特征字段输出栅格数据。如H1基础及其影响范围为例，如下图4所示。
本发明还提出了一种地下空间开发适宜性评价系统，包括存储器和处理器，所述处理器和所述存储器通过总线完成相互间的通信，所述存储器存储有可被所述处理器执行的程序指令，所述处理器调用所述程序指令能够执行上述的方法对地下空间开发适宜性进行评价。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。

Claims (10)

1.一种地下空间开发适宜性评价方法，其特征在于，包括以下步骤：
S1,建立评价区域的地下空间三维体素坐标系∑0，并将坐标系∑0的X、Y、Z方向分别以Rx、Ry和Rz的单位长度进行划分，在评价空间区域内形成空间体素集合V，记v(x,y,z)为空间体素集合V中的体素，记s为v对应的地表体素，则s的坐标为(x,y,0)，即s是v在Z方向的原点，s(x,y,0)地表体素构成集合S；
S2,在2000国家大地坐标系∑1中查得坐标系∑0的原点在坐标系∑1的坐标，并记为(x0,y0,z0)，建立坐标系∑0和坐标系∑1仿射变换关系；
S3,建立二维地表栅格坐标系∑s，该坐标系∑s与坐标系∑0的原点、X方向、Y方向以及各方向划分的单位长度均相同，在评价区域内形成栅格集合G，记栅格g(x,y)∈G，则s(x,y,0)∈S是g(x,y)对应的三维地表体素，其地下空间体素构成集合Vxy；
S4,基于g(x,y)建立数据结构T(x,y,I)用以描述地下空间体素集合Vxy各评价指标数据，其中x和y描述地表位置信息，I为城市地下空间开发适宜性评价指标集合，确定坐标系∑0中体素单元v(x,y,z)的地下空间适宜性评价指标Ii的属性值p,p值即为该体素单元参与评价运算的指标数据值；
S5,建立适宜性评价指标体系：对评价指标进行敏感性分类，分为一般指标和敏感指标；
S6,确定一般指标常权权重；
S7，结合指标数据值，对所述一般指标通过定量指标和定性指标进行综合评价，计算一般指标的评价值；
S8,提高评价值过低或极高的指标的权重，调整得到该评价单元的综合评价值，并根据该值划分地下空间适宜性评价等级。
2.根据权利要求1所述的地下空间开发适宜性评价方法，其特征在于，所述地下空间三维体素坐标系∑0的X轴指向BIH1984.0定义的零子午面与协议赤道的交点，Z轴指向历元2000.0的地球参考极的方向，Y轴按右手坐标系确定，以评价区域地表西南边界上的一点作为X和Y方向原点，取地表体素中心点作为Z方向原点，并将坐标系的X、Y、Z方向分别以Rx、Ry和Rz的单位长度进行划分。
3.根据权利要求1所述的地下空间开发适宜性评价方法，其特征在于，所述坐标系∑0和坐标系∑1仿射变换关系为：
地下空间点P在∑1中的坐标是(xp,yp,zp)，对应的地表点坐标是(xp,yp,zs)，则点P在∑0中的坐标(x,y,z)计算公式为：
Figure FDA0002816537330000021
其中，floor是向下取整函数。
4.根据权利要求1所述的地下空间开发适宜性评价方法，其特征在于，所述步骤S4采用链式存储结构存储评价指标集合，根据评价需求调整评价指标体系：
对于每一评价指标Ii，T(x,y,Ii)采用链式存储结构存放三元组集合，每个三元组[t,b,p]以游程编码的方式描述地下空间体素集合Vxy中一个深度区间的评价指标数据，其中t为区间顶部埋深、b为区间底部埋深、p为该深度区间的评价指标属性；
确定坐标系∑0中体素单元v(x,y,z)的地下空间适宜性评价指标Ii的属性值p,p＝{pk|z·Rz∈[tk,bk)}，p值即为该体素单元参与评价运算的指标数据值，其中Rz为Z轴划分的单位长度。
5.根据权利要求1所述的地下空间开发适宜性评价方法，其特征在于，步骤S7的具体方法为：
将定量指标中对评价结果影响方向为优的作为正向指标，对评价结果影响方向为劣的作为反向指标，对正向指标和反向指标进行标准化处理，具体如下：
正向指标标准值
Figure FDA0002816537330000031
反向指标标准值
Figure FDA0002816537330000032
其中：x为评价指标数据原始值；amax、amin分别为各原始评价等级区间的最大值及最小值，di为各评价等级的区间标准长度，D为各区间标准长度之和；
利用定量指标数据标准值计算该评价指标关于四个评价等级的隶属关系，计算函数为：
Figure FDA0002816537330000033
Figure FDA0002816537330000034
Figure FDA0002816537330000035
Figure FDA0002816537330000036
其中，当该定量指标为正向指标时，x0为正向指标标准值x01，当该定量指标为反向指标时，x0为反向指标标准值x02；
对隶属函数计算结果进行归一化，得到定量指标隶属度μⅣ、μⅢ、μⅡ、μⅠ,如下式：
Figure FDA0002816537330000037
其中，(i＝Ⅰ,Ⅱ,Ⅲ,Ⅳ)，；
对于定性指标，采用模糊定性方法确定定性指标隶属度；
计算单个一般指标的评价值
Figure FDA0002816537330000041
6.根据权利要求1所述的地下空间开发适宜性评价方法，其特征在于，步骤S8的具体方法为：
计算评价单元适宜性综合得分，计算公式为：评价单元适宜性综合得分
Figure FDA0002816537330000042
其中：
Figure FDA0002816537330000043
为变权综合决策模型，Wj(X')＝wj(x1',...,xm')是与指标评价值有关的变权向量；j＝1,2,…,m，X'＝(x1',...,xm')为步骤S7得到的评价指标评价值向量；W＝(w1,...,wm)为评价指标常权向量；S(X')＝(S1(X'),…,Sm(X'))为评价指标状态变权向量；bs为与敏感指标相关的权重参数，若计算单元体内敏感指标的属性值p＝1，则bs＝0；若计算单元体内敏感指标的属性值p＝0，则bs＝1；
状态变权向量通过构造惩罚-激励协调均衡函数来获得，具体如下：
Figure FDA0002816537330000044
其中，0<xj'≤0.1为重度惩罚区，0.1<xj'≤0.3为中度惩罚区，0.3<xj'≤0.4为轻度惩罚区，0.4<xj'≤0.8为合格区，0.8<xj'<1为激励区；
根据评价单元适宜性综合得分M划分地下空间适宜性评价等级。
7.根据权利要求1所述的地下空间开发适宜性评价方法，其特征在于，还包括步骤S9，根据三维评价单元体的计算结果数据在软件中进行可视化展视。
8.根据权利要求1所述的地下空间开发适宜性评价方法，其特征在于，还包括建立地下空间开间开发适宜性评价指标体系的步骤，该体系中包含评价指标要素，将所述评价指标要素标准化为shp文件，所述shp文件属性包含：区间顶部埋深t、区间底部埋深b、区间属性p；对所述评价指标要素进行栅格化，并以评价指标要素的区间顶部埋深t、区间底部埋深b、区间属性p作为特征字段输出栅格数据。
9.根据权利要求1所述的地下空间开发适宜性评价方法，其特征在于，所述地下空间开间开发适宜性评价指标体系包括一级评价指标和二级评价指标，所述一级评价指标包括区域地质A、环境地质B、岩体性质C、水文地质D、地形地貌E、建设现状F、区位价值G、地下空间开发现状H；所述二级评价指标是对一级评价指标的细化，包括地震震级A1、洪水A2、地质构造A3、地质灾害B1、绿地B2、岩石饱和抗压强度C1、岩体完整性C2、松散岩类裂隙水D1、基岩裂隙水D2、地面标高E1、地面坡度E2、文物保护区F1、广场F2、道路F3、区域GDPG1、区域人口G2、土地等级G3、基础及其影响范围H1、既有地下建筑物H2、城市轨道交通H3、既有隧道H4、人防硐室H5。
10.一种地下空间开发适宜性评价系统，其特征在于，包括存储器、处理器和显示器，所述处理器和所述存储器通过总线完成相互间的通信，所述处理器与显示器通信连接，所述存储器存储有可被所述处理器执行的程序指令，所述处理器调用所述程序指令能够执行如权利要求1至9任一所述的方法对地下空间开发适宜性进行评价，所述处理器将评价结果于显示器上进行可视化展视。
CN202011408822.5A 2020-12-04 2020-12-04 一种地下空间开发适宜性评价方法及系统 Active CN112508399B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011408822.5A CN112508399B (zh) 2020-12-04 2020-12-04 一种地下空间开发适宜性评价方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011408822.5A CN112508399B (zh) 2020-12-04 2020-12-04 一种地下空间开发适宜性评价方法及系统

Publications (2)

Publication Number Publication Date
CN112508399A true CN112508399A (zh) 2021-03-16
CN112508399B CN112508399B (zh) 2023-11-14

Family

ID=74971739

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011408822.5A Active CN112508399B (zh) 2020-12-04 2020-12-04 一种地下空间开发适宜性评价方法及系统

Country Status (1)

Country Link
CN (1) CN112508399B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113610974A (zh) * 2021-07-28 2021-11-05 南京师范大学 城市地下空间资源的三维评价模型构建方法及装置
CN115030729A (zh) * 2022-04-30 2022-09-09 上海建科工程咨询有限公司 基于空间坐标系的钻爆法隧道围岩结构面数字化构建方法
CN115099675A (zh) * 2022-07-14 2022-09-23 海门市三德体育用品有限公司 一种健身场所布局规划方法
CN117611759A (zh) * 2023-11-30 2024-02-27 博雅达勘测规划设计集团有限公司 基于三维模型的评分地图生成方法、装置、终端及存储介质
CN118313695A (zh) * 2024-06-05 2024-07-09 广州海洋地质调查局三亚南海地质研究所 海洋空间适宜性评价方法、系统、电子设备及存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101659780B1 (ko) * 2016-02-23 2016-09-27 (주)이지스 3차원 지리 공간 시스템에서 2차원 도면 데이터를 지구 구면 공간에 표현하는 방법
CN106203760A (zh) * 2015-05-08 2016-12-07 中国人民解放军理工大学 一种城市地下空间规划方案评价方法
KR20170080256A (ko) * 2015-12-31 2017-07-10 한국국토정보공사 지하 공간 정보를 이용한 종합 관리 시스템
CN107730426A (zh) * 2017-09-12 2018-02-23 中山大学 一种基于ArcGIS的地下空间适宜性评价模型
CN109102180A (zh) * 2018-07-30 2018-12-28 北京大学 致密砂岩储层双甜点评价的综合参数评价方法
CN110472882A (zh) * 2019-08-21 2019-11-19 河南大学 基于主成分分析的城市开发用地适宜性评价方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106203760A (zh) * 2015-05-08 2016-12-07 中国人民解放军理工大学 一种城市地下空间规划方案评价方法
KR20170080256A (ko) * 2015-12-31 2017-07-10 한국국토정보공사 지하 공간 정보를 이용한 종합 관리 시스템
KR101659780B1 (ko) * 2016-02-23 2016-09-27 (주)이지스 3차원 지리 공간 시스템에서 2차원 도면 데이터를 지구 구면 공간에 표현하는 방법
CN107730426A (zh) * 2017-09-12 2018-02-23 中山大学 一种基于ArcGIS的地下空间适宜性评价模型
CN109102180A (zh) * 2018-07-30 2018-12-28 北京大学 致密砂岩储层双甜点评价的综合参数评价方法
CN110472882A (zh) * 2019-08-21 2019-11-19 河南大学 基于主成分分析的城市开发用地适宜性评价方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘顺昌;李黎;徐德馨;: "基于地质环境的地下空间开发利用研究", 城市勘测, no. 01, pages 198 - 202 *
张晶晶;马传明;匡恒;周爱国;夏友;: "郑州市地下空间开发地质环境适宜性变权评价", 水文地质工程地质, no. 02, pages 128 - 135 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113610974A (zh) * 2021-07-28 2021-11-05 南京师范大学 城市地下空间资源的三维评价模型构建方法及装置
CN115030729A (zh) * 2022-04-30 2022-09-09 上海建科工程咨询有限公司 基于空间坐标系的钻爆法隧道围岩结构面数字化构建方法
CN115099675A (zh) * 2022-07-14 2022-09-23 海门市三德体育用品有限公司 一种健身场所布局规划方法
CN117611759A (zh) * 2023-11-30 2024-02-27 博雅达勘测规划设计集团有限公司 基于三维模型的评分地图生成方法、装置、终端及存储介质
CN117611759B (zh) * 2023-11-30 2024-07-09 博雅达勘测规划设计集团有限公司 基于三维模型的评分地图生成方法、装置、终端及存储介质
CN118313695A (zh) * 2024-06-05 2024-07-09 广州海洋地质调查局三亚南海地质研究所 海洋空间适宜性评价方法、系统、电子设备及存储介质

Also Published As

Publication number Publication date
CN112508399B (zh) 2023-11-14

Similar Documents

Publication Publication Date Title
CN112508399A (zh) 一种地下空间开发适宜性评价方法及系统
Guo et al. Three-dimensional geological modeling and spatial analysis from geotechnical borehole data using an implicit surface and marching tetrahedra algorithm
CN113283802B (zh) 一种复杂艰险山区滑坡危险性评估方法
Zhao et al. Slope unit-based landslide susceptibility mapping using certainty factor, support vector machine, random forest, CF-SVM and CF-RF models
CN111898315B (zh) 基于分形—机器学习混合模型的滑坡易发性评估方法
Dou et al. 3D geological suitability evaluation for urban underground space development based on combined weighting and improved TOPSIS
CN113554345B (zh) 地震滑坡灾害区域及灾害链风险评估方法
CN112529267A (zh) 一种基于模糊综合评价法的超前地质预报方法及系统
CN112652066B (zh) 一种基于三维地质模型的地质表征情况的评价方法及系统
CN109118004B (zh) 一种工程构筑选址适宜区预测方法
Grunwald et al. Soil layer models created with profile cone penetrometer data
Su et al. Engineering application of fuzzy evaluation based on comprehensive weight in the selection of geophysical prospecting methods
CN114359713A (zh) 一种基于遥感地质条件的道路通行性分析方法
Dong et al. Visualization of foundation evaluation for urban rail transit based on CGB technology integration
Chitsazan et al. Land subsidence susceptibility mapping using PWRSTFAL framework and analytic hierarchy process: fuzzy method (case study: Damaneh-Daran Plain in the west of Isfahan Province, Iran)
Liu et al. Soil polygon disaggregation through similarity-based prediction with legacy pedons
Aldungarova et al. Geotechnical interpolation methodology for determining intermediate values of soil properties
Kim et al. Identification of geomorphological features of landslides using airborne laser altimetry
Han et al. Analysis on Influencing Factors of urban comprehensive carrying capacity: a case study of five provinces in Northwest China
CN108132486A (zh) 一种重磁梯度与地震数据联合界面反演的优化模拟退火法
CN118840003A (zh) 基于三维地质模型和层次评价法的地下环境定量评价方法
Li et al. Study on slope deformation partition and monitoring point optimization considering spatial correlation
Liu et al. Rainfall-induced Landslide Hazard Analysis Based on GIS in Yanbian Region, China
CN117541736A (zh) 一种扫描地质三维模型的地灾风险分布图生成方法及系统
Reis et al. Endokarst Potential Assessment in the Northern Sector of Santo António Plateau (Estremadura Limestone Massif, Central Portugal)

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