CN110288144B - 一种基于效用理论和d-s证据推理的深基坑安全性评估方法 - Google Patents

一种基于效用理论和d-s证据推理的深基坑安全性评估方法 Download PDF

Info

Publication number
CN110288144B
CN110288144B CN201910534172.XA CN201910534172A CN110288144B CN 110288144 B CN110288144 B CN 110288144B CN 201910534172 A CN201910534172 A CN 201910534172A CN 110288144 B CN110288144 B CN 110288144B
Authority
CN
China
Prior art keywords
evaluation
foundation pit
deep foundation
safety
evaluation index
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
CN201910534172.XA
Other languages
English (en)
Other versions
CN110288144A (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.)
Nanjing Deep Intelligent Construction Technology Research Institute Co ltd
Original Assignee
Nanjing Deep Intelligent Construction Technology Research Institute 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 Nanjing Deep Intelligent Construction Technology Research Institute Co ltd filed Critical Nanjing Deep Intelligent Construction Technology Research Institute Co ltd
Priority to CN201910534172.XA priority Critical patent/CN110288144B/zh
Publication of CN110288144A publication Critical patent/CN110288144A/zh
Application granted granted Critical
Publication of CN110288144B publication Critical patent/CN110288144B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0635Risk analysis of enterprise or organisation activities
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393Score-carding, benchmarking or key performance indicator [KPI] analysis

Landscapes

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

Abstract

本发明公开了一种基于效用理论和D‑S证据推理的深基坑安全性评估方法,该方法依托深基坑项目检测数据,基于三阶Volterra级数对深基坑未来状态进行预测,提高方案对项目潜在风险的预测能力;采用变异系数和信息熵对评估指标融合权重进行确定,考虑指标变化幅度和包含信息量对安全评估的表征作用,提高评估的准确性;利用效用理论和D‑S证据推理进行安全性评估,实现项目安全性的定性和定量分析。本方法利用有限长单位冲激响应滤波器对预测模型进行简化,减少模型系数,提高了状态预测速度。本发明提供的方法无需建立复杂的评估模型,且不需要大量的项目风险数据进行模型训练,具备良好的方案可实现性和评估准确性,无需太多专家人员的介入,方便后期方案的有效使用和推广。

Description

