本发明属于冻土领域的比热计算领域，特别是一种考虑相变潜热的冻土等效比热计算方法，可用于将冻土的相变潜热转换为土体的等效比热。
背景技术：
冻土的相变潜热，是指土中的液相水与固相冰相态转换过程中释放或吸收的热量。某温度下所对应的物质比热值即为真比热。工程上常采用某温度内比热的平均值作为物体比热，即为物体的平均比热。单位质量的带相变材料经相变过程达到某温度吸收(释放)的热量与温度区间之比，即为等效比热。真比热、平均比热和等效比热统称为比热。冻土比热的测试方法一般分为两种，即直接测试和间接测试。直接测试是将冻土当成一整体物质，进而测试其比热，常见的直接测试方法有混合量热法。间接测试方法是指依据冻土中土骨架、未冻水、冰体和气体的含量和各相的比热值，按照加权计算方法确定冻土的比热。直接法测试冻土比热的过程中，冻土试样在量热器中融化，导致冻土中土颗粒、冰体、未冻水和气体的升温，同时也导致冻土中的部分冰体融化。由于冻土中土颗粒、冰、未冻水和气的比热容不同，升高至同一温度过程中存在外热源与冻土中各相的热传导，以及冻土中各相间的热传导。h2o(水或冰)是晶体材料，其在相变过程中吸收热量而温度不产生变化，依据温度改变所吸收的总热量无法剔除相变热对测试结果的影响，也无法获取某温度点的冻土比热。因此，直接法仅可直观得到材料的平均比热。由加权计算原理可知，依据某温度下各相的比热和各相含量间接确定的为材料真比热。可见直接法和间接法均无法直接获取冻土的等效比热。然而，在冻土地基、冻土区水合物沉积物的温度场计算中常常需要将潜热转换为等效比热，以满足通用温度场计算软件对热参数的要求，提升温度场计算效率。因而，迫切需要提供一种能够将冻土相变潜热转换为等效比热的计算方法，该方法对于冻土温度场的精确计算具有重要意义。技术实现要素：本发明目是提供一种将冻土相变潜热转换为等效比热的计算方法，以利于冻土等效比热的计算。为实现上述目的，本发明提出一种将冻土相变潜热转换为等效比热的计算方法，该方法包括以下步骤：1)分别将冻土在[tj℃～tk℃]、[tk℃～ti℃](0≥tj＞tk＞ti)的潜热qj～k、qk～i用式(1)和式(2)表示出来，式(1)和式(2)为：式中，qj～k、qk～i分别为冻土在[tj℃～tk℃]、[tk℃～ti℃](0≥tj＞tk＞ti)的潜热，kj/kg；l为水凝结成冰释放的热量，取333.51kj/kg；ρd为土样的干密度，kg/m3；wuj、wuk、wui分别为温度tj、tk、ti时的冻土未冻水含量，％；ρ为土样的密度，kg/m3。2)分别将冻土在[tj℃～tk℃]、[tk℃～ti℃](0≥tj＞tk＞ti)的潜热转化成的土体比热用式(3)和式(4)表示出来，式(3)和式(4)为：式中，cj'～k、c'k～i分别为冻土在[tj℃～tk℃]、[tk℃～ti℃](0≥tj＞tk＞ti)的潜热转化成的土体比热，kj/(kg·℃)；l为水凝结成冰释放的热量，取333.51kj/kg；ρd为土样的干密度，kg/m3；wuj、wuk分别为温度tj、tk时刻的冻土未冻水含量，％；ρ为土样的密度，kg/m3；tj、tk、ti分别为冻土的温度，℃。3)则tk℃的冻土比热可由公式(5)计算出来，公式(5)为：式中，c'k为tk℃的冻土等效比热，kj/(kg·℃)；ck为tk℃的冻土真比热，kj/(kg·℃)；cj'～k、c'k～i分别为冻土在[tj℃～tk℃]、[tk℃～ti℃](0≥tj＞tk＞ti)的潜热转化成的土体比热，kj/(kg·℃)。依据该方法可将冻土的真比热转化为冻土的等效比热，为冻土温度场的计算奠定基础。本发明的效果是提供了能够将冻土相变潜热转换为等效比热的计算方法，提供的方法考虑了冻土等效比热与平均比热、真比热的不同。通过实测和计算表明，该计算方法的预测精度较传统方法提升了84.7％。精度的提高能够最大限度的将冻土相变潜热转换为等效比热，为冻土温度场的计算提供便利。具体实施方式对本发明的一种将冻土相变潜热转换为等效比热的计算方法加以说明。本发明的一种将冻土相变潜热转换为等效比热的计算方法原理：通过求取两个相邻温度区间的冻土相变潜热，从而计算获取两个相邻温度区间的冻土比热，在此基础上，通过对两个相邻温度区间的比热求平均值并与两个相邻温度区间分界点的真比热相加的方法，获取两个相邻温度区间分界点的冻土的等效比热。本发明提出一种将冻土相变潜热转换为等效比热的计算方法，该方法包括以下步骤：1)分别将冻土在[tj℃～tk℃]、[tk℃～ti℃](0≥tj＞tk＞ti)的潜热qj～k、qk～i用式(6)和式(7)表示出来，式(6)和式(7)为：式中，qj～k、qk～i分别为冻土在[tj℃～tk℃]、[tk℃～ti℃](0≥tj＞tk＞ti)的潜热，kj/kg；l为水凝结成冰释放的热量，取333.51kj/kg；ρd为土样的干密度，kg/m3；wuj、wuk、wui分别为温度tj、tk、ti时的冻土未冻水含量，％；ρ为土样的密度，kg/m3。2)分别将冻土在[tj℃～tk℃]、[tk℃～ti℃](0≥tj＞tk＞ti)的潜热转化成的土体比热用式(8)和式(9)表示出来，式(8)和式(9)为：式中，cj'～k、c'k～i分别为冻土在[tj℃～tk℃]、[tk℃～ti℃](0≥tj＞tk＞ti)的潜热转化成的土体比热，kj/(kg·℃)；l为水凝结成冰释放的热量，取333.51kj/kg；ρd为土样的干密度，kg/m3；wuj、wuk分别为温度tj、tk时刻的冻土未冻水含量，％；ρ为土样的密度，kg/m3；tj、tk、ti分别为冻土的温度，℃。3)则tk℃的冻土比热可由公式(10)计算出来，公式(10)为：式中，c'k为tk℃的冻土等效比热，kj/(kg·℃)；ck为tk℃的冻土真比热，kj/(kg·℃)；cj'～k、c'k～i分别为冻土在[tj℃～tk℃]、[tk℃～ti℃](0≥tj＞tk＞ti)的潜热转化成的土体比热，kj/(kg·℃)。依据该方法可将冻土的真比热转化为冻土的等效比热，为冻土温度场的计算奠定基础。实施例通过试验确定干密度为1.8g/cm3，含水量为15.8％，ds为2.72的粉质粘土试样在-1℃、-2℃、-3℃的未冻水含量，通过混合量热试验确定所述干密度为1.8g/cm3，ds为2.72的粉质粘土试样中的土颗粒的比热和含量，以及土中各相的质量含量。依据加权平均法获取冻土真比热ck为1.338kj/kg·℃，即为传统方法计算值。利用提出的一种将冻土相变潜热转换为等效比热的计算方法，计算见表1所示。计算获取了-2℃摄氏度下的冻土比热c'k为2.094kj/kg·℃。表1某冻土相变潜热转换为等效比热的计算项wu1.5wu2wu2.5q1.5～2q2～2.5c’1.5～2c’2～2.5ckc’k单位％％％kj/kgkj/kgkj/(kg·℃)kj/(kg·℃)kj/(kg·℃)kj/(kg·℃)数值9.6919.6909.6880.2630.7010.5261.4021.1382.102经对比传统方法的计算值ck和本发明的一种将冻土相变潜热转换为等效比热的计算方法的计算值c'k，发现本发明提供的计算方法较传统方法的预测精度提升了84.7％。以上所述仅为结合本次计算过程进行说明，对于本领域的技术人员来说，本发明可有各种变化和更改。凡依据本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。当前第1页12