宋利红 杨 宇 李海福 孙晗森 郑丽婧 徐春阳

1. 成都理工大学能源学院 2. 中海油研究总院有限责任公司 3. 中国石油辽河油田公司勘探开发研究院

0 引言

地热储层温度的确定在地热资源的开发评价过程中有着举足轻重的地位，众多学者[1-4]针对不同地区利用阳离子Na-K温标、K-Mg温标、Na-K-Ca温标、Na-K-Mg温标等对热储温度进行了计算。但地热储层温度通常会受到很多因素的影响：不同储层的地热流体可能在采出前发生混合，对流体温度进行分析时存在一定的误差；相邻层段的层间窜流、循环对流等会使地热流体在开采前就与其他液体发生混合，从而影响储层温度的测量。特别是地面冷水与热流体的混合会使大量的Cl-、Na+、K+等离子混入，这些离子之间的比值会大大影响地热温标方法计算的准确性。因此，针对这种情况，常常利用混合模型计算热储层温度、分析地下热流体深部温度特征和冷热水混合过程[5]。

在1974年，Fournier和Truesdell根据硅-焓（即二氧化硅-水的焓值）平衡提出了混合模型并成功对储层原始流体的温度进行了估算；王晓翠[6]对云南沘江沿岸温泉水进行了水化学分析，采用硅-焓图解法对地下热水在采出前未发生蒸汽和热量损失的情况进行了热储温度计算。但前人的实践应用表明，硅-焓图解法对热储温度的计算已经广泛应用，但绘图过程繁琐。基于此，在硅-焓图解法的基础上，提出改进图解法以及Newton-Raphson迭代算法，并对热储温度进行了计算。

1 混合模型的提出

地表的大气降水通常会与地热水混合， 通常的混合情况如图1所示。……