一种基于效用理论和D-S证据推理的深基坑安全性评估方法
技术领域
本发明涉及深基坑项目安全评估技术领域,具体涉及一种基于效用理论和D-S证据推理的深基坑安全性评估方法。
背景技术
深基坑在建设和后期使用过程中,具有安全因素复杂多变的特性,对其进行及时有效的安全性评估进而采取有效的应对措施,对实现深基坑安全风险的有效降低有着重要意义。由于深基坑存在一定的复杂性和差异性,目前常用的深基坑项目安全评估方法,大多根据专家经验或现场调查方法实现定性分析,结合部分定量分析方法,做出粗略的安全评估分级,或针对部分重要指标设置一定的安全报警阈值,避免发生严重事故。上述安全评估方法存在一定的主观性、不确定性和滞后性,难以实现及时有效的深基坑安全评估,无法为深基坑项目的智能管理和长期风险方法提供支持。
随着现代传感器检测、智能评估等领域的不断发展,为深基坑项目的安全性评估方案提供了新的思路和解决方案。根据深基坑项目特点,利用现代技术手段结合先进的安全评估方案,实现深基坑安全性及时有效评估,对提高项目安全建设和有效防范远期工程风险,及时采取有效应对措施有着重要的工程实践意义。
深基坑项目安全因素众多,难以建立复杂完整的物理评估模型,同时可用于安全性评估的样本数据较少,因此难以采用传统的贝叶斯网络、故障数、模糊理论以及需要大量样本数据的神经网络等方法实现准确有效的安全性评估。D-S证据推理理论结合了决策理论和证据理论,是一种有效的不确定推理方案,可以根据相关数据对复杂的评估对象进行有效的定性和定量分析。基于D-S证据推理对深基坑项目进行安全性评估,依托项目检测数据,无需建立复杂的评估模型,且不需要大量的项目风险数据进行模型训练,具备较高的方案可实施性。通过对深基坑未来状态信息进行预测,可以提高评估模型的预测性能,对未来可能出现的风险信息进行提前响应,及时预测风险。
发明内容
技术问题:本发明提供一种基于效用理论和D-S证据推理的深基坑安全性评估方法,在评估数据方面,利用历史和当前检测数据,基于三阶Volterra级数对深基坑未来状态进行预测,能够提高方案对项目潜在风险的预测能力。采用变异系数和信息熵对评估指标融合权重进行确定,有效考虑了指标变化幅度和包含信息量对安全评估的表征作用,提高评估的准确性。利用效用理论和D-S证据推理进行安全性评估,可以有效地实现项目安全性的定性和定量分析。
技术方案:一种基于效用理论和D-S证据推理的深基坑安全性评估方法,其特征在于,该方法包括以下步骤:通过深基坑现场传感器设备获取用于安全性评估的评估指标数据后,
步骤1:使用基于三阶Volterra级数建立的未来状态非线性预测模型,通过历史和当前数据对深基坑未来状态进行预测,最终得到状态预测值y(k):
Figure BSA0000184791190000011
步骤2:基于效用理论,对评估指标数据进行处理,得到评估指标关于各评估等级的置信分布S(zi):
S(zi)={(Fj,βj,i),j=1,…N}
其中,zi表示指标i的数值,Fj表示第j级评估等级,共N个评估等级。
步骤3:基于变异系数法确定评估指标融合权重,得到评估指标xi在k时刻的变异系数融合权重ω1i(k):
Figure BSA0000184791190000021
步骤4:基于信息熵法确定评估指标融合权重,得到评估指标xi在k时刻的信息熵融合权重ω2i(k):
Figure BSA0000184791190000022
步骤5:基于卡方距离实现评估指标组合赋权,得到评估指标最终融合权重向量ω:
ω=α1ω12ω2
步骤6:基于D-S证据推理建立评估指标融合模型,得到k时刻深基坑安全性评估结果的置信分布:
O(g(k))={(Fj,βj(k)),j=1,...,N}
其中,输出转换模型O(g(k))将深基坑安全性评估结果g(k)转换为置信分布形式,βj(k)表示k时刻深基坑安全性对于评估等级Fj的置信度。
步骤7:基于D-S证据推理建立多时刻深基坑安全性评估结果融合模型,融合历史时刻、当前时刻和未来预测时刻的安全性评估结果,得到k时刻深基坑综合评估结果的置信分布:
O(G(k))={(Fj,β′j(k)),j=1,...,N}
其中,βj′(k)表示k时刻下系统安全性对于等级Fj的置信度。
步骤8:基于期望效用理论,将深基坑安全性评估置信分布变换为安全性评估数值,进而得到评估等级。
进一步地,所述步骤1中的深基坑未来状态预测方法包括如下步骤:
步骤1.1:以深基坑安全性指标的历史数据和当前观测值作为模型输入,如下式所示:
Z(k)=[z(k),z(k-1),…,z(k-T+1)]T
其中,z(k)表示评估指标在k时刻的采集数值,T表示用于预测未来状态的时间序列长度。
步骤1.2:基于三阶Volterra级数建立未来状态预测模型,得到深基坑未来状态预测值y(k):
Figure BSA0000184791190000023
步骤1.3:引入有限长单位冲激响应滤波器,并利用乘积耦合减少预测模型系数,得到简化系数后的模型输出预测值y(k):
Figure BSA0000184791190000031
其中,hm,r表示滤波器系数。
步骤1.4:利用改进的非线性归一化LMS滤波算法(NNLMS),达到提高滤波器收敛速度的目的,得到近似的滤波器系数为:
Figure BSA0000184791190000032
Figure BSA0000184791190000033
Figure BSA0000184791190000034
其中,A(k),B(k),C(k)分别表示三个有限长单位冲激响应滤波器系数,Z′(k)=[1,z(k)...z(k-T+1)]T
Figure BSA0000184791190000035
A(k)=[a0,α1…aT]T,B(k)=[b0,b1…bT]T,C(k)=[c0,c1…cT]T,并且μ1≠μ2≠μ3
进一步地,所述步骤2中的评估指标关于评估等级置信分布求取过程包括如下步骤:
步骤2.1:针对深基坑安全性评估,定义N个安全性等级,得到深基坑评估等级的完备集F:
F={F1,…,Fj,…,FN}
步骤2.2:根据效用理论,得到第i项评估指标的采集数据zi的置信分布S(zi):
S(zi)={(Fj,βj,i),j=1,…N}
其中,βj,i表示评估指标i的数值对等级Fj的置信度,满足0≤βj,i
Figure BSA0000184791190000036
步骤2.3:确定定值U(Fj)(j=1,2...N)作为评估等级Fj的评估效用值,假设等级Fj+1效用值高于Fj的效用值,可以求得评估指标i的数值对等级Fj的置信度为:
Figure BSA0000184791190000037
βj+1,i(zi)=1-βj,i(zi),U(Fj)≤zi≤U(Fj+1)
βp,i(zi)=0,p=1,...,N,p≠j,j+1
其中,j=1,2,...,N-1。
进一步地,所述步骤3中的基于变异系数法求取评估指标融合权重的过程包括如下步骤:
步骤3.1:基于变异系数法,求得指标xi在k时刻的融合权重ω1i(k):
Figure BSA0000184791190000041
其中,
Figure BSA0000184791190000042
表示评估指标i的变异系数。
步骤3.2:评估指标i的变异系数可由下式求取:
Figure BSA0000184791190000043
其中,
Figure BSA0000184791190000044
表示评估指标xi在预测时间序列T内的样本均方差值,
Figure BSA0000184791190000045
表示样本均值。
步骤3.3:样本均方差值和样本均值可分别由下式求取:
Figure BSA0000184791190000046
Figure BSA0000184791190000047
进一步地,所述步骤4中的基于信息熵法求取评估指标融合权重的过程包括如下步骤:
步骤4.1:根据信息熵特点,定义评估指标i在时刻k的特征比重pi(k):
Figure BSA0000184791190000048
步骤4.2:根据信息熵理论,求得评估指标i的信息熵值ei
Figure BSA0000184791190000049
其中,b为正常数。
步骤4.3:定义评估指标i的差异性系数为1-ei,对其进行归一化处理,得到基于信息熵法求得的评估指标i在时刻k的归一化融合权重ω2i(k):
Figure BSA00001847911900000410
其中,f=1,2,...,n,n表示评估指标数量。
进一步地,所述步骤5中的基于卡方距离实现评估指标组合赋权的过程包括如下步骤:
步骤5.1:定义评估指标融合的组合权重ω=(ω1,ω2,...,ωn),基于变异系数法和信息熵法得到的单一融合权重分别为ω1=(ω11,ω12,...,ω1n)和ω2=(ω21,ω22,...,ω2n)。
步骤5.2:根据卡方距离思想,定于组合赋权的加权系数分别为α1和α2,得到组合权重ω:
ω=α1ω12ω2
步骤5.3:基于卡方距离,对组合赋权方案进行优化,得到优化模型如下:
Figure BSA0000184791190000051
步骤5.4:求解步骤5.3中的优化问题,得到评估指标融合的组合权重ω。
进一步地,所述步骤6中的基于评估指标融合模型求取k时刻深基坑安全性评估结果的过程包括如下步骤:
步骤6.1:基于D-S证据推理理论,融合评估指标信息,得到k时刻深基坑安全性评估结果置信分布O(g(k)):
O(g(k))={(Fj,βj(k)),j=1,...,N}
其中,输出转换模型O将深基坑安全性评估结果g(k)转换为置信分布形式,βj(k)表示k时刻深基坑安全性对于评估等级Fj的置信度。
步骤6.2:根据ER解析算法,求解得到k时刻深基坑安全性对于评估等级Fj的置信度βj(k)。
Figure BSA0000184791190000052
其中,ωi(k)表示评估指标xi在k时刻的融合权重,由步骤6获得;n表示评估指标数量。
进一步地,所述步骤7中的基于多时刻评估结果融合模型求取k时刻深基坑安全性综合评估结果的过程包括如下步骤:
步骤7.1:基于D-S证据推理理论,融合多时刻评估结果,得到k时刻深基坑安全性综合评估结果置信分布O(G(k)):
O(G(k))={(Fj,β′j(k)),j=1,...,N}
其中,βj′(k)表示k时刻下深基坑安全性对评估等级Fj的置信度。
步骤7.2:根据ER解析算法,得到k时刻下深基坑安全性对评估等级Fj的置信度β′j(k):
Figure BSA0000184791190000053
其中,ω′(k-1),ω′(k),ω′(k+1)分别为历史、当前、未来预测时刻下,根据步骤6求得的深基坑安全性评估结果的融合权重,三者满足如下条件:
Figure BSA0000184791190000061
进一步地,所述步骤8中的基于期望效用理论求取深基坑安全性评估数值的过程包括如下步骤:
步骤8.1:引入期望效用,将深基坑安全性评估置信分布变换为最终的安全性综合评估数值,得到评估数值Q(k):
Figure BSA0000184791190000062
其中,Uj表示等级Fj的效用。
步骤8.2:根据深基坑安全评估需求,使安全性综合评估数值对应表征N个评估等级。
有益效果:本发明与现有技术相比,具有以下有点:
在现有的深基坑项目安全评估方法中,大多根据专家经验或现场调查方法实现定性分析,结合部分定量分析方法,做出粗略的安全评估分级,或针对部分重要指标设置一定的安全报警阈值,避免发生严重事故。安全评估方法存在一定的主观性、不确定性和滞后性,难以实现及时有效的深基坑安全评估,无法为深基坑项目的智能管理和长期风险方法提供支持。本发明中基于D-S证据推理对深基坑项目进行安全性评估,依托项目检测数据,无需建立复杂的评估模型,且不需要大量的项目风险数据进行模型训练,具备较高的方案可实施性。通过对深基坑未来状态信息进行预测,有效提高了评估模型的预测性能,对未来可能出现的风险信息进行提前响应。基于变异系数法和信息熵法对评估指标融合权重进行求取,实质上是考虑评估指标变化幅度和包含有效信息量多少对深基坑项目安全性评估的表征能力,有效提高评估的准确性和全面性。
本发明在深基坑安全性评估过程中,利用历史和当前检测数据,基于三阶Volterra级数对深基坑未来状态进行预测,能够提高方案对项目潜在风险的预测能力。同时,利用有限长单位冲激响应滤波器对预测模型进行简化,减少模型系数,提高了项目预测速度。采用变异系数和信息熵对评估指标融合权重进行确定,有效考虑了指标变化幅度和包含信息量对安全评估的表征作用,提高评估的准确性。利用效用理论和D-S证据推理进行安全性评估,有效地实现项目安全性的定性和定量分析。与其他深基坑项目安全性评估方法相比,本发明提供的方法具备更好的方案可实现性和评估准确性,无需太多专家人员的介入,方便后期方案的有效使用和推广。
附图说明
图1为本发明的深基坑安全性评估整体框架。
图2为本发明的深基坑安全性评估流程图。
具体实施方式
下面结合实施例和说明书对本发明作进一步的说明。
本发明的一种基于效用理论和D-S证据推理的深基坑安全性评估方法,其特征在于,该方法包括以下步骤:
步骤1:使用基于三阶Volterra级数建立的未来状态非线性预测模型,通过历史和当前数据对深基坑未来状态进行预测。
其中,深基坑未来状态预测方法包括如下步骤:
步骤1.1:以深基坑安全性指标的历史数据和当前观测值作为模型输入,如下式所示:
Z(k)=[z(k),z(k-1),…,z(k-T+1)]T
其中,z(k)表示评估指标在k时刻的采集数值,T表示用于预测未来状态的时间序列长度。
步骤1.2:基于三阶Volterra级数建立未来状态预测模型,得到深基坑未来状态预测值y(k):
Figure BSA0000184791190000071
步骤1.3:引入有限长单位冲激响应滤波器,并利用乘积耦合减少预测模型系数,得到简化系数后的模型输出预测值y(k):
Figure BSA0000184791190000072
其中,hm,r表示滤波器系数。
步骤1.4:利用改进的非线性归一化LMS滤波算法(NNLMS),达到提高滤波器收敛速度的目的,得到近似的滤波器系数为:
Figure BSA0000184791190000073
Figure BSA0000184791190000074
Figure BSA0000184791190000075
其中,A(k),B(k),C(k)分别表示三个有限长单位冲激响应滤波器系数,Z′(k)=[1,z(k)...z(k-T+1)]T
Figure BSA0000184791190000076
A(k)=[a0,a1…aT]T,B(k)=[b0,b1…bT]T,C(k)=[c0,c1…cT]T,并且μ1≠μ2≠μ3
步骤2:基于效用理论,对评估指标数据进行处理,得到评估指标关于各评估等级的置信分布
其中,评估指标关于评估等级置信分布求取过程包括如下步骤:
步骤2.1:针对深基坑安全性评估,定义N个安全性等级,得到深基坑评估等级的完备集F:
F={F1,…,Fj,…,FN}
步骤2.2:根据效用理论,得到第i项评估指标的采集数据zi的置信分布S(zi):
S(zi)={(Fj,βj,i),j=1,…N}
其中,βj,i表示评估指标i的数值对等级Fj的置信度,满足0≤βj,i
Figure BSA0000184791190000077
步骤2.3:确定定值U(Fj)(j=1,2...N)作为评估等级Fj的评估效用值,假设等级Fj+1效用值高于Fj的效用值,可以求得评估指标i的数值对等级Fj的置信度为:
Figure BSA00001847911900000810
βj+1,i(zi)=1-βj,i(zi),U(Fj)≤zi≤U(Fj+1)
βp,i(zi)=0,p=1,...,N,p≠j,j+1
其中,j=1,2,...,N-1。
步骤3:基于变异系数法确定评估指标融合权重,得到评估指标xi在k时刻的变异系数融合权重。
其中,基于变异系数法求取评估指标融合权重的过程包括如下步骤:
步骤3.1:基于变异系数法,求得指标xi在k时刻的融合权重ω1i(k):
Figure BSA0000184791190000082
其中,
Figure BSA0000184791190000089
表示评估指标i的变异系数。
步骤3.2:评估指标i的变异系数可由下式求取:
Figure BSA0000184791190000083
其中,
Figure BSA0000184791190000084
表示评估指标xi在预测时间序列T内的样本均方差值,
Figure BSA0000184791190000085
表示样本均值。
步骤3.3:样本均方差值和样本均值可分别由下式求取:
Figure BSA0000184791190000086
Figure BSA0000184791190000087
步骤4:基于信息熵法确定评估指标融合权重,得到评估指标xi在k时刻的信息熵融合权重。
其中,基于信息熵法求取评估指标融合权重的过程包括如下步骤:
步骤4.1:根据信息熵特点,定义评估指标i在时刻k的特征比重pi(k):
Figure BSA0000184791190000088
步骤4.2:根据信息熵理论,求得评估指标i的信息熵值ei
Figure BSA0000184791190000091
其中,b为正常数。
步骤4.3:定义评估指标i的差异性系数为1-ei,对其进行归一化处理,得到基于信息熵法求得的评估指标i在时刻k的归一化融合权重ω2i(k):
Figure BSA0000184791190000092
其中,f=1,2,...,n,n表示评估指标数量。
步骤5:基于卡方距离实现评估指标组合赋权,得到评估指标最终融合权重向量。
其中,基于卡方距离实现评估指标组合赋权的过程包括如下步骤:
步骤5.1:定义评估指标融合的组合权重ω=(ω1,ω2,...,ωn),基于变异系数法和信息熵法得到的单一融合权重分别为ω1=(ω11,ω12,...,ω1n)和ω2=(ω21,ω22,...,ω2n)。
步骤5.2:根据卡方距离思想,定于组合赋权的加权系数分别为α1和α2,得到组合权重ω:
ω=α1ω12ω2
步骤5.3:基于卡方距离,对组合赋权方案进行优化,得到优化模型如下:
Figure BSA0000184791190000093
步骤5.4:求解步骤5.3中的优化问题,得到评估指标融合的组合权重ω。
步骤6:基于D-S证据推理建立评估指标融合模型,得到k时刻深基坑安全性评估结果的置信分布。
其中,基于评估指标融合模型求取k时刻深基坑安全性评估结果的过程包括如下步骤:
步骤6.1:基于D-S证据推理理论,融合评估指标信息,得到k时刻深基坑安全性评估结果置信分布O(g(k)):
O(g(k))={(Fj,βj(k)),j=1,...,N}
其中,输出转换模型O将深基坑安全性评估结果g(k)转换为置信分布形式,βj(k)表示k时刻深基坑安全性对于评估等级Fj的置信度。
步骤6.2:根据ER解析算法,求解得到k时刻深基坑安全性对于评估等级Fj的置信度βj(k)。
Figure BSA0000184791190000094
其中,ωi(k)表示评估指标xi在k时刻的融合权重,由步骤6获得;n表示评估指标数量。
步骤7:基于D-S证据推理建立多时刻深基坑安全性评估结果融合模型,融合历史时刻、当前时刻和未来预测时刻的安全性评估结果,得到k时刻深基坑综合评估结果的置信分布。
其中,基于多时刻评估结果融合模型求取k时刻深基坑安全性综合评估结果的过程包括如下步骤:
步骤7.1:基于D-S证据推理理论,融合多时刻评估结果,得到k时刻深基坑安全性综合评估结果置信分布O(G(k)):
O(G(k))={(Fj,β′j(k)),j=1,...,N}
其中,βj′(k)表示k时刻下深基坑安全性对评估等级Fj的置信度。
步骤7.2:根据ER解析算法,得到k时刻下深基坑安全性对评估等级Fj的置信度β′j(k):
Figure BSA0000184791190000101
其中,ω′(k-1),ω′(k),ω′(k+1)分别为历史、当前、未来预测时刻下,根据步骤6求得的深基坑安全性评估结果的融合权重,三者满足如下条件:
Figure BSA0000184791190000102
步骤8:基于期望效用理论,将深基坑安全性评估置信分布变换为安全性评估数值,进而得到评估等级。
其中,基于期望效用理论求取深基坑安全性评估数值的过程包括如下步骤:
步骤8.1:引入期望效用,将深基坑安全性评估置信分布变换为最终的安全性综合评估数值,得到评估数值Q(k):
Figure BSA0000184791190000103
其中,Uj表示等级Fj的效用。
步骤8.2:根据深基坑安全评估需求,使安全性综合评估数值对应表征N个评估等级。
上述实施例仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和等同替换,这些对本发明权利要求进行改进和等同替换后的技术方案,均落入本发明的保护范围。

