论文部分内容阅读
油气资源开采、核废料处置等工程的安全稳定与岩体渗流及传热特性密切相关。天然岩体中大量的节理、裂隙、断层是地下水流动的主要通道，对工程岩体渗透特性起到重要作用。因此，深入研究结构面及其构成裂隙网络的渗流及传热行为有着重要的理论及工程意义。
为了系统地研究裂隙岩体渗流及传热特性，本文以裂隙岩体为研究对象，运用理论、试验和数值模拟方法，开展了粗糙裂隙和相交裂隙渗流试验，建立了裂隙水热耦合模型，提出了三维裂隙网络水热耦合计算方法并应用于地热储层开采问题。主要研究内容如下：
(1)自主设计了一套粗糙单裂隙渗流可视化模拟试验装置，并预制不同粗糙度的裂隙试样。通过开展一系列不同水压下单裂隙渗流试验，结合数值方法研究了粗糙裂隙的非线性渗流特性。结果表明，Forchheimer方程可以描述粗糙裂隙非线性渗流行为。非线性参数与裂隙几何形状有关，非线性参数的范围从6e10~5e15不等，相差5个数量级。同时，建立了非线性参数与水力开度、接触率和分形维度的半经验方程，并基于此，得到了一个新的临界雷诺数的半经验公式。分析可知，临界雷诺数随着接触率和分形维度的增加而减少，取值范围在0~60之间。另外，基于裂隙单元流速与其亚裂隙组单元开度成正比的假定，建立了三维粗糙裂隙非线性渗流数学模型并得到了渗流试验验证。
(2)基于严格的数学理论，利用楔形裂隙流体流动和传热半解析解，得到了考虑粗糙度和开度变化的流体与岩石裂隙之间的水热耦合模型。根据三维光学扫描仪获得真实的裂隙，基于该模型研究了裂隙粗糙度和注入水压力对水流传热特性的影响。结果表明，裂隙表面粗糙度对局部传热系数分布影响较大，但对温度分布影响较小。水压会改变温度和局部传热系数大小，而不改变其分布规律。同时也到了裂隙分形维度、水压和水力开度与平均传热系数的关系。随着分形维度的增加，平均传热系数减小；水压的增加导致平均传热系数增加；而水力开度与平均传热系数呈负相关。
(3)分别预制了不同交叉形式和交叉角度的岩石裂隙试样，包括平直裂隙、弯折裂隙、交叉裂隙和分叉裂隙。对四种相交裂隙进行了一系列不同水压下的渗流试验，研究了相交裂隙非线性渗流行为，同时进行了相应的数值模拟分析。试验和的数值模拟结果表明，对于弯折裂隙，随着弯折角α的增加，Forchheimer系数β减少而临界水力梯度Jc增加。相反，对于交叉裂隙，随着交叉角θ的增加，Forchheimer系数β呈现非线性增加而临界水力梯度Jc呈现非线性减小。对于分叉裂隙，非线性渗流特性取决于分叉角γ1和γ2的值。|γ1-γ2|越大，Forchheimer系数β越大的而临界水力梯度Jc越小。另外，分别得到了三种类型的相交裂隙Forchheimer方程非线性系数半经验关系式。
(4)基于三维离散裂隙网络，提出了雷诺方程和Forchheimer方程相结合的有限元计算方法，并用直接迭代法求解大规模非线性方程组。利用该方法模拟了三维裂隙网络渗流行为，探讨了裂隙粗糙度和连通性对裂隙网络非线性渗流特性的影响。结果表明，裂隙粗糙度越大，即使离散裂隙网络连通性较好，也容易产生优势的流动路径，进一步发生“沟槽流动效应”。随着水力梯度的增大，离散裂隙网络的导水系数逐渐减小。另外，导水系数随连通度的增大而增大，随裂隙表面粗糙度的增大而减小。最后分析了粗糙度和连通度对裂隙网络的临界水力梯度的影响，结果表明，在裂隙表面较光滑、连通性较差的模型中，临界水力梯度较低。临界水力梯度范围为0.001~0.5之间。
(5)基于热传递模型和非线性渗流模拟方法，构建了一个基于裂隙网络模型的水热耦合计算方法并验证了其可靠性。将该方法应用于澳大利亚Habanero地热开采工程案例，计算结果与现场监测成果较吻合。在注入井和生产井位置，出现了高流速区域，大约是其他位置速度的十倍以上，该位置可能出现湍流现象。在靠近注入井的位置，温度急剧增加；当距离达到100m左右时，温度平稳，达到岩石温度483.15K。注入压力会影响储层地热开采寿命，压力越大时，生产井维持初始温度的时间增加，储层开采寿命越大。
为了系统地研究裂隙岩体渗流及传热特性，本文以裂隙岩体为研究对象，运用理论、试验和数值模拟方法，开展了粗糙裂隙和相交裂隙渗流试验，建立了裂隙水热耦合模型，提出了三维裂隙网络水热耦合计算方法并应用于地热储层开采问题。主要研究内容如下：
(1)自主设计了一套粗糙单裂隙渗流可视化模拟试验装置，并预制不同粗糙度的裂隙试样。通过开展一系列不同水压下单裂隙渗流试验，结合数值方法研究了粗糙裂隙的非线性渗流特性。结果表明，Forchheimer方程可以描述粗糙裂隙非线性渗流行为。非线性参数与裂隙几何形状有关，非线性参数的范围从6e10~5e15不等，相差5个数量级。同时，建立了非线性参数与水力开度、接触率和分形维度的半经验方程，并基于此，得到了一个新的临界雷诺数的半经验公式。分析可知，临界雷诺数随着接触率和分形维度的增加而减少，取值范围在0~60之间。另外，基于裂隙单元流速与其亚裂隙组单元开度成正比的假定，建立了三维粗糙裂隙非线性渗流数学模型并得到了渗流试验验证。
(2)基于严格的数学理论，利用楔形裂隙流体流动和传热半解析解，得到了考虑粗糙度和开度变化的流体与岩石裂隙之间的水热耦合模型。根据三维光学扫描仪获得真实的裂隙，基于该模型研究了裂隙粗糙度和注入水压力对水流传热特性的影响。结果表明，裂隙表面粗糙度对局部传热系数分布影响较大，但对温度分布影响较小。水压会改变温度和局部传热系数大小，而不改变其分布规律。同时也到了裂隙分形维度、水压和水力开度与平均传热系数的关系。随着分形维度的增加，平均传热系数减小；水压的增加导致平均传热系数增加；而水力开度与平均传热系数呈负相关。
(3)分别预制了不同交叉形式和交叉角度的岩石裂隙试样，包括平直裂隙、弯折裂隙、交叉裂隙和分叉裂隙。对四种相交裂隙进行了一系列不同水压下的渗流试验，研究了相交裂隙非线性渗流行为，同时进行了相应的数值模拟分析。试验和的数值模拟结果表明，对于弯折裂隙，随着弯折角α的增加，Forchheimer系数β减少而临界水力梯度Jc增加。相反，对于交叉裂隙，随着交叉角θ的增加，Forchheimer系数β呈现非线性增加而临界水力梯度Jc呈现非线性减小。对于分叉裂隙，非线性渗流特性取决于分叉角γ1和γ2的值。|γ1-γ2|越大，Forchheimer系数β越大的而临界水力梯度Jc越小。另外，分别得到了三种类型的相交裂隙Forchheimer方程非线性系数半经验关系式。
(4)基于三维离散裂隙网络，提出了雷诺方程和Forchheimer方程相结合的有限元计算方法，并用直接迭代法求解大规模非线性方程组。利用该方法模拟了三维裂隙网络渗流行为，探讨了裂隙粗糙度和连通性对裂隙网络非线性渗流特性的影响。结果表明，裂隙粗糙度越大，即使离散裂隙网络连通性较好，也容易产生优势的流动路径，进一步发生“沟槽流动效应”。随着水力梯度的增大，离散裂隙网络的导水系数逐渐减小。另外，导水系数随连通度的增大而增大，随裂隙表面粗糙度的增大而减小。最后分析了粗糙度和连通度对裂隙网络的临界水力梯度的影响，结果表明，在裂隙表面较光滑、连通性较差的模型中，临界水力梯度较低。临界水力梯度范围为0.001~0.5之间。
(5)基于热传递模型和非线性渗流模拟方法，构建了一个基于裂隙网络模型的水热耦合计算方法并验证了其可靠性。将该方法应用于澳大利亚Habanero地热开采工程案例，计算结果与现场监测成果较吻合。在注入井和生产井位置，出现了高流速区域，大约是其他位置速度的十倍以上，该位置可能出现湍流现象。在靠近注入井的位置，温度急剧增加；当距离达到100m左右时，温度平稳，达到岩石温度483.15K。注入压力会影响储层地热开采寿命，压力越大时，生产井维持初始温度的时间增加，储层开采寿命越大。