一种基于蒙特卡罗法的地源热泵竖直地埋管的可靠性分析方法
【技术领域】
[0001] 本发明属于地源热栗系统的可靠性分析技术领域,涉及一种基于蒙特卡罗法的地 源热栗竖直地埋管的可靠性分析方法。
【背景技术】
[0002] 在地源热栗竖直地埋管的研究中,传热分析一直是许多学者研究的重点,也因此 产生了许多不同的传热模型。基于传热分析与传热模型的建立,地源热栗竖直地埋管的设 计也出现了相应的计算公式和方法。现有的地埋管计算公式和设计方法中,极少部分有考 虑相关参数的动态变化和随机性。Jeffrey D. Spitler等人提到地下导热率的不确定性和 分布范围,王勇做了建筑负荷的动态变化对地源热栗性能影响的相关研究,管昌生等人考 虑了岩土热阻的随机性,并做了相关研究,这些在进行地源热栗地埋管管长计算或可靠性 分析时,都没有全面考虑影响地源热栗的随机因素,只是在假定其他因素不变的前提下,某 一特定因素变化所导致的评价指标(本发明采用竖直地埋管管长这一评价指标)的变化程 度,无法反映各随机因素同时变化引起评价指标的变化程度。
【发明内容】
[0003] 本发明的目的在于克服上述技术存在的缺陷,提供一种基于蒙特卡罗法的地源热 栗竖直地埋管的可靠性分析方法,在计算机环境下,借助EXCEL软件提供的函数功能,拟以 竖直地埋管管长计算为例,应用蒙特卡罗(Monte Carlo)模拟法对随机因素影响下的竖直 地埋管管长进行成千上万次的模拟运算,最后利用SPSS软件对所得到的样本进行处理,对 地埋管的可靠性进行分析。
[0004] 其具体技术方案为:
[0005] -种基于蒙特卡罗法的地源热栗竖直地埋管的可靠性分析方法,包括以下步骤:
[0006] (1)分析评价参数的特性,确定随机变量的某些概率分布律;
[0007] (2)按照一定的统计分布规律,在计算机上产生随机数,在此基础上,建立相应的 数学模型;
[0008] (3)通过成千上万次计算机模拟运算,得到足够数量的参数样本值,根据这些参数 样本值,进行统计特征分析,求出地源热栗竖直地埋管管长的概率分布;
[0009] (4)根据上述所求的概率分布,分析评价预测管长的可靠性。
[0010] 优选地,在进行模拟运算时,主要采用的函数包括:RAND:是一个Excel随机数字发 生器,产生〇到1之间的一个随机数;N0RMINV:是一个返回给定概率、正态分布的均值和标准 差的区间点;AVERAGE:计算指标的平均值;STDEV:估算基于样本的标准偏差;MAX:计算样本 的最大值;MIN:计算样本的最小值。
[0011] 与现有技术相比,本发明的有益效果为:
[0012] (1)基于EXCEL环境下,采用蒙特卡罗(Monte Carlo)模拟法,可以解决现行可靠性 计算中仅对特定因素变化引起评价指标变化的可靠性分析问题,解决了各随机因素同时变 化引起评价指标的变化的可靠性分析问题。
[0013] (2)本发明得到的固定样本是成千上万次实验模拟的结果,具有一定的物理意义 和数学意义,是可信的。
[0014] (3)根据系统不同的保证率,可以得到不同的计算管长。就本发明的实例而言,与 现有技术相比,根本无法保证其置信度达到95.44% ( "2〇"范围内);若按照100m的井深计 算,要想保证地源热栗系统的可靠性在"2〇"范围内,则需要额外至少增加1 口井,不超过5口 井(否则会增加初投资);要想保证地源热栗系统的可靠性在V'范围内,则需要额外至少增 加2 口井;不超过4 口井(否则会增加初投资)。
[0015] (4)本发明实例所涉及随机参数及其分布类型,是在假设的基础上提出来的,为工 程实际提供一定的参考。
【附图说明】
[0016] 图1是管长计算基本资料;
[0017] 图2是参数随机取数后一次计算的结果;
[0018] 图3是模拟运算5000次的部分数据;
[0019]图4是固定样本的统计参数计算;
[0020] 图5是5000个样本频率分布直方图。
【具体实施方式】
[0021]为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结 合附图和具体实例,进一步阐述本发明。
[0022] 1蒙特卡罗模拟法原理及步骤
[0023] 1.1蒙特卡罗模拟法的基本原理
[0024]蒙特卡罗(Monte Carlo)方法,是一种基于随机数的统计模拟方法。在一些统计问 题中,有时候很难给出完美的理论结果,这时可以通过随机模拟来给出近似的结果。一般需 要根据已知的概率分布构造出相应的模拟样本,用它们的样本频率代替相应的概率作统计 分析和推断。
[0025]蒙特卡罗分析是情景分析法中的一种,在该分析方法下,被分析的项目将面临大 量的情景。以竖直地埋管管长为评价指标,在每一种情景下,计算机为每个变量随机选择一 个数值,然后运用这些数据计算地埋管管长,并将这些管长值存储到计算机中。接下来,再 为输入变量随机选取第二组数据,并计算第二个管长值。通过成千上万次的模拟,可以获得 相应的概率分布,通过对大量实验样本进行统计规律分析,得到满足一定可信度的结果。通 过蒙特卡罗模拟,可以对地源热栗竖直地埋管管长计算中的多个不确定性因素进行有效分 析,克服常规可靠性分析中仅对单一随机因素分析的弊端。
[0026] 1.2蒙特卡罗模拟法的分析步骤
[0027]在工程技术问题中,常常很难获取实验数据,而且会耗费很多人力、物力。对此采 用计算机随机模拟的方法无疑显得简单、实用而经济。以地源热栗竖直地埋管管长计算为 例,蒙特卡罗模拟法的分析步骤如下:
[0028] (1)分析评价参数的特性,并根据历史资料或专家意见,确定随机变量的某些概率 分布律;
[0029] (2)按照一定的统计分布规律,在计算机上产生随机数,在此基础上,建立相应的 数学模型;
[0030] (3)通过成千上万次计算机模拟运算,得到足够数量的参数样本值,根据这些参数 样本值,进行统计特征分析,求出地源热栗竖直地埋管管长的概率分布;
[0031] (4)根据上述所求的概率分布,分析评价预测管长的可靠性。
[0032] 1.3基于Excel的蒙特卡罗模拟所需函数
[0033] Excel软件内含有大量的统计函数,可以实时反映各种变量之间的变化,获得各随 机因素下目标函数的变动规律。在进行蒙特卡罗模拟运算时,主要采用的函数包括:RAND: 是一个Excel随机数字发生器,产生0到1之间的一个随机数;N0RMINV:是一个返回给定概 率、正态分布的均值和标准差的区间点;AVERAGE:计算指标的平均值;STDEV:估算基于样本 的标准偏差;MAX:计算样本的最大值;MIN:计算样本的最小值。
[0034] 2管长计算公式
[0035]由于地埋管换热器中的传热过程是三维非稳态的传热,影响因素非常复杂,很难 用简单的公式加以描述和概括,因此在实际工程中广泛采用以半经验公式为主的设计计算 方法,主要根据最大冷、热负荷估算地埋管所需埋管的长度。
[0036]地埋管换热器钻孔总长度:
[0037]
[0038]
[0039]
[0040]
[0041 ]其中,下标Η代表供热,C代表制冷;L是钻孔长度;Q是热栗在设计进水温度下的额 定出力(供热或制冷量);C0P是热栗的性能系数;F是运行份额;RP和Rs分别是管壁的热阻和 岩土的热阻;t〇〇是地下未受干扰时的平均温度;别为最高和最低进水温度;Τη是 最热月中的运行小时数;Tc是最冷月中的运行小时数;d H是最热月的天数;dc是最冷月的天 数。
[0042]由公式(1)、(2)可知,Q、RP、Rs、F、C0P、t等因素都具有随机性,于是L也是随机变化 的。为同时满足供热和供冷的需求,应采用Lh和U中的较大者作为设计钻孔总长度。
[0043] 3管长计算公式
[0044] K-S检验法称为柯尔莫哥洛夫一斯米尔洛夫检验法(Kolmogorov-Smirnove Test),这个检验通过比较样本经验分布函数?"(1)与总体分布函数F(x)的差异构造检验统 计量,克服了 X2检验依赖于区间划分的缺点。
[0045] K-S统计量为
[0046]
(5)
[0047] 其中,经骀分布函数Fn(x):
[0048]
[0049] 显然,拒绝域的形式为Dn 2 k,临界值k由检验水平α和样本容量η确定,记为Dn,a,满 足P(Dn2D n,a)=a,其中,Dn,a的值可由现有技术查得。
[0050] 4实例分析
[0051 ] 4.1基本资料
[0052]武汉市某花园小区总建筑面积38000m2,由于整个建筑采取了节能措施,经过计 算,最大冷负荷为1560kW,最大热负荷为1000kW。地源热栗系统拟采用DN32单U型管,管材为 HDPE。假设各随机变量服从正态分布,且各随机变量的概率分布资料见表1所不。
[0053]表1随机
变量的统计参数及分布类型
[0054]
[0055] 4.2变量表现形式
[0056]由于上述随机变量均服从正态分布,故采用N0RMINV、RAND函数来产生随机数,相 关参数的表现形式见表2所示。
[0057] 表2各随机变量的表现形式
[0058]
[0059] ?.ο vi
[0060] (1)在Excel环境下,建立如图1所示的工作表,输入随机参数的基本资料,并计算 相应的管长Lc。
[0061 ] (2)在C2~C8中分别输入如表2所示的符合问题要求的随机变量。输入完上述变量 的计算公式后,在单元格C9(如图2所示)得到上述随机变量一次实验后的管长值。
[0062] (3)将上述模拟计算5000次,并且将5000次的计算结果存储在Excel工作表中,如 图3所示(由于样本数据太多,只截取部分数据)。
[0063] (4)上述5000次模拟运算的数据是随机产生的,而且这5000个管长值是"活"的,只 要工作表中做一次修改,这些数据就会改变一次,相当于又重新做了 5000次实验。因而需要 对这5000个管长数据进行统计分析,寻找其统计规律。将上述工作表中的数据E2:E5001复 制新的工作表中(固定实验样本),分别由相应的函数求管长的统计参数,如图4所示。
[0064] (5)利用SPSS软件,对所得的固定样本进行K-S检验,其检验结果如表3所示,频率 分布直方图如表4所示,根据"3〇"法则,管长落在(μ-3〇#+3〇)(99.74%)内几乎是肯定的 事,计算三种不同置信度下管长的置信区间,如表5所示。
[0065] 表3 固定样本Kolmogorov-Smirnov检验
[0066]
[0067] a.检验分布为正态分布。
[0068] b.根据数据计算得到。
[0069] 表4不同的置信度1 一α下管长的置信区间
[0070]
[0071] 4.4模拟结果分析
[0072] 如图1得计算管长LC = 26350m,可知根本无法保证其置信度达到95.44% ("2σ"范 围内);若按照l〇〇m的井深计算,要想保证地源热栗系统的可靠性在"2〇"范围内,则需要额 外至少增加1 口井,不超过5 口井(否则会增加初投资);要想保证地源热栗系统的可靠性在 V'范围内,则需要额外至少增加2口井;不超过4口井(否则会增加初投资)。
[0073]以上所述,仅为本发明最佳实施方式,任何熟悉本技术领域的技术人员在本发明 披露的技术范围内,可显而易见地得到的技术方案的简单变化或等效替换均落入本发明的 保护范围内。
【主权项】
1. 一种基于蒙特卡罗法的地源热栗竖直地埋管的可靠性分析方法,其特征在于,包括 以下步骤: (1) 分析评价参数的特性,确定随机变量的某些概率分布律; (2) 按照一定的统计分布规律,在计算机上产生随机数,在此基础上,建立相应的数学 丰旲型; (3) 通过成千上万次计算机模拟运算,得到足够数量的参数样本值,根据这些参数样本 值,进行统计特征分析,求出地源热栗竖直地埋管管长的概率分布; (4) 根据上述所求的概率分布,分析评价预测管长的可靠性。2. 根据权利要求1所述的基于蒙特卡罗法的地源热栗竖直地埋管的可靠性分析方法, 其特征在于,在进行模拟运算时,主要采用的函数包括:RAND:是一个Excel随机数字发生 器,产生0到1之间的一个随机数;NORMINV:是一个返回给定概率、正态分布的均值和标准差 的区间点;AVERAGE:计算指标的平均值;STDEV:估算基于样本的标准偏差;MAX:计算样本的 最大值;MIN:计算样本的最小值。
【专利摘要】本发明公开了一种基于蒙特卡罗法的地源热泵竖直地埋管的可靠性分析方法,提出在EXCEL环境下,采用蒙特卡罗模拟法,研究各随机变量同时变化时埋管长度的变化程度,通过成千上万次模拟运算,得到大量的管长样本,再通过SPSS软件对样本进行处理。结合工程实例计算,与确定性设计方法进行对比分析,根据系统所要求不同的保证率,可以得到不同的计算管长,为地埋管的可靠性分析提供一定的参考。
【IPC分类】G06Q10/06
【公开号】CN105488612
【申请号】CN201510824733
【发明人】符永正, 鲁芬豹, 刘腊美, 石磊
【申请人】武汉科技大学
【公开日】2016年4月13日
【申请日】2015年11月24日...