第３８卷第７期 开发工程 ６５ 考虑多重滤失效应的前置液酸压有效缝长模拟 党录瑞 周长林 黄媚 蒋德生 １．中国石油西南油气田公司２．中国石油西南油气田公司工程技术研究院 摘要为了提高天然裂缝较发育的碳酸盐岩储层酸化压裂设计的准确性和有效性，在采用经典裂缝拟三维延伸数学模型模拟 前置液造缝过程中裂缝几何尺寸动态变化的基础上，结合液相反应平衡原理和局部反应平衡原理，建立了考虑天然裂缝、蚓孔及基 质多重滤失效应的酸液流动反应模型，以现场实例井的酸化压裂改造为例，使用所建模型模拟前置液酸压注酸过程中裂缝内酸液渗滤 过程及酸岩动态刻蚀形态，通过酸液沿水力裂缝长度方向的浓度变化和残酸极限浓度综合确定酸蚀裂缝有效缝长，并与压力恢复试井的 解释结果进行了对比。结果表明：①裂缝一孔隙型储层酸化压裂过程中，沿裂缝长度方向的酸液滤失速度非定值，滤失速度曲线呈 锯齿状波动变化，在酸蚀蚓孔与天然裂缝相遇时酸液滤失速度普遍大于基质处的滤失速度；②裂缝一孔隙型储层酸化压裂过程中酸 液滤失严重，酸液有效作用距离大幅度降低；③验证井酸化压裂解堵增产效果显著，该井试井解释结果与所建立的数学模型模拟解 释结果相吻合，说明所建立的前置液酸化压裂数学模型可靠。结论认为，考虑多重滤失效应的前置液酸化压裂数学模型更适用于裂缝～ 孔隙型储层的酸化压裂模拟。 关键词碳酸盐岩储集层 天然裂缝蚓孔 前置液酸化压裂酸液滤失有效缝长 裂缝一孔隙型 酸化压裂模拟 Ｄ０Ｉ：１０．３７８７￣．ｉｓｓｎ．１０００．０９７６．２０１８．０７．００９ Ｓｉｍｕｌａｔｉｏｎ ｏｆ ｅｆｆｅｃｔｉｖｅ ｆｒａｃｔｕｒｅ ｌｅｎｇｔｈ ｏｆ ｐｒｅｐａｄ ａｃｉｄ ｆｒａｃｔｕｒｉｎｇ ｃｏｎｓｉｄｅｒｉｎｇ ｍｕｌｔｉ— ｐｉｅ ｌｅａｋ—Ｏｆ１ｆｅｆｆｅｃｔ Ｄａｎｇ Ｌｕｍｉ ，Ｚｈｏｕ Ｃｈａｎｇｌｉｎ ，Ｈｕａｎｇ Ｍｅｉ ＆Ｊｉａｎｇ Ｄｅｓｈｅｎｇ ｛１．ＰｅｔｒｏＣｈｉｎａ Ｓｏｕｔｈｗｅｓｔ Ｏｉｌ＆Ｇａｓｉｅｌｆｄ Ｃｏｍｐａｎｙ，Ｃｈｅｎｇｄｕ，Ｓｉｃｈｕａｎ ６１００５１ Ｃｈｉｎａ；２．Ｅｎｇｉｎｅｅｒｉｎｇ Ｔｅｃｈｎｏｌｏｇｙ Ｒｅ— ｓｅａｒｃｈ Ｉｎｓｔｉｔｕｔｅ，ＰｅｔｒｏＣｈｉｎａ Ｓｏｕｔｈｗｅｓｔ Ｏｉｌ＆Ｇａｓｉｆｅｌｄ Ｃｏｍｐａｎｙ，Ｃｈｅｎｇｄｕ，Ｓｉｃｈｕａｎ ６１００３１，Ｃｈｉｎａ） ＮＡＴｕＲ．ＧＡＳ ＩＮＤ．ＶＯＬＵＭＥ ３８，ＩＳＳＵＥ ７，ＰＰ．６５—７２，７／２５／２０ １ ８．（ＩＳＳＮ ｌ０００．０９７６；Ｉｎ Ｃｈｉｎｅｓｅ） Ａｂｓｔｒａｃｔ：Ｆｏｒ ｔｈｅ ｐｕｒｐｏｓｅ ｏｆ ｉｍｐｒｏｖｉｎｇ ｔｈｅ ａｃｃｕｒａｃｙ ａｎｄ ｅｆｆｅｃｔｉｖｅｎｅｓｓ ｏｆ ａｃｉｄ ｆｒａｃｔｕｒｉｎｇ ｄｅｓｉｇｎ ｆｏｒ ｃａｒｂｏｎａｔｅ ｒｅｓｅｒｖｏｉｒｓ ｗｉｔｈ ｄｅｖｅｌｏｐｅｄ ｎａｔｕｒａｌ ｆｒａｃｔｕｒｅｓ，ａｎ ａｃｉｄｉｚｉｎｇ ｆｌｕｉｄ ｆｌｏｗ ａｎｄ ｒｅａｃｔｉｏｎ ｍｏｄｅｌ ｔａｋｉｎｇ ｔｈｅ ｍｕｌｔｉｐｌｅ ｌｅａｋ．ｏｆｆ ｅｆｆｅｃｔ ｏｆ ｎａｔｕｒａｌ ｆｒａｃｔｕｒｅｗｏｒｍｈｏｌｅ ａｎｄ ｍａｔｒｉｘ ．ｉｎｔｏ ａｃｃｏｕｎｔ ｗａｓ ｅｓｔａｂｌｉｓｈｅｄ ａｃｃｏｒｄｉｎｇ ｔｏ ｔｈｅ ｌｉｑｕｉｄ ｐｈａｓｅ ｒｅａｃｔｉｏｎ ｅｑｕｉｌｉｂｒｉｕｍ ｐｒｉｎｃｉｐｌｅ ａｎｄ ｔｈｅ ｌｏｃａｌ ｒｅａｃｔｉｏｎ ｅｑｕｉｌｉｂｒｉｕｍ ｐｒｉｎｃｉｐｌｅ ａｆｔｅｒ ｔｈｅ ｄｙｎａｍｉｃ ｃｈａｎｇｅ ｏｆ ｆｒａｃｔｕｒｅ ｇｅｏｍｅｔｒｙ ｉｎ ｔｈｅ ｐｒｏｃｅｓｓ ｏｆ ｆｒａｃｔｕｒｅ ｃｒｅａｔｉｎｇ ｂｙ ｐｒｅｐａｄ ｆｌｕｉｄ ｗａｓ ｓｉｍｕｌａｔｅｄ ｉｎ ｔｈｅ ｃｌａｓｓｉｃａｌ ｐｓｅｕｄｏ—ｔｈｒｅｅ ｄｉｍｅｎｓｉｏｎａｌ ｍａｔｈｅｍａｔｉｃａｌ ｍｏｄｅｌ ｏｆ ｆｒａｃｔｕｒｅ ｐｒｏｐａｇａｔｉｏｎ．Ｔｈｅｎ，ｔｈｅ ａｃｉｄ ｆｒａｃｔｕｒｉｎｇ ｓｔｉｍｕｌａｔｉｏｎ ｏｆ ａ ｃａｓｅ ｗｅｌｌ ｏｎ ｓｉｔｅ ｗａｓ ｔａｋｅｎ ａｓ ａｎ ｅｘａｍｐｌｅ．Ｔｈｅ ｎｅｗｌｙ ｄｅｖｅｌｏｐｅｄ ｍｏｄｅｌ ｗａｓ ｕｓｅｄ ｔｏ ｓｉｍｕｌａｔｅ ｔｈｅ ｆｉｌｔｒａｔｉｏｎ ｐｒｏｃｅｓｓ ｏｆ ａｃｉｄｉｚｉｎｇ ｆｌｕｉｄ ｉｎ ｆｒａｃｔｕｒｅｓ ａｎｄ ｔｈｅ ｄｙｎａｍｉｃ ｅｔｃｈｉｎｇ ｍｏｒｐｈｏｌｏｇｙ ｏｆ ａｃｉｄｉｃ ｒｏｃｋｓ ｄｕｒｉｎｇ ｔｈｅ ａｃｉｄｉｚｉｎｇ ｆｌｕｉｄ ｉｎｊｅｃｔｉｏｎ ｏｆ ｐｒｅｐａｄ ａｃｉｄ ｆｒａｃｔｕｒｉｎｇ．Ｔｈｅ ｅｆｆｅｃｔｉｖｅ ｌｅｎｇｔｈ ｏｆ ｅｔｃｈｅｄ ｆａｃｔｕｒｅｓ ｗａｓ ｄｅｔｅｒｍｉｎｅｄ ｂｙ ａｎａｌｙｚｉｎｇ ｔｈｅ ｃｏｎｃｅｎｔｒａｔｉｏｎ ｃｈａｎｇｅ ｏｆ ａｃｉｄｉｚｉｎｇ ｆｌｕｉｄ ａｌｏｎｇ ｔｈｅ ｄｉｒｅｃｔｉｏｎ ｏｆ ｈｙｄｒａｕｌｉｃ ｆｒａｃｔｕｒｅ ｌｅｎｇｔｈ ａｎｄ ｔｈｅ ｔｈｒｅｓｈｏｌｄ ｃｏｎｃｅｎｔｒａｔｉｏｎ ｏｆ ｒｅｓｉｄｕａｌ ａｃｉｄｉｚｉｎｇ ｆｌｕｉｄ ｃｏｍｐｒｅｈｅｎｓｉｖｅ１ ａｎｄ ｔｈｅｎ ｉｔ ｗａｓ ｃｏｍｐａｒｅｄ ｗｉｔｈ ｔｈｅ ｉｎｔｅｒｐｒｅｔａｔｉｏｎ ｒｅｓｕｌｔ ｏｆ ｐｒｅｓｓｕｒｅ ｂｕｉｌｄｕｐ ｔｅｓｔ．Ａｎｄ ｔｈｅ ｆｏｌｌｏｗｉｎｇ ｒｅｓｅａｒｃｈ ｒｅｓｕｌｔｓ ｗｅｒｅ ｏｂｔａｉｎｅｄ．Ｆｉｒｓｔ，ｉｎ ｔｈｅ ｐｒｏｃｅｓｓ ｏｆ ａｃｉｄ ｆｒａｃｔｕｒｉｎｇ ｉｎ ｆｒａｃｔｕｒｅｄ－ｐｏｒｏｕｓ ｒｅｓｅｒｖｏｉｒｓ，ｔｈｅ ａｃｉｄｉｚｉｎｇ ｌｕｉｆｄ ｆｉｌｔｒａｔｉｏｎ ｖｅｌｏｃｉｔｙ ａｌｏｎｇ ｔｈｅ ｄｉｒｅｃｔｉｏｎ ｏｆ ｆｒａｃｔｕｒｅ ｌｅｎｇｔｈ ｉｓ ｎｏｔ ｃｏｎｓｔａｎｔ ａｎｄ ｔｈｅ ｆｉｌｔｒａｔｉｏｎ ｖｅｌｏｃｉｔｙ ｃｕｒｖｅ ｆｌｕｃｔｕａｔｅｓ ｉｎ ａ ｓｅｒｒａｔｅｄ ｓｈａｐｅ． Ａｎｄ ｔｈｅ ａｃｉｄｉｚｉｎｇ ｆｌｕｉｄ ｆｉｌｔｒａｔｉｏｎ ｖｅｌｏｃｉｔｙ ｗｈｅｒｅ ｅｔｃｈｅｄ ｗｏｒｍｈｏｌｅｓ ｍｅｅｔ ｎａｔｕｒａｌ ｆｒａｃｔｕｒｅｓ ｉｓ ｃｏｍｍｏｎｌｙ ｈｉｇｈｅｒ ｔｈａｎ ｔｈａｔ ｉｎ ｍａｔｒｉｘＳｅｃｏｎｄ， ．ａｃｉｄｉｚｉｎｇ ｆｌｕｉｄ ｉｓ ｌｏｓｔ ｓｅｒｉｏｕｓｌｙ ａｎｄ ｔｈｅ ｅｆｆｅｃｔｉｖｅ ｄｉｓｔａｎｃｅ ｏｆ ａｃｉｄｉｚｉｎｇ ｆｌｕｉｄ ｇｅｔｓ ｓｈｏｒｔ ｓｉｇｎｉｆｉｃａｎｔｌｙ ｄｕｒｉｎｇ ｔｈｅ ａｃｉｄ ｆｒａｃｔｕｒｉｎｇ ｏｆ ｆｒａｃｔｕｒｅｄ— ｐｏｒｏｕｓ ｒｅｓｅｒｖｏｉｒｓ．Ｔｈｉｒｄ，ａｃｉｄ ｆｒａｃｔｕｒｉｎｇ ｉｎ ｔｈｅ ｖｅｒｉｆｉｃａｔｉｏｎ ｗｅｌｌ ｉｓ ｒｅｍａｒｋａｂｌｅ ｉｎ ｂｌｏｃｋａｇｅ ｒｅｍｏｖｉｎｇ ａｎｄ ｓｔｉｍｕｌａｔｉｏｎ，ａｎｄ ｉｔｓ ｗｅｌｌ ｔｅｓｔ ｉｎｔｅｒｐｒｅｔａｔｉｏｎ ｒｅｓｕｌｔｓ ａｒｅ ｃｏｎｓｉｓｔｅｎｔ ｗｉｔｈ ｔｈｅ ｓｉｍｕｌａｔｉｏｎ ｉｎｔｅｒｐｒｅｔａｔｉｏｎ ｒｅｓｕｌｔｓ ｐｒｏｖｉｄｅｄ ｂｙ ｔｈｅ ｎｅｗｌｙ ｄｅｖｅｌｏｐｅｄ ｍａｔｈｅｍａｔｉｃａｌ ｍｏｄｅ１Ｉｔ ｉｓ ．ｉｎｄｉｃａｔｅｄ ｔｈａｔ ｔｈｉｓ ｎｅｗｌｙ ｄｅｖｅｌｏｐｅｄ ｍｏｄｅｌ ｉｓ ｒｅｌｉａｂｌｅ．Ｉｎ ｃｏｎｃｌｕｓｉｏｎ，ｔｈｅ ｍａｔｈｅｍａｔｉｃａｌ ｍｏｄｅｌ ｏｆ ｐｒｅｐａｄ ａｃｉｄ ｆｒａｃｔｕｒｉｎｇ ｗｈｉｃｈ ｃｏｎｓｉｄｅｒｓ ｍｕｌｔｉｐｌｅ ｌｅａｋ－ｏｆｆ ｅｆｆｅｃｔ ｉｓ ｍｏｒｅ ｓｕｉｔａｂｌｅ ｆｏｒ ａｃｉｄ ｆｒａｃｔｕｒｉｎｇ ｓｉｍｕｌａｔｉｏｎ ｏｆ ｒａｃｔｕｒｅｄ－ｐｏｒｏｕｓ ｒｅｓｅｒｖｏｉｆｒｓ． Ｋｅｙｗｏｒｄｓ：Ｃａｒｂｏｎａｔｅ ｒｏｃｋ；Ｒｅｓｅｒｖｏｉｒ；Ｎａｔｕｒａｌ ｆａｃｔｕｒｅ；Ｗｏｒｍｈｏｌｅ；Ｐｒｅｐａｄ ａｃｉｄ ｆｒａｃｔｕｒｉｎｇ；Ａｃｉｄｉｚｉｎｇ ｆｌｕｉｄ ｆｉｌｔｒａｔｉｏｎ；Ｅｆｆｅｃｔｉｖｅ ｆｒａｃｔｕｒｅ ｌｅｎｇｔｈ；Ｆｒａｃｔｕｒｅｄ—ｐｏｒｏｕｓ；Ａｃｉｄ ｆｒａｃｔｕｒｉｎｇ ｓｉｍｕｌａｔｉｏｎ 基金项目：国家科技重大专项“四川盆地大型碳酸盐岩气田开发示范工程”（编号：２０１６ＺＸ０５０５２）。 作者简介：党录瑞，１９６８年生，高级工程师：主要从事油气地质研究及科技管理工作。地址：（６１００５１）四川省成都市府青路一段３号。 ＯＲＣＩＤ：００００·０００１—６７９０－３０７７。Ｅ－ｍａｉｌ：ｄａｎｇ＿ｌｒ￣２ｐｅｔｒｏｃｈｉｎａ．ｃｏｍ．ｃｎ ６６ 天 然气 工业 ２０１８年７月 ０ 引言 酸液有效作用距离是指酸液在裂缝内的流动反 方向浓度的变化和残酸极限浓度综合确定酸蚀裂缝 的有效缝长。 应过程中，当酸液浓度降低到一定程度，基本失去 １考虑多重滤失效应的前置液酸压数 溶蚀能力时所流经的距离［１－２］。通常取酸液初始注入 浓度的１０％作为残酸极限浓度，以此确定酸液有效 作用距离［３－４］。准确预测酸液有效作用距离对掌控酸 压施工效果，精准设计酸压施工参数及泵注程序有 着十分重要的影响［５－７１。 以往的研究学者针对酸液有效作用距离做了 大量的研究，普遍都是根据酸液浓度分布模型结合 残酸极限浓度来确定酸液有效作用距离。国外学者 Ｒｏｂｅｒｔｓ和Ｇｕｉｎ【８ 将酸液滤失速度考虑为常数，建立 二维稳态酸液浓度分布模型。Ｒｏｍｅｒｏ等 建立了三 维非稳态酸液浓度分布模型，并耦合ＫＧＤ几何尺寸 模型，与二维模型相比，其酸液有效作用距离的计 算结果更准确。国内学者以任书泉［１０１赵碧华［１１－１２］、 李平［１ 等为代表着重讨论了同离子效应、温度、反 应生成热、物性参数及施工参数等对酸液流动反应 的影响。任书泉在酸液流动反应影响因素研究的基 础上，建立了无滤失酸液有效作用距离计算图版， 提出增大酸液有效作用距离的措施方法［Ｉ 。侯博恒 在前人研究的数学模型基础上，进一步考虑了酸压 滤失及酸岩反应产物ＣＯ，等对酸液有效作用距离的 影响，并改进了酸液有效作用距离预测方法，提高 了计算的合理性ｕ 。 国内外的研究学者尽管针对酸液有效作用距离 做了大量研究，但大都是针对孔隙性储层提出的， 很少考虑酸液经天然裂缝和酸蚀蚓孔所引起的滤失。 在采用通用的商业软件进行酸压模拟时，由于未考 虑酸液在蚓孑Ｌ与天然裂缝中滤失的影响，模拟出的 酸液有效距离通常比实际结果大，不能很好地指导 生产实践。且试验表明在裂缝一孔洞型储层中酸液 滤失量会急剧增加［１６－１８］对酸液有效作用距离的影 响是不可忽视的。因此，在酸压施工设计中必须考 虑到酸液在天然裂缝及蚓孔中的滤失对酸液有效作 用距离的影响，才能针对天然裂缝较发育的碳酸盐 岩储层提高酸压设计的准确性。为此，在采用经典 的裂缝拟三维延伸数学模型模拟前置液造缝过程中 裂缝的几何尺寸动态变化基础上，结合液相反应平 衡原理和局部反应平衡原理，建立了考虑天然裂缝、 蚓孔及基质多重滤失效应的酸液流动反应模型，精 细模拟前置液酸压注酸过程中裂缝内酸液渗滤过程 及酸岩动态刻蚀形态，并通过酸液沿水力裂缝缝长 学模型 假设条件如下：①储层、盖层、底层均为连续、 各向同性、均质的线弹性介质；②注入排量恒定； ③裂缝内流体沿一维方向流动；④不考虑温度对酸 液性质的影响。 １．１水力裂缝延伸模型 前置液酸压工艺是采用高黏非反应性前置压裂 液压开储层形成水力裂缝，然后高压挤入酸液刻蚀裂 缝表面，形成非均匀溶蚀缝面来增大裂缝的导流能 力。因此，计算酸蚀有效缝长，首先需要模拟前置液 造缝过程。本文采用经典的裂缝拟三维延伸数学模 型模拟前置液造缝过程中裂缝的几何尺寸动态变化， 并以此为基础开展下一步研究。 应力对称下裂缝拟三维延伸数学模型 ９Ｊ的连续 性方程为： 一 Ｏｑ（ｘ，ｔ）： ）＋ （１） 式中ｑ（ｘ， 表示ｔ时刻缝长 处的流量，ｍ ／ｍｉｎ；ｘ 表示缝长方向位置，１ＴＩ；ｆ表示施工时刻，ｒａｉｎ； ｌ 一 ，ｔ）表示ｔ时刻缝长 处的流体滤失量，ｍ２／ｍｉｎ；Ａ ，ｆ） 表示ｆ时刻缝长 处裂缝的横截面积，ｍ 。 流体滤失量（１，．）的计算式为： ２ｈｐＣ：———— ＝＝一（—Ｘ，ｔ） （２）芦一 Ｌ √ｆ—ｆ 式中ｈ 表示油层厚度，１１１；Ｃ（ｘ，力表示 时刻缝长 处的流体滤失系数，ｒｎ／ｍｉｎ ；，表示流体到达缝长 处的时间，ｍｉｎ。 简化的缝中流体压降方程为： ｆ（ ，０，ｔ）一 ６４／．ｉｆｑ（ｘ，ｆ） … （ ，ｔ）ｗｆ（ ，０， ） 式中ｐ ，０，ｔ）表示ｔ时刻缝长 、裂缝中心处压力， ＭＰａ； 表示压裂液黏度，ｍＰａ·Ｓ；ｈ（ｘ， 表示ｔ时 刻缝长 处的缝高，ｍ；ｗ ，０， 表示 时刻缝长 处 裂缝半宽度，ｍ。 裂缝宽度方程为： ： 寿 糌㈩ 第３８卷第７期 开发工程 其中 ，ｌ，ｒ一 其中 Ｖ ，， ：一 ＯＸ ＝ 式中Ｗ ，Ｚ，ｆ）表示在ｔ时刻缝长Ｘ、缝高ｚ处的裂缝 宽度，１１１；ｚ表示缝高方向位置，Ｉｌｌ；ｖ表示泊松比， 无量纲； 表示弹性模量，ＭＰａ； 表示无因次坐标， 无量纲；Ｕ表示高度积分变量，无量纲； 高ｚ处的净压力，ＭＰａ。 裂缝高度方程为： 一５×１０ ｗａ 印 ｏｌ２ ＾ Ｚｌ Ｋ Ｐｆ—Ｐ。 ，． ）表示缝 式中 ｆ表示裂缝Ｘ方向上的流速，ｍ／ｍｉｎ；Ｗａ表示 酸蚀裂缝宽度，ｍ； 表示酸液黏度，ｍＰａ·Ｓ； 表示裂缝ｚ方向上的流速，ｍ／ｍｉｎ；ｖ１ 表示垂直于裂 缝壁面的酸液动态滤失速度，ｒｒｇｍｉｎ；Ｋ表示岩心渗 篆 南ｃ（ ） Ｋ２ ｃ一 寿 ｄ冗～ ｘ㈦ 式中 表示应力强度因子，ＭＰａ·ｍ“ ；Ａｓ表示应 力差，ＭＰａ；ｆ表示油层厚度ｈ。与缝高ｈ（ｘ， ）之比， 无量纲。 １．２考虑多重滤失效应的酸液流动反应模型 基于水力裂缝几何尺寸，模拟注酸过程中裂缝 的动态刻蚀形态，并通过酸液浓度分布图和残酸极 限浓度可以确定酸蚀裂缝有效缝长。但是在高压挤 酸过程中，由于酸液会在垂直于裂缝壁面的方向上 滤失进入基质储层，形成酸蚀蚓孔，蚓孔的形成造 成酸液滤失，从而降低有效缝长［２ 。另外，由于碳 酸盐岩储层通常微裂缝发育，水力裂缝沟通微裂缝 后，酸液会进入微裂缝，触发酸液与微裂缝壁面发生 反应，从而造成酸液滤失。因此，注酸过程的模拟 应综合考虑天然裂缝及蚓孔中酸液的共同滤失效应， 才能更为准确地预测酸蚀裂缝有效作用距离。 在裂缝中酸液的流动反应物理模型如图１所示。 图１ 裂缝中酸液流动反应物理模型示意图 首先，通过物质平衡推导得到酸液在裂缝中的 流动、反应及滤失控制方程为 ： 纵 ＋ Ｏｚ ＋２ ：一掣 Ｕｔ ｃ６ 透率，ｍＤ；Ｐ，表示裂缝中压力，ＭＰａ：Ｐ。表示地层 压力，ＭＰａ；Ｗ 表示酸液滤失深度，ｍ。 酸压裂缝内的反应平衡方程为： ＯＸ ＋ ＋２Ｃｅ ＋２ （ｃ 一ｃｗ）： … ａ（Ｃ￣ｗａ） 式中Ｃｆ表示裂缝内的酸液浓度，ｋｍｏｌ／ｍ ： 表示 反应速度，ｍ／ｍｉｎ；Ｃｗ表示裂缝壁面的酸液浓度， ｋｍｏｌ／ｍ 。 裂缝壁面上的局部反应方程为： （Ｃｆ—Ｃｗ） （Ｃｗ） （８） 式中 （ ）表示单步不可逆反应的溶蚀速度，ｋｍｏｌ／ （ｍ ·ｍｉｎ）。 酸蚀裂缝宽度变化方程为： 喜 Ｅ２ Ｃｆ （Ｃｗ）］＝ 式中ｉ表示不同岩矿类型数； 表示ＨＣ１与灰岩／白 云岩间的溶解能力，ｋｇ／ｋｍｏｌ；Ｐ 表示灰岩／白云岩密 度，ｋｇ／ｍ。； 表示孔隙度，无因次； 表示滤失酸液 中与缝壁岩石发生反应的酸液占滤失酸液的百分数， 大多数情况下 ≈０。 １．３酸液滤失控制模型 在常规模型中，假设裂缝中酸液滤失速度仅与 基质渗透率有关，同时假设蚓孔长度为一恒定值口 。 但随着酸液不断注入，蚓孔不断延伸，加之天然裂 缝不断被刻蚀，酸液滤失加剧。因此，酸液滤失速 度与储层基质渗透率、天然裂缝分布及蚓孑Ｌ扩展密 切相关。综合考虑以上因素的影响，酸液滤失物理 模型如图２所示。 这里，引入双尺度模型旺 研究酸压过程中裂缝 壁面上酸液滤失引起的基质及天然裂缝中酸岩动态 反应过程，计算瞬态条件下人工裂缝壁面上不同位 ６８－ 天然气工业 ２０１８年７月 酸压裂缝 岩石溶解除了导致孔隙度发生变化，其他属性 参数也相应发生改变，Ｇａｒｍａｎ－－Ｋｏｚｅｎｙ模型口 给 出了孔隙度、渗透率、比表面和孔道半径间的微观 裂缝 图２酸液滤失物理模型示意图 一 Ｏ： ＋ （１０） ｔ ａｘ 其中 ：一 一ｏ（ＯＣｍ＇￣厂 ０（Ｏ【 Ｃｍ ）一 ㈦（ ）㈩～ 挈： （１３ 变化关系，此不赘述。 １．４初始条件和边界条件 １．４．１ 水力裂缝中酸液流动初始条件和边界条件 初始条件为： Ｐｆ（ ，ｚ，０）：Ｐｆ（ ，０，ｔｅｎｄ） （ ，ｚ，０）＝０ （１４） Ｗａ（ ，ｚ，０）＝Ｗｅ（ ，ｚ， ） 式中Ｐｆ ，０， ）表示前置液造缝结束时压力，ＭＰａ； 表示水力裂缝造缝阶段的完成时刻，ｍｉｎ；Ｗ ，Ｚ， ，ｅ ）表示前置液造缝结束时裂缝宽度，１ＴＩ。 边界条件为： Ｊｆ旦 ｌ出 ｈ１２ ／－／ａ Ｉ１．ｏ Ｐ （ ，ｚ）＝Ｐ 堕Ｉ ：０ ｒ １５） ０ｚ Ｉｚ＝－ｈ 监『：０ ｌ Ｃｆ（０，ｚ，ｔ）＝ 式中ｑｉ ｉ表示注酸排量，ｍ。／ｍｉｎ；Ｐｆ ， 表示裂缝尖 端处的压力，ＭＰａ；Ｌ表示裂缝长度，ｍ； 表示裂 缝入口处酸液浓度，ｋｍｏｌ／ｍ 。 １．４．２基质中酸液流动初始条件和边界条件 初始条件为： ｛Ｉ ＰＣｍ （ ｙ ，ｚ， Ｏ０ ０）＝ Ｐ。 （… １６） 边界条件为： （ ，１）＝ Ｐ （ ，１）＝Ｐ； Ｐ （ ，Ｍ）＝Ｐ。 ｌ ：０ （１７） ＆ Ｃｍ（ ，１）＝ ０ 式中１， 表示不同时刻，裂缝壁面上的滤失速度，ｍ／ 第３８卷第７期 ６９ ｍｉｎ：Ｐ ， 表示基质边界网格处压力，ＭＰａ；Ｍ 表１ Ｘ井酸压模拟输入参数表 参数 产层厚度／ｍ 表示 方向的网格节点个数；ｐ／表示不同时刻ｆ裂缝 内的净压力，ＭＰａ；ｃ／表示不同时刻ｔ裂缝内的酸液 浓度，ｋｍｏｌ／ｍ 。 数值 尺 产层应力／ＭＰａ 如 ５ ｋ ｌ÷ ０ ●２实例模拟与分析 以四川盆地高石梯一磨溪地区上震旦统灯影组 气藏ｘ井为例，该井酸压改造层为灯影组四段下亚 段，射孔井段为：５ ３８５～５ ３９９ ｍ、５ ４０５～５ ４０８ ｍ、 产层断裂韧性／（ＭＰａ·ｍ ） 产层弹性模量／ＭＰａ 产层泊松比 盖／底层应力／ＳＰａ 盖／底层弹性模量／ＭＰａ 盖／底层泊松比 渗透率／ｍＤ 发 育 ∞ ５ ４ｌ３．０～５ ４２０．５ ｍ，酸压改造段测井资料显示气层 为３层、差气层为１层，产层厚度２３．５ ｍ，孔隙度 介于２．０％～１３．６％，其中１号储层，常规测井资料 显示孔隙欠发育，电成像资料显示溶孔、裂缝欠发育： ２号、３号、４号储层，常规测井资料显示孔隙较发 育，电成像资料显示裂缝较发育。综合解释认为该 井酸压改造段储层为典型的裂缝一孔隙型储层。由 于ｘ井漏失严重，采用常规基质酸化技术难以解除 钻井液的伤害，酸压是该储层改造增产的必要手段。 采用射孔一酸压～测试联作工艺，共注入胶凝酸 ２３９．３６ ｍ ，其中高挤段（酸压段）用酸量为２２２．８６ ｍ ，施工排量介于３．８～４．９ ｍ ／ｍｉｎ，泵注压力（油 储层温度／￣Ｃ 酸液排量／（ｍ ·ｍｉｎ ） 酸液用量／ｍ３ 胶凝酸浓度 反应级数 溶蚀能力数／（ｋｇ·ｋｇ ） 反应活化能／（Ｊ·ｍｏｌ ） 传质系数／（ｒｎ ·Ｓ ） 频率因子／（ｍｏｌ·Ｌ ·ｓ ） 压）介于７２～９６ ＭＰａ。根据测井资料、岩心分析数据、 岩石力学参数、地应力室内实验成果及酸岩反应动 液４４ ｍｉｎ时，在天然裂缝不发育的情况下酸压有效 力学数据，建立ｘ井裂缝模型，输入参数如表１所示。 通过求解本文建立的数学模型得到不同时刻酸 液的浓度分布（图３）以及有效裂缝长度（图４）。注 ５Ｏ 缝长为４２ ｍ；在天然裂缝发育的情况下酸压有效缝 长为３０ ｍ。可见，酸液在天然裂缝中的滤失对酸压 效果的影响不可忽视。 ５Ｏ ３０ ３０ ｅ－，、一 一 ５Ｏ ３０ １０ 拗 一 ３墓ｌ 一 ６ １０ （ｍｏｌ·Ｌ ） 浓度／ ５０ ７０ ９０ 、１０ 要一１０ —蒌一１０ —娄一１０ —３０ ５０ ３０ ５０ ３０ ５０ ———ｌ０ ３０ ５０ ７Ｏ ９０ １０ ３０ 缝长／ｍ 缝 ／ｍ ｂ人然裂缝 ：发育（，：１０ ｍｉｎ） ５０ ３０ ５Ｏ ３０ １０ ５Ｏ ３０ １０ ｃ．天然裂缝不发育（ｒ：４４ ｍｉｎ） 三１０ 皇 一１０ ３０ ５０ 妻一１０ —ｉｚＥ－一１０ ３０ ５０ １０ ３Ｏ ５Ｏ ７Ｏ ９Ｏ —３０ ５０ １０ ３Ｏ ５Ｏ ７Ｏ ９０ ————缝Ｋ／ｍ 缝长／ｍ ｆ天然裂缝发育ｆ，＝４４ｍｉｎ） ｅ天然裂缝发育（ １０ｍｉｎ） 图３不同时刻酸液浓度分布云图 ７０ 对比水力裂缝形态和酸蚀裂缝形态（ ５） 矢【】， 酸化后裂缝宽度增加并且酸蚀裂缝非均匀展布。沿 缝长方向，由于酸液不断渗滤反应消耗，酸液浓度 逐渐降低、酸岩刻蚀能 也相应变差，表现出酸蚀 缝宽度随裂缝长度的增加而急剧降低。 峻液流过一 定距离之后，鲜酸变为残酸，不能继续：；ｆ１Ｊ蚀岩石壁埘。 本次在裂缝长度为４０ ＩＴＩ处，酸液基本火 反应能力， 酸作用后缝宽仍为水力裂缝宽度， ｉ 液时 ／ｍｉｎ 停 后朱经酸 液刻蚀的缝宽会逐渐闭合最后消失。 图４不同Ｂ，－ｔ￣ｌｌ有效裂缝长度曲线图 缝侥／Ｉｌｌ 距井筒２０ ｍ范围内酸岩刻蚀沟椭较深，酸液刻 ” Ｉ ｌ…）Ｉ ５ １ｌ Ｉ｝ ……ｌ 心｜＝＝＝＝１．１水力 缝缝 分ｎ ｂ腋蚀 绁 ， 分ｆ 图５ Ｘ井裂缝形态模拟云图 蚀增加缝宽明显，缝口处最大缝宽约３ ｃｍ，缝长２０ ＩＴＩ以外裂缝有一定程度刻蚀，但由于酸液有效浓度 较低，刻蚀程度逐渐变弱，结合酸液浓度分布图（图 ６）综合分析可知酸蚀裂缝长度３０ ｎｌ（残酸浓度取值 ÷ 为０．５ ｍｏｌ－Ｌ ）。 如图７所示，在酸蚀蚓孔及天然裂缝共同作用下， 酸压过程中沿缝长方向的酸液动态滤失速度并不为 定值，滤失速度曲线呈锯齿状波动变化（非光滑曲 线），在酸蚀蚓孔与天然裂缝相遇时酸液滤失速度普 遍大于基质处的滤失速度，裂缝壁面滤失速度介于 １０ ～１０＿。ｍ／ｍｉｎ，滤失速度较小，主要为面溶蚀 （）０２５ 量 缱Ｋ／ｎｌ 图６ Ｘ井酸液浓度分布图 （）（）２（） （）０１ ５ （）０ｌＯ （）００５ 缝【 ／Ｉ１１ Ｋ／［１１ ａ 虑失 ＩＪ【ｍ ｂ；０ｇ 、ｉ丛Ｊ堑 ｂ ｌ 图７ ｘ井裂缝壁面天然裂缝 蚓孑Ｌ滤失剖面及酸液滤失速度曲线图 第３８卷第７期 开发工程 形态。 ｘ井酸压后，测试产量为５４．２７×１０ ｎｌ ／ｄ，测 试油压为３２．８１ ＭＰａ，该井压力恢复试井解释选用有 限导流＋无限大边界模型，解释结果显示缝长为３１．９ ｍ，表皮系数为一５．５２，酸压解堵增产效果显著。与 建立的数学模型模拟解释结果吻合度较高，说明考 虑蚓孔、天然裂缝及基质多重滤失效应的前置液酸 压数学模型预测结果可靠。 ３结论 １）裂缝一孔隙型储层酸压过程中，沿裂缝长度 方向的酸液滤失速度并不为定值，滤失速度曲线呈 锯齿状波动变化，在酸蚀蚓孔与天然裂缝相遇时酸 液滤失速度普遍大于基质处的滤失速度。 ２）裂缝一孔隙型储层由于天然裂缝普遍发育， 加之酸蚀蚓孔增长及天然裂缝扩宽，酸液滤失严 重，酸液有效作用距离大幅降低。因此，在酸压模 拟中，酸液在酸蚀蚓孔及天然裂缝中的共同滤失不 可忽略。 ３）考虑蚓孔、天然裂缝及基质多重滤失效应的 前置液酸压数学模型解释的缝长与压力恢复试井解 释的结果吻合度较高，更适用于裂缝一孔隙型储层 的酸压模拟。 参考文献 【ｌ】薛衡，黄祖熹，赵立强，蒋卫东，刘平礼，梁冲．考虑岩矿非 均质性的前置液酸压模拟研究［Ｊ］．天然气工业，２０１８，３８（２）： ５９．６６． Ｘｕｅ Ｈｅｎｇ，Ｈｕａｎｇ Ｚｕｘｉ，Ｚｈａｏ Ｌｉｑｉａｎｇ，Ｊｉａｎｇ Ｗｅｉｄｏｎｇ，Ｌｉｕ Ｐｉｎｇｌｉ ＆Ｌｉａｎｇ Ｃｈｏｎｇ．Ａ ｓｉｍｕｌａｔｉｏｎ ｓｔｕｄｙ ｏｎ ｔｈｅ ｐｒｅｆｌｕｓｈ ａｃｉｄ ｆｒａｃ— ｔｕｒｉｎｇ ｃｏｎｓｉｄｅｒｉｎｇ ｒｏｃｋ ｈｅｔｅｒｏｇｅｎｅｉｔｙ［Ｊ】．Ｎａｔｕｒａｌ Ｇａｓ Ｉｎｄｕｓｔｒｙ， ２０ １ ８，３８（２）：５９－６６． ［２】赵立强，邹洪岚，刘平礼．对酸液有效作用距离预测方法的改 进［Ｊ】．天然气工业，１９９８，１８（５）：５９．６３． Ｚｈａｏ Ｌｉｑｉａｎｇ．Ｚｏｕ Ｈｏｎｇｌａｎ＆Ｌｉｕ Ｐｉｎｇｌｉ．Ａｎ ｉｍｐｒｏｖｅｍｅｎｔ ｏｎ ｔｈｅ ｍｅｔｈｏｄ ｆｏｒ ｆｏｒｅｃａｓｔｉｎｇ ｔｈｅ ｅｆｆｅｃｔｉｖｅ ａｃｔｉｏｎ ｄｉｓｔａｎｃｅ ｏｆ ａｃｉｄｉｚ－ ｉｎｇ ｌｆｕｉｄ［Ｊ］．Ｎａｔｕｒａｌ Ｇａｓ Ｉｎｄｕｓｔｒｙ，１９９８，１８（５）：５９—６３． ［３］Ｌｙｏｎｓ Ｊ，Ｎａｓｒａｂａｄｉ Ｈ＆Ｎａｓｒ—Ｅ１·Ｄｉｎ ＨＡ．Ａ ｎｏｖｅｌ ｐｏｒｅ—ｓｃａｌｅ ｔｈｅｒｍａｌ—ｆｒａｃｔｕｒｅ—ａｃｉｄｉｚｉｎｇ ｍｏｄｅｌ ｗｉｔｈ ｈｅｔｅｒｏｇｅｎｅｏｕｓ ｒｏｃｋ ｐｒｏｐ－ ｅｒｔｉｅｓ［Ｊ］．ＳＰＥ Ｊｏｕｒｎａｌ，２０１６，２１（１）：２８０—２９２． ［４】王荣．Ｄ油田灰岩储层酸化技术研究［Ｄ】．西南石油大学， ２Ｏ１４． Ｗａｎｇ Ｒｏｎｇ．Ｓｔｕｄｙ ｏｎ ａｃｉｄｉｉｆｃａｔｉｏｎ ｔｅｃｈｎｏｌｏｇｙ ｏｆ ｌｉｍｅｓｔｏｎｅ ｒｅｓ— ｅｒｖｏｉｒ ｉｎ Ｄ Ｏｉｌｉｆｅｌｄ［Ｄ］．Ｓｏｕｔｈｗｅｓｔ Ｐｅｔｒｏｌｅｕｍ Ｕｎｉｖｅｒｓｉｔｙ，２０１４． ［５】Ｙｏｓｈｉｍｕｒａ Ｋ，Ｍａｔｓｕｉ Ｈ＆Ｍｏｒｉｔａ Ｎ．Ｄｅｖｅｌｏｐｍｅｎｔ ｏｆ ｐｏｌｙｇｌｙ－ ｃｏｌｉｃ ａｎｄ ｐｏｌｙａｃｔｉｃ ａｃｉｄ ｆｌｕｉｄ—ｌｏｓｓ ｃｏｎｔｒｏｌ ｍａｔｅｒｉａｌｓ ｆｏｒ ｆｒａｃｔｕｒｉｎｇ ｌｆｕｉｄｓ［Ｊ］．ＳＰＥ Ｄｒｉｌｌｉｎｇ＆Ｃｏｍｐｌｅｔｉｏｎ，２０１６，３０（４）：２９５—３０９． ［６】ＺｈａｎｇＹｅ，Ｙａｎｇ Ｓｈｅｎｇｌａｉ，Ｚｈａｎｇ Ｓｈｉｃｈｅｎｇ＆Ｍｏｕ Ｊｉａｎｙｅ．Ｗｏｒｍ— ｈｏｌｅ ｐｒｏｐａｇａｔｉｏｎ ｂｅｈａｖｉｏｒ ａｎｄ ｉｔｓ ｅｆｆｅｃｔ ｏｎ ａｃｉｄ ｌｅａｋｏｆｆ ｕｎｄｅｒ ｉｎ ｓｉｕｔ ｃｏｎｄｉｔｉｏｎｓ ｉｎ ａｃｉｄ ｆｒａｃｔｕｒｉｎｇ［Ｊ］．Ｔｒａｎｓｐｏｒｔ ｉｎ Ｐｏｒｏｕｓ Ｍｅｄｉａ， ２０１４，１０１（１）：９９—１ １４． ［７］Ｋａｌｉａ Ｎ＆Ｂａｌａｋｏｔａｉａｈ Ｖ．Ｍｏｄｅｌｉｎｇ ａｎｄ ａｎａｌｙｓｉｓ ｏｆ ｗｏｒｍｈｏｌｅ ｆｏｒｍａｔｉｏｎ ｉｎ ｒｅａｃｔｉｖｅ ｄｉｓｓｏｌｕｔｉｏｎ ｏｆ ｃａｒｂｏｎａｔｅ ｒｏｃｋｓ［Ｊ］．Ｃｈｅｍｉｃａｌ Ｅｎｇｉｎｅｅｒｉｎｇ Ｓｃｉｅｎｃｅ，２００７，６２（４）：９１９—９２８． 【８］Ｒｏｂｅｒｔｓ ＬＤ＆Ｇｕｉｎ ＪＡ．Ａ ｎｅｗ ｍｅｔｈｏｄ ｆｏｒ ｐｒｅｄｉｃｔｉｎｇ ａｃｉｄ ｐｅｎｅ— ｔｒａｔｉｏｎ ｄｉｓｔａｎｃｅ［Ｊ］．ＳＰＥ Ｊｏｕｒｎａｌ，１９７５，１５（４）：２７７－２８６． ［９］Ｒｏｍｅｒｏ Ｊ，Ｇｕ Ｈ＆Ｇｕｌｒａｊｎａｉ ＳＮ．Ｔｈｒｅｅ—ｄｉｍｅｎｓｉｏｎａｌ ｔｒｎａｓｐｏｒｔ ｉｎ ａｃｉｄ ｆｒａｃｔｕｒｉｎｇ ｔｒｅａｔｍｅｎｔｓ：Ｔｈｅｏｒｅｔｉｃａｌ ｄｅｖｅｌｏｐｍｅｎｔ ａｎｄ ｃｏｎｓｅ— ｑｕｅｎｃｅｓ ｆｏｒ ｈｙｄｒｏｃａｒｂｏｎ ｐｒｏｄｕｃｔｉｏｎＩｓ］．ＳＰＥ Ｐｒｏｄｕｃｔｉｏｎ＆Ｆａｃｉｌ— ｉｔｉｅｓ，２００１，１６（２）：１２２—１３０， ［１０］任书泉，熊宏杰．同离子效应和温度场对裂缝中酸液有效作 用距离的影响［Ｊ］．西南石油大学学报（自然科学版），１９８６， ８（４）：２３－３５． Ｒｅｎ Ｓｈｕｑｕａｎ，Ｘｉｏｎｇ Ｈｏｎｇｊｉｅ．Ｅｆｆｅｃｔｓ ｏｆ ｔｅｍｐｅｒａｔｕｒｅ ｉｆｅｌｄ ａｎｄ ｃｏｍｍｏｎ ｉｏｎ ｅｆｆｅｃｔ ｏｎｅｆｆｅｃｔｉｖｅ ａｃｉｄ ｐｅｎｅｔｒａｔｉｏｎ ｉｎ ｆｒａｃｔｕｒｅ ｓｙｓ— ｔｅｍｓ［Ｊ］．Ｊｏｕｒｎａｌ ｏｆ Ｓｏｕｔｈｗｅｓｔ Ｐｅｔｒｏｌｅｕｍ Ｉｎｓｔｉｔｕｔｅ（Ｅｄｉｔｉｏｎ ｏｆ Ｎａｔｕｒａｌ Ｓｃｉｅｎｃｅ），１９８６，８（４）：２３－３５． ［１１】赵碧华，任书泉．同离子效应对酸液有效作用距离的影响［Ｊ］． 石油学报，１９８６，７（１１：７４．８５． Ｚｈａｏ Ｂｉｈｕａ＆Ｒａｎ Ｓｈｕｑｕａｎ．Ｔｈｅ ｉｎｆｌｕｅｎｃｅ ｏｆ ｃｏｍｍｏｎ ｉｏｎ ｅｆｆｅｃｔ ｏｎｅｆｆｅｃｔｉｖｅ ａｃｉｄ ｐｅｎｅｔｒａｔｉｏｎ［Ｊ］．Ａｃｔａ Ｐｅｔｒｏｌｅｉ Ｓｉｎｉｃａ，１９８６，７（１）： ７４—８５． ［１２］赵碧华，任书泉．考虑同离子效应酸液有效作用距离计算图 版的研制［Ｊ】．石油钻采工艺，１９８５，７（３）：４５．５３． Ｚｈａｏ Ｂｉｈｕａ＆Ｒｅｎ Ｓｈｕｑｕａｎ．Ｄｅｖｅｌｏｐｍｅｎｔ ｏｆ ｃａｌｃｕｌａｔｉｎｇ ｐｌａｔｅ ｏｆ ｅｆｆｅｃｔｉｖｅ ｄｉｓｔａｎｃｅ ｏｆ ａｃｉｄ ｅｆｆｅｃｔ［Ｊ］．Ｏｉｌ Ｄｒｉｌｌｉｎｇ＆Ｐｒｏｄｕｃｔｉｏｎ Ｔｅｃｈｎｏｌｏｇｙ，１ ９８５，７（３）：４５－５３． ［１３】李平．前置液酸压缝中酸液浓度分布的全隐式数值计算方法 ［Ｊ］．天然气工业，１９８８，８（４）：５６—６３． Ｌｉ Ｐｉｎｇ．Ｆｕｌｌｙ ｉｍｐｌｉｃｉｔ ｎｕｍｅｒｉｃａｌ ｃａｌｃｕｌａｔｉｏｎ ｍｅｔｈｏｄ ｆｏｒ ａｃｉｄ ｓｏｌｕｔｉｏｎ ｃｏｎｃｅｎｔｒａｔｉｏｎ ｄｉｓｔｒｉｂｕｔｉｏｎ ｉｎ ａｃｉｄｉｚｅｄ ａｎｄ ｆｒａｃｔｕｒｅｄ ｆｉｓ． ｓｕｒｅｓ ｉｎ ｃａｓｅ ｏｆ ｕｓｉｎｇ ｐｒｅｐａｄ［Ｊ］．Ｎａｔｕｒａｌ Ｇａｓ Ｉｎｄｕｓｔｒｙ，１ ９８８，８（４） ５６—６３． ［１４】任书泉，李联奎．酸在裂缝中流动有效作用距离的计算［Ｊ］．石 油勘探与开发，１９８１，３（６）：６５．７５． Ｒｅｎ Ｓｈｕｑｕａｎ＆Ｌｉ Ｌｉａｎｋｕｉ．Ｃａｌｃｕｌａｔｉｏｎ ｏｆ ｅｆｆｅｃｔｉｖｅ ｄｉｓｔａｎｃｅ ｏｆ ａｃｔｉｏｎ ｏｆ ａｃｉｄ ｉｎ ｃｒａｃｋ［Ｊ］．Ｐｅｔｒｏｌｅｕｍ Ｅｘｐｌｏｒａｔｉｏｎ ａｎｄ Ｄｅｖｅｌｏｐ— ｍｅｎｔ，１９８１，３（６）：６５—７５． ［１５］侯博恒，施松杉．碳酸盐岩酸化压裂ＣＯ：量的计算［Ｊ】＿断块 油气田，２０１３，２Ｏ（５）：６５６．６５８． Ｈｏｕ Ｂｏｈｅｎｇ＆Ｓｈｉ Ｓｏｎｇｓｈａｎ．Ｃａｌｃｕｌａｔｉｏｎ ｏｆ ｃａｒｂｏｎ ｄｉｏｘｉｄｅ ｖｏ１． ｕｍｅ ｉｎ ａｃｉｄ ｒｆａｃｔｕｒｉｎｇ ｏｆ ｃａｒｂｏｎａｔｅ ｒｏｃｋ［Ｊ］．Ｆａｕｌｔ—Ｂｌｏｃｋ Ｏｉｌ ａｎｄ Ｇａｓ Ｆｉｅｌｄ，２０１３，２Ｏ（５）：６５６—６５８． ［１６］李勇明，郭建春，赵金洲，吴小庆，李娅．裂缝性气藏压裂 液滤失模型的研究及应用［Ｊ］．石油勘探与开发，２００４，３ｌ（５）： １２Ｏ．１２２． Ｌｉ Ｙｏｎｇｍｉｎｇ，Ｇｕｏ Ｊｉａｎｃｈｕｎ，Ｚｈａｏ Ｊｉｎｚｈｏｕ，Ｗｕ Ｘｉａｏｑｉｎｇ＆Ｌｉ Ｙａ．Ｎｅｗ ｍｏｄｅｌ ｆｏｒ ｆｒａｃｔｕｒｉｎｇ ｌｆｕｉｄ ｌｅａｋｏｆｆ ｉｎ ｎａｔｕｒａｌｌｙ ｆｒａｃｔｕｒｅｄ ｇａｓ ｉｆｅｌｄｓ ａｎｄ ｉｔｓ ａｐｐｌｉｃａｔｉｏｎ［Ｊ］．Ｐｅｔｒｏｌｅｕｍ Ｅｘｐｌｏｒａｔｉｏｎ ａｎｄ Ｄｅ－ ７２ 天然气工业 ２０１８年７月 ｖｅｌｏｐｍｅｎｔ，２００４，３１（５）：１２０—１２２． ｂｌｏｃｋ ｃｏｐｏｌｙｍｅｒ［Ｊ］．Ｊｏｕｒｎａｌ ｏｆ Ｐｏｒｏｕｓ Ｍａｔｅｒｉａｌｓ，２０１６，２３（６）： １４３１—１４３８． 【１７】Ｌｉｕ Ｐｉｙａｎｇ，Ｙａｏ Ｊｕｎ，Ｃｏｕｐｌｅｓ ＧＤ，Ｈｕａｎｇ Ｚｈａｏｑｉｎ，Ｓｕｎ Ｈａｉ＆ Ｍａ Ｊｉｎｇｓｈｅｎｇ．Ｎｕｍｅｒｉｃａｌ ｍｏｄｅｌｌｉｎｇ ａｎｄ ａｎａｌｙｓｉｓ ｏｆ ｒｅａｃｔｉｖｅ ｌｆｏｗ ［２１］王强．天然裂缝性碳酸盐岩储层酸化模拟研究［Ｄ】．成都：西 南石油大学，２００６． Ｗａｎｇ Ｑｉａｎｇ．Ｓｔｕｄｙ ｏｎ ｔｈｅ ａｃｉｄｉｉｆｃａｔｉｏｎ ｏｆ ｎａｔｒａｌｕ ｆｒａｃｔｕｒｅｄ Ｃａｒ－ ａｎｄ ｗｏｒｍｈｏｌｅ ｆｏｒｍａｔｉｏｎ ｉｎ ｆｒａｃｔｕｒｅｄ ｃａｒｂｏｎａｔｅ ｒｏｃｋｓ［Ｊ］．Ｃｈｅｍｉ— ｃａｌ Ｅｎｇｉｎｅｅｒｉｎｇ Ｓｃｉｅｎｃｅ，２０１７，１７２：１４３—１５７． 【１ ８】Ｊａｈａｎｉ Ｎ，Ｈａｕｇｅｎ Ｂ＆Ｂｅｒｇｅ Ｇ．Ｃｏｕｐｌｅｄ ｆｌｕｉｄ ｆｌｏｗ ａｎｄ ｅｌａｓ— ｔｏ－·ｐｌａｓｔｉｃ ｄａｍａｇｅ ａｎａｌｙｓｉｓ ｆｏｒ ｆｒａｃｔｕｒｅｄ ｐｏｒｏｕｓ ｃｈａｌｋ ｗｉｔｈ ｉｎ— ｂｏｎａｔｅ ｒｅｓｅｒｖｏｉｒ［Ｄ］．Ｃｈｅｎｇｄｕ：Ｓｏｕｔｈｗｅｓｔ Ｐｅｔｒｏｌｅｕｍ Ｕｎｉｖｅｒｓｉｔｙ， ２００６． ｄｕｃｅｄ ｗｏｒｍｈｏｌｅ［Ｊ］．Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｏｆ Ｒｏｃｋ Ｍｅｃｈａｎｉｃｓ＆ Ｍｉｎｉｎｇ Ｓｃｉｅｎｃｅｓ，２０１５，８０（１）：１２９—１３６． 【２２］Ｍｏｕ Ｊｉｎｙｅ，Ｚｈａａｎｇ Ｓｈｉｃｈｅｎｇ＆Ｚｈａｎｇ Ｙｅ．Ａｃｉｄ ｌｅａｋｏｆｆ ｍｅｃｈａ— ｎｉｓｍ ｉｎ ａｃｉｄ ｆｒａｃｔｕｒｉｎｇ ｏｆ ｎａｔｕｒａｌｌｙ ｆｒａｃｔｕｒｅｄ ｃａｒｂｏｎａｔｅ Ｏ－１ ｒｅｓｅ￣ 【１９］Ｋａｌｉａ Ｎ＆Ｂａｌａｋｏｔａｉａｈ Ｖ Ｅｆｆｅｃｔ ｏｆ ｍｅｄｉｕｍ ｈｅｔｅｒｏｇｅｎｅｍｅｓ ｏｎ ｒｅａｃｔｉｖｅ ｄｉｓｓｏｌｕｔｉｏｎ ｏｆ ｃａｒｂｏｎａｔｅｓ［Ｊ］．Ｃｈｅｍｉｃａｌ Ｅｎｇｉｎｅｅｒｉｎｇ Ｓｃｉ— ｅｎｃｅ，２００９，６４（２）：３７６—３９０． ｖｏｉｒｓ［Ｊ］．Ｔｒａｎｓｐｏｎ ｉｎ Ｐｏｒｏｕｓ Ｍｅｄｉａ，２０１２，９１（２）：５７３—５８４． ［２３］Ｓｃｈｅｉｄｅｇｇｅｒ ＡＥ．Ｔｈｅ ｐｈｙｓｉｃｓ ｏｆ ｌｆｏｗ ｔｈｒｏｕｇｈ ｐｏｒｏｕｓ ｍｅｄｉａ［Ｊ］． Ｓｏｉｌ Ｓｃｉｅｎｃｅ，１９５８，８６（６）：３５５． ［２０】Ｊｉａｏ Ｊｉａｎ，Ｘｉａｎｇ Ｙｉｑｉｎ，Ｃａｏ Ｊｉｎｇ＆Ｘｉａ Ｙ Ｗｏｒｍ－ｈｏｌｅ ｓｔｒｕｃｔｕｒｅｄ ｍｅｓｏｐｏｒｏｕｓ ｃａｒｂｏｎ ｍｏｎｏｌｉｔｈｓ ｓｙｎｔｈｅｓｉｚｅｄ ｗｉｔｈ ａｍｐｈｉｐｈｉｌｉｃ ｔｒｉ— （修改回稿日期２０１８—０６．０５编辑孔玲） 国内首座橇装化天然气处理厂建成投运 ２０１８年７月６日，随着长庆神木气田佳县天然气处理厂放空火炬的点燃，标志着国内首座橇装一体化天然气处理厂在长庆气田建成 并顺利投运。一座年处理天然气５×１０ ｍ ，日处理能力达到１５０×１０ ｍ 的大型联合一体化天然气处理厂屹立在长庆气田的最东部。 长庆神木气田佳县天然气处理厂是由中国石油天然气股份有限公司（以下简称股份公司）选树的重点工程，该处理厂及其配套工程 由西安长庆科技工程有限责任公司采用ＥＰＣ（工程总承包）管理模式承建，并成立佳县天然气处理站工程ＥＰＣ项目部，主要对设计、采购、 施工过程进行统一协调组织和管理。此项工程建设工期从２０１７年７月１日动工，有效工期缩短为７个月建成投产使用。 目前，橇装技术成为采油、采气厂站技术升级的主攻方向。佳县天然气处理厂工程设计建设主要采用先丙烷制冷脱油脱水后增压工 艺主体技术，全面推行标准化施工技术，严格按照股份公司标准化施工要求组织建设，确保将工程建成示范工程、标杆工程。佳县天然 气处理厂工程ＥＰＣ项目部坚持“一业为主、两翼齐飞、创新驱动”的发展战略，以做精ＥＰＣ工程总承包工程为目标，以项目全过程管 理为核心，以标准化施工为主线，以质量、安全监管为重点，提升ＥＰＣ项目管理水平，在气田项目建设中，形成了全过程、全要素、全 流程的项目管理新模式。 国内首座橇装一体化天然气处理厂在长庆气田的建成，在油气行业中具有示范引领作用，应用模块化集成缩短建设周期、节约用地 近５０亩（１亩＝６６６．６７ ｍ ），并能有效减少对当地环境的污染，节能减排；同时设备可重复利用，并最终实现无人值守（可减少劳动用 工人数２３人），降低了人员和运营的成本费用。 西安长庆科技工程有限责任公司ＥＰＣ项目部主要负责．人介绍说：“佳县天然气处理站的建成，实现了气源互补、灵活调配的枢纽转 输点功能。佳县天然气处理站的天然气经榆林天然气处理厂配气站转输至陕京天然气管线后，有效缓解了首都北京及周边地区天然气供 应紧张的压力。” 据悉，长庆神木气田今年将建成达到１９×１０ ｍ 天然气的处理能力，届时将更好地缓解神木气田周边输气管网的运行压力，为全国 冬供保气增添新的动力。 （天工摘编自天然气工业网）