Claims (9)

1.一种基于效用理论和D-S证据推理的深基坑安全性评估方法,其特征在于,该方法包括以下步骤:通过深基坑现场传感器设备获取用于安全性评估的评估指标数据后;
步骤1:使用基于三阶Volterra级数建立的未来状态非线性预测模型,通过历史和当前数据对深基坑未来状态进行预测,最终得到状态预测值y(k):
Figure FDA0003989379900000011
其中,k表示时刻k,hm,r表示滤波器系数;
步骤2:基于效用理论,对评估指标数据进行处理,得到评估指标关于各评估等级的置信分布S(zi):
S(zi)={(Fjj,i),j=1,···N}
其中,zi表示评估指标i的数值,Fj表示第j级评估等级,共N个评估等级;
步骤3:基于变异系数法确定评估指标融合权重,得到评估指标i在k时刻的变异系数融合权重ω1i(k):
Figure FDA0003989379900000012
其中,
Figure FDA0003989379900000013
表示评估指标i的变异系数;
步骤4:基于信息熵法确定评估指标融合权重,得到评估指标i在k时刻的信息熵融合权重ω2i(k):
Figure FDA0003989379900000014
其中,ei表示评估指标i的信息熵值,f=1,2,…,n,n表示评估指标数量
Figure FDA0003989379900000015
步骤5:基于卡方距离实现评估指标组合赋权,得到评估指标最终融合权重向量ω:
ω=α1ω12ω2
其中,α1和α2表示组合赋权的加权系数;
步骤6:基于D-S证据推理建立评估指标融合模型,得到k时刻深基坑安全性评估结果的置信分布:
O(g(k))={(Fjj(k)),j=1,...,N}
其中,输出转换模型O(g(k))将深基坑安全性评估结果g(k)转换为置信分布形式,βj(k)表示k时刻深基坑安全性对于评估等级Fj的置信度;
步骤7:基于D-S证据推理建立多时刻深基坑安全性评估结果融合模型,融合历史时刻、当前时刻和未来预测时刻的安全性评估结果,得到k时刻深基坑综合评估结果的置信分布:
O(G(k))={(Fjj'(k)),j=1,...,N}
其中,βj'(k)表示k时刻下系统安全性对于等级Fj的置信度;
步骤8:基于期望效用理论,将深基坑安全性评估置信分布变换为安全性评估数值,进而得到评估等级。
2.根据权利要求1所述的基于效用理论和D-S证据推理的深基坑安全性评估方法,其特征在于,所述步骤1中的深基坑未来状态预测方法包括如下步骤:
步骤1.1:以深基坑安全性指标的历史数据和当前观测值作为模型输入,如下式所示:
Z(k)=[z(k),z(k-1),···,z(k-T+1)]T
其中,z(k)表示评估指标在k时刻的采集数值,T表示用于预测未来状态的时间序列长度;
步骤1.2:基于三阶Volterra级数建立未来状态预测模型,得到深基坑未来状态预测值y(k):
Figure FDA0003989379900000021
其中,k表示时刻k,
Figure FDA0003989379900000027
表示评估指标在k+1时刻的状态预测值;
步骤1.3:引入有限长单位冲激响应滤波器,并利用乘积耦合减少预测模型系数,得到简化系数后的模型输出预测值y(k):
Figure FDA0003989379900000022
其中,hm,r表示滤波器系数;
步骤1.4:利用改进的非线性归一化LMS滤波算法(NNLMS),达到提高滤波器收敛速度的目的,得到近似的滤波器系数为:
Figure FDA0003989379900000023
Figure FDA0003989379900000024
Figure FDA0003989379900000025
其中,A(k),B(k),C(k)分别表示三个有限长单位冲激响应滤波器系数,Z′(k)=[1,z(k)…z(k-T+1)]T
Figure FDA0003989379900000026
A(k)=[a0,a1···aT]T,B(k)=[b0,b1···bT]T,C(k)=[c0,c1···cT]T,并且μ1≠μ2≠μ3
3.根据权利要求1所述的基于效用理论和D-S证据推理的深基坑安全性评估方法,其特征在于,所述步骤2中的评估指标关于评估等级置信分布求取过程包括如下步骤:
步骤2.1:针对深基坑安全性评估,定义N个安全性等级,得到深基坑评估等级的完备集F:
F={F1,···,Fj,···,FN}
步骤2.2:根据效用理论,得到第i项评估指标的采集数据zi的置信分布S(zi):
S(zi)={(Fjj,i),j=1,···N}
其中,βj,i表示评估指标i的数值对等级Fj的置信度,满足0≤βj,i
Figure FDA0003989379900000031
步骤2.3:确定定值U(Fj)(j=1,2…N)作为评估等级Fj的评估效用值,假设等级Fj+1效用值高于Fj的效用值,可以求得评估指标i的数值对等级Fj的置信度为:
Figure FDA0003989379900000032
βj+1,i(zi)=1-βj,i(zi),U(Fj)≤zi≤U(Fj+1)
βp,i(zi)=0,p=1,...,N,p≠j,j+1
其中,j=1,2,…,N-1。
4.根据权利要求1所述的基于效用理论和D-S证据推理的深基坑安全性评估方法,其特征在于,所述步骤3中的基于变异系数法求取评估指标融合权重的过程包括如下步骤:
步骤3.1:基于变异系数法,求得评估指标i在k时刻的融合权重ω1i(k):
Figure FDA0003989379900000033
其中,
Figure FDA00039893799000000310
表示评估指标i的变异系数;
步骤3.2:评估指标i的变异系数可由下式求取:
Figure FDA0003989379900000034
其中,
Figure FDA0003989379900000035
表示评估指标i在预测时间序列T内的样本均方差值,
Figure FDA0003989379900000036
表示样本均值;
步骤3.3:样本均方差值和样本均值可分别由下式求取:
Figure FDA0003989379900000037
Figure FDA0003989379900000038
5.根据权利要求1所述的基于效用理论和D-S证据推理的深基坑安全性评估方法,其特征在于,所述步骤4中的基于信息熵法求取评估指标融合权重的过程包括如下步骤:
步骤4.1:根据信息熵特点,定义评估指标i在时刻k的特征比重pi(k):
Figure FDA0003989379900000039
步骤4.2:根据信息熵理论,求得评估指标i的信息熵值ei
Figure FDA0003989379900000041
其中,b为正常数;
步骤4.3:定义评估指标i的差异性系数为1-ei,对其进行归一化处理,得到基于信息熵法求得的评估指标i在时刻k的归一化融合权重ω2i(k):
Figure FDA0003989379900000042
其中,f=1,2,…,n,n表示评估指标数量。
6.根据权利要求1所述的基于效用理论和D-S证据推理的深基坑安全性评估方法,其特征在于,所述步骤5中的基于卡方距离实现评估指标组合赋权的过程包括如下步骤:
步骤5.1:定义评估指标融合的组合权重ω=(ω12,…,ωn),基于变异系数法和信息熵法得到的单一融合权重分别为ω1=(ω1112,…,ω1n)和ω2=(ω2122,…,ω2n);
步骤5.2:根据卡方距离思想,定于组合赋权的加权系数分别为α1和α2,得到组合权重ω:
ω=α1ω12ω2
步骤5.3:基于卡方距离,对组合赋权方案进行优化,得到优化模型如下:
Figure FDA0003989379900000043
步骤5.4:求解步骤5.3中的优化问题,得到评估指标融合的组合权重ω。
7.根据权利要求1所述的基于效用理论和D-S证据推理的深基坑安全性评估方法,其特征在于,所述步骤6中的基于评估指标融合模型求取k时刻深基坑安全性评估结果的过程包括如下步骤:
步骤6.1:基于D-S证据推理理论,融合评估指标信息,得到k时刻深基坑安全性评估结果置信分布O(g(k)):
O(g(k))={(Fjj(k)),j=1,...,N}
其中,输出转换模型O将深基坑安全性评估结果g(k)转换为置信分布形式,βj(k)表示k时刻深基坑安全性对于评估等级Fj的置信度;
步骤6.2:根据ER解析算法,求解得到k时刻深基坑安全性对于评估等级Fj的置信度βj(k);
Figure FDA0003989379900000051
其中,ωi(k)表示评估指标xi在k时刻的融合权重,由步骤6获得;n表示评估指标数量。
8.根据权利要求1所述的基于效用理论和D-S证据推理的深基坑安全性评估方法,其特征在于,所述步骤7中的基于多时刻评估结果融合模型求取k时刻深基坑安全性综合评估结果的过程包括如下步骤:
步骤7.1:基于D-S证据推理理论,融合多时刻评估结果,得到k时刻深基坑安全性综合评估结果置信分布O(G(k)):
O(G(k))={(Fjj'(k)),j=1,...,N}
其中,βj'(k)表示k时刻下深基坑安全性对评估等级Fj的置信度;
步骤7.2:根据ER解析算法,得到k时刻下深基坑安全性对评估等级Fj的置信度βj'(k):
Figure FDA0003989379900000052
其中,ω'(k-1),ω'(k),ω'(k+1)分别为历史、当前、未来预测时刻下,根据步骤6求得的深基坑安全性评估结果的融合权重,三者满足如下条件:
Figure FDA0003989379900000053
9.根据权利要求1所述的基于效用理论和D-S证据推理的深基坑安全性评估方法,其特征在于,所述步骤8中的基于期望效用理论求取深基坑安全性评估数值的过程包括如下步骤:
步骤8.1:引入期望效用,将深基坑安全性评估置信分布变换为最终的安全性综合评估数值,得到评估数值Q(k):
Figure FDA0003989379900000054
其中,Uj表示等级Fj的效用;
步骤8.2:根据深基坑安全评估需求,使安全性综合评估数值对应表征N个评估等级。
CN201910534172.XA 2019-06-19 2019-06-19 一种基于效用理论和d-s证据推理的深基坑安全性评估方法 Active CN110288144B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910534172.XA CN110288144B (zh) 2019-06-19 2019-06-19 一种基于效用理论和d-s证据推理的深基坑安全性评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910534172.XA CN110288144B (zh) 2019-06-19 2019-06-19 一种基于效用理论和d-s证据推理的深基坑安全性评估方法

Publications (2)

Publication Number Publication Date
CN110288144A CN110288144A (zh) 2019-09-27
CN110288144B true CN110288144B (zh) 2023-04-11

Family

ID=68004021

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910534172.XA Active CN110288144B (zh) 2019-06-19 2019-06-19 一种基于效用理论和d-s证据推理的深基坑安全性评估方法

Country Status (1)

Country Link
CN (1) CN110288144B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111504963B (zh) * 2020-04-10 2023-07-07 上海蓝长自动化科技有限公司 一种应用于叶绿素、蓝绿藻荧光检测的数据时空融合方法
CN111861094A (zh) * 2020-05-28 2020-10-30 南京邮电大学 一种分布式电源接入配电网的安全性态势综合评价方法
CN111967791B (zh) * 2020-08-28 2021-08-31 中国人民解放军火箭军工程大学 一种考虑扰动与故障阈值的设备性能评估方法及系统
CN112101757A (zh) * 2020-09-03 2020-12-18 湖北文理学院 一种深基坑施工安全风险评估方法及系统
CN112288308A (zh) * 2020-11-11 2021-01-29 华能澜沧江水电股份有限公司 基于多重变异系数赋权的ds证据理论滑坡评价方法
CN115638808B (zh) * 2022-12-23 2023-03-21 中国人民解放军火箭军工程大学 综合时空域信息的惯导系统质量状态确定方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106156343A (zh) * 2016-07-21 2016-11-23 中国地质大学(武汉) 一种深基坑施工方案安全评审知识库以及自动安全评审方法
CN206775646U (zh) * 2017-05-12 2017-12-19 华中科技大学 一种工程施工安全监测机构

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104408324B (zh) * 2014-12-11 2017-06-13 云南师范大学 基于d‑s证据理论的多传感器信息融合方法
CN104899682A (zh) * 2015-05-19 2015-09-09 上海市建工设计研究院有限公司 临海地区深基坑防渗止水帷幕体施工风险评价方法
CN105373963A (zh) * 2015-09-14 2016-03-02 中国电力科学研究院 一种基于组合权重electre评价模型的发电计划评估方法
CN108537448B (zh) * 2018-04-12 2022-04-29 西南交通大学 基于集对分析和证据理论的高铁接触网健康状态评估方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106156343A (zh) * 2016-07-21 2016-11-23 中国地质大学(武汉) 一种深基坑施工方案安全评审知识库以及自动安全评审方法
CN206775646U (zh) * 2017-05-12 2017-12-19 华中科技大学 一种工程施工安全监测机构

Also Published As

Publication number Publication date
CN110288144A (zh) 2019-09-27

Similar Documents

Publication Publication Date Title
CN110288144B (zh) 一种基于效用理论和d-s证据推理的深基坑安全性评估方法
Tian et al. Multicriteria decision‐making approach based on gray linguistic weighted Bonferroni mean operator
CN108230049A (zh) 订单的预测方法及系统
CN109214581B (zh) 一种考虑风向和置信区间的铁路沿线风速预测方法
CN115169479A (zh) 污水处理过程远程监控方法、系统及存储介质
CN113156917A (zh) 基于人工智能的电网设备故障诊断方法及系统
CN110110425B (zh) 一种基于信度规则推理的边坡滑动力预测方法
CN105825271A (zh) 基于证据推理的卫星故障诊断与预测方法
CN112149976B (zh) 一种基于人工智能的电网精准投资项目决策方法
Hu et al. A particle filter and long short-term memory fusion technique for lithium-ion battery remaining useful life prediction
CN108961460B (zh) 基于稀疏esgp与多目标优化的故障预测方法及装置
CN115588161A (zh) 一种基于数字孪生技术的数字化系统及其构建方法
CN115330203A (zh) 一种基于博弈证据网络的液化天然气储罐健康状态评估方法
CN116933946A (zh) 一种基于客流去向结构的轨道交通od客流预测方法及系统
Zhao et al. Monte Carlo method and quantile regression for uncertainty analysis of wind power forecasting based on Chaos-LS-SVM
KR20200091508A (ko) 특허와 논문 데이터를 활용한 국가 및 기업들의 과학 기술력 진단 및 예측 방법
CN102095953A (zh) 一种蓄电池充电机性能在线检测方法
CN112541636A (zh) 输电线路覆冰风险预警方法、装置、介质和电子设备
CN114282658B (zh) 一种针对流量序列分析及预测的方法、设备和介质
CN116451322A (zh) 基于贝叶斯优化的lstm深度学习网络的力学预测方法
CN114692692B (zh) 一种基于微波衰减信号融合核极限学习机的降雪识别方法
CN116070368A (zh) 一种海上风电机组高速轴承剩余寿命预测方法
CN110986873B (zh) 高速铁路工程服役状态预警指标获取方法
CN110556595B (zh) 储能用磷酸铁锂电池安全管理系统的质态失衡预测方法
CN114065624A (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
GR01 Patent grant
GR01 Patent grant