二元复合型泡排剂COG的研制及现场试验

王维恒; 陆俊华; 韩倩

doi:10.11911/syztjs.2022053

二元复合型泡排剂COG的研制及现场试验

中国石化华东油气分公司石油工程技术研究院，江苏南京 210004

详细信息

作者简介:
王维恒（1989—），男，江苏南京人，2011年毕业于中国地质大学（北京）材料化学专业，2014年获中国地质大学（北京）应用化学专业工学硕士学位，助理研究员，从事油田化学助剂开发工作。E-mail：beiton@163.com

中图分类号: TE357.46⁺2
计量
- 文章访问数: 273
- HTML全文浏览量: 181
- PDF下载量: 35
出版历程
- 收稿日期: 2021-06-14
- 修回日期: 2022-04-18
- 网络出版日期: 2022-05-04
- 刊出日期: 2022-06-08

Preparation and Field Test of Binary Complex Foam Agent COG

Petroleum Engineering Technology Research Institute, Sinopec East China Oil & Gas Company, Nanjing, Jiangsu, 210004, China

摘要

摘要:
针对平桥南区块高温、高矿化度井所用泡排剂效果不佳的问题，研制了适用于该区块的泡排剂。优选椰油酰胺丙基甜菜碱（CAB）、十二烷基二甲基氧化铵（OB）和羟丙基胍胶（Guar）作为原料，复配制得了二元复合型泡排剂COG；采用正交试验方法，分析了CAB、OB和Guar质量比对COG起泡性能的影响，并确定了最佳加量。室内试验评价了泡排剂COG的发泡性能、稳泡性能及抗温抗盐性能，在温度95 ℃、矿化度10×10⁴ mg/L条件下，其综合性能明显优于现场常用的2种泡排剂。在平桥南区块3口井进行了现场试验，试验结果表明，泡排剂COG具有较好的排水增产效果，用后单井产气量比应用原有泡排剂提高了10%以上。研究结果表明，二元复合型泡排剂COG适用于平桥南区块页岩气井的泡沫排水，具有较好的推广应用价值。
- 页岩气井 /
- 泡排剂 /
- 研制 /
- 表面活性剂 /
- 平桥南区块
Abstract:
The foaming agents applied to the high-temperature and high-salinity wells in South Pingqiao Block prove to be ineffective, a foaming agent readily applicable to this block was thus prepared to address this problem. Cocamidopropyl betaine (CAB), dodecyldimethylamine oxide (OB), and hydroxypropyl guanidine gum (Guar) were adopted as raw materials to prepare a binary composite foaming agent COG. The influences of the mass ratios of CAB, OB, and Guar on the foaming performance of COG were analyzed by the orthogonal test method, and their optimal dosages was determined. The foaming properties, foaming stability, temperature resistance, and salt resistance of the proposed foaming agent COG were evaluated in the laboratory, and its comprehensive performance was significantly better than that of the two foaming agents commonly used in the field at the temperature of 95 ℃ and the salinity of 10 × 10⁴ mg/L. The test results of three wells in South Pingqiao Block showed that the foaming agent COG provided a favorable effect of drainage production stimulation, with the single-well gas production more than 10% higher than that obtained with the previous foaming agents. In research results show that the proposed binary composite foaming agent COG is suitable for the foam drainage of the shale gas wells in South Pingqiao Block, and it has high popularization and application value.
- shale gas well /
- foaming agent /
- preparation /
- surfactant /
- South Pingqiao Block

HTML全文

具有低渗透、非均质和各向异性特征的复杂油气藏逐渐成为油气勘探开发的重点，无论从储层评价还是钻井工程来看，对近钻头地质导向系统均有迫切需求^[1–2]。目前，近钻头伽马成像测井已成为国外随钻地质导向作业的常用技术，在钻井过程中首先利用安装在钻铤侧面的单个或多个伽马探测器、定向探测地层中某一扇区内总的自然伽马射线强度，然后随着井下钻具旋转动态扫描获取全井眼不同方位扇区地层的伽马成像图，最后根据成像信息来估算地层与井眼的变化趋势，据此调整钻头前进方向，以实现地质导向作业^[3–5]。D. C. Minette^[6]最早介绍了方位伽马测井，在LWD仪器随钻头旋转过程中，将测量的伽马强度分组到相应的各个扇区，得到各个扇区的伽马测量值；后来，R. L. Spross^[7]根据仪器与井壁的环空间隙对钻进中各个方位的采样数据进行了加权，得到了各个扇区的伽马计数，实现了扇区伽马成像；M. Bittar等人^[8]在扇区伽马成像测井中引入了“成像深度”的概念，以精确计算地层倾角，并准确进行地质导向作业；K. McKinny等人^[9]介绍了正弦曲线自动拾取方法，并采用蒙特卡罗模拟和相关数学方法，进行了成像计算倾角的误差分析；袁超等人^[10]进行了扇区伽马成像正演模拟，分析了井眼条件对方位伽马成像的影响；Wang Jiaxin等人^[11]采用平均值拟合法实现了正弦曲线的自动拾取，利用蒙特卡罗模拟分析了地层密度、井筒环境、钻井液密度和仪器偏心等因素对成像深度的影响，并给出了计算倾角的方法。

地质导向作业中，倾斜地层界面在近钻头伽马成像图中呈现出正弦曲线特征，提取该曲线的特征参数，可以分析地层走向和井眼轨迹。但是，现有随钻测井数据遥测系统的传输能力有限，无法将大量的扇区伽马成像数据实时上传至地面测控系统，因而不能及时地指导司钻进行井眼轨迹调整，给钻井作业带来了风险。近年来，大规模集成器件（如数字信号处理器DSP）在随钻测井仪器中得到了广泛应用，可以实现对井下复杂情况的控制，并且具有较强的计算处理能力^[12–13]。其中，DSP中实时多扇区伽马测量值拟合是近钻头伽马成像测井仪器的关键技术之一。为了满足现场实时获取伽马成像图的迫切需求，笔者提出了一种基于扇区伽马成像的快速正弦曲线拟合方法，并将其应用于近钻头伽马成像测井中井下固件算法的设计，以解决成像数据与传输速率不匹配的问题；在此基础上，试验研究了用标准岩样构造分层倾斜地层模拟井眼采集的伽马数据，反演得到8扇区伽马成像图，并验证了快速正弦曲线拟合方法的可行性和准确性。

1. 倾斜界面近钻头伽马成像特征

为分析倾斜界面近钻头伽马成像特征，建立了充满液体井眼倾斜界面地层模型，如图1所示。模型中，地层界面与井眼轴线的夹角为α，井眼直径R为215.9 mm，井眼内充满淡水；钻头附近，钻铤外径r_o为172.0 mm，钻铤侧面开槽放置1个伽马探测器，钻铤水眼中充满淡水，水眼直径r_i为57.0 mm；上部地层为低放射性的纯砂岩目的层，放射性强度为0；下部地层为高放射性的泥岩地层。模拟计算时，使近钻头伽马成像测井仪器旋转，从上部低放射性砂岩层穿过下部高放射性泥岩，记录伽马探测器在8或16个方位的伽马计数率，并采用样条插值方法获取全井眼8扇区或16扇区伽马成像图。

图 1 充满液体井眼倾斜界面地层模型

Figure 1. Inclined interface stratigraphic model with liquid-filled wellbore

下载: 全尺寸图片幻灯片

近钻头伽马成像测井仪器从上部低放射性目的层穿过倾斜界面进入高放射性泥岩层，利用蒙特卡罗方法模拟获取了8扇区伽马成像图（见图2）。

图 2 倾斜界面近钻头8扇区伽马成像图

Figure 2. Inclined interface near-bit 8-sector gamma image

下载: 全尺寸图片幻灯片

从图2可以看出，当近钻头伽马成像测井仪器随着钻进方向测量时，伽马探测器在井眼内逐渐从顶部靠近地层边界，在方位角0°处首先探测到高放射性泥岩层；当该仪器随着钻进继续前移，其他方位开始探测到高放射性泥岩层；最后，该仪器完全通过层界面，并全部进入到高放射性泥岩层中；近钻头伽马成像测井仪器穿过倾斜界面进入高放射性泥岩层时，地层界面在扇区伽马成像图中呈现正弦曲线特征。成像图可直观地显示高放射性地层，亮度越大，表示该位置的放射性强度越高。

根据扇区伽马成像图的正弦曲线特征，可计算出地层倾斜界面与井眼相交的相对倾角，计算公式为：

$\alpha \;{\rm{ = }}\;\arctan {\frac{H}{{D + 2d}}}$

(1)

式中：H为伽马成像图正弦波高度（地层边界正弦曲线的波峰与波谷的垂直高度），mm；D为井眼直径，mm；d为在垂直井眼轴线方向上观测到伽马成像的地层深度（也称伽马成像深度），mm。

H和d可以从伽马成像图中得到。但由于H和d受井眼和地层条件的影响，计算得到的相对倾角也会受到井眼和地层条件的影响。而近钻头伽马成像测井仪器以一定井眼相对倾角穿过地层界面进入高放射性泥质围岩层，可以测量得到呈正弦曲线特征的扇区成像图。利用扇区伽马成像图可以精确反演井眼相对倾角和围岩地层厚度，从而确定地层边界，实现井眼轨迹精确控制。

2. 快速正弦曲线拟合方法

基于扇区伽马成像图，通过拟合正弦曲线可以提取成像图中特征参数，利用式（1）可以计算井眼相对倾角；与之相似，扇区伽马成像测量值与方位角的关系也可采用正弦曲线来表征。前人已对正弦曲线的拟合方法进行了大量研究，例如四参数正弦波拟合算法属于非线性迭代拟合，若拟合初始值距离目标值“太远”，则容易导致迭代过程发散或收敛到局部最优点而不是总体最优点，致使拟合结果错误^[14]；并且，利用四参数正弦波拟合算法进行拟合耗时较长，会影响算法的时效性和适用性。另外，近钻头伽马成像测井仪器井下电子系统的处理能力有限，无法在井下进行大量的数学迭代运算。为此，针对实时获取井下伽马成像图的需求以及考虑随钻测井数据传输速率的限制，基于扇区伽马成像图呈现正弦曲线的特征，提出了一种快速正弦曲线拟合方法：先对扇区成像图中正弦曲线的频率进行估计，再在频率已知的条件下进行三参数正弦曲线快速拟合，以获得正弦曲线的幅度、相位和直流分量等3个参数，最后利用正弦曲线的幅度、频率、相位和直流分量等4个参数来拟合反演形成扇区伽马成像图。

多扇区（8或16扇区）伽马成像图中某一正弦曲线可以看成均匀采样后获得的等间隔采样序列（如图3所示），记为 $G{a_1},\;G{a_2},\; \cdots ,\;G{a_n}$ 。设方位采样间隔为 $\Delta Az$ ，方位角采样数据序列表示为 $A{z_i} = (i - 1)\Delta Az,$ $\;i = 1,\;2,\; \cdots ,\;n$ ，数字角频率表示为 $\omega = 2{\text{π}}f\Delta Az$ ，则理想正弦信号的离散形式为：

图 3 多扇区伽马成像数据序列

Figure 3. Multi-sector gamma imaging data sequence

下载: 全尺寸图片幻灯片

$\begin{array}{l} Ga(i) = E\cos [\omega (i - 1) + \varphi ] + Q \end{array}$

(2)

式中： $E$ 为信号幅值，API； $\varphi$ 为信号的初相位，rad； $Q$ 为信号的直流分量，API。

令 $x(i) = E\cos [\omega (i - 1) + \varphi ],\;i = 1,\;2,\; \cdots ,\;n$ ，记伽马测量值 $G{a_i}$ 的误差为 ${\gamma _i}$ ，并令 $g = 2\cos \omega$ ，则：

$G{a_i} = Ga(i) + {\gamma _i} = x(i) + Q + {\gamma _i}$

(3)

$x(i) + x(i - 2) = 2\cos \omega E\cos [\omega (i - 2) + \varphi ] = gx(i - 1)$

(4)

$Ga(i) - Q + Ga(i - 2) - Q = g [Ga(i - 1) - Q]$

(5)

$G{a_i} - Q - {\gamma _i} + G{a_{i - 2}} - Q - {\gamma _{i - 2}} = g (G{a_{i - 1}} - Q - {\gamma _{i - 1}})$

(6)

$G{a_{i + 1}} - Q - {\gamma _{i + 1}} + G{a_{i - 1}} - Q - {\gamma _{i - 1}} = g(G{a_i} - Q - {\gamma _i})$

(7)

再令 ${z_i} = G{a_i} - G{a_{i - 1}}$ 、 ${\zeta _i} = {\gamma _i} - {\gamma _{i - 1}}$ ，式（6）和式（7）相减，得：

${z_{i + 1}} + {z_{i - 1}} - g{z_i} = {\zeta _{i + 1}} + {\zeta _{i - 1}} - g{\zeta _i}$

(8)

由于伽马测量值的误差 ${\gamma _i}$ 为随机误差，故式（8）右侧同样可以认为是随机误差。通过选取 $g$ 的取值，使 $\,\rho = \sum\limits_{i\;{\rm{ = }}\;3}^{n - 1} {{{({z_{i + 1}} + {z_{i - 1}} - g{z_i})}^2}}$ 最小，即 ${{{\rm{d}}\rho } / {{\rm{d}}g}} = 0$ ，则：

$g = \frac{{{z_n}{z_{n - 1}} + {z_3}{z_2} + 2\sum\limits_{i = 4}^{n - 1} {{z_i}{z_{i - 1}}} }}{{\sum\limits_{i = 3}^{n - 1} {z_i^2} }}$

(9)

$\tilde \omega = \arccos \frac{g }{2}$

(10)

式中： $\tilde \omega$ 为数字角频率的最小二乘估计值。

基于最小二乘估计获得正弦曲线的频率 $f$ 后，利用三参数正弦拟合算法构造残差平方和 $\varepsilon$ ，即：

$\varepsilon = \sum\limits_{i = 1}^n {{{\left[ {G{a_i} - {A_1}\cos (2{\text{π}}fA{z_i}) - {B_1}\sin (2{\text{π}}fA{z_i}) - C} \right]}^2}}$

(11)

将式（11）改写成矩阵形式：

$\varepsilon = {({{y}} - {{{D}}_0}{{{s}}_0})^{\rm{T}}}({{y}} - {{{D}}_0}{{{s}}_0})$

(12)

$\!{\text{其中}}\quad{{{D}}_0} = \left[ {\begin{array}{*{20}{c}} {\cos (2{\text{π}}fA{z_1})}&{\sin (2{\text{π}}fA{z_1})}&1\\ {\cos (2{\text{π}}fA{z_2})}&{\sin (2{\text{π}}fA{z_2})}&1\\ \vdots & \vdots & \vdots \\ {\cos (2{\text{π}}fA{z_n})}&{\sin (2{\text{π}}fA{z_n})}&1 \end{array}} \right] \!\!\!\!\!\!$

(13)

${{y}} = \left[ {\begin{array}{*{20}{c}} {G{a_1}}\\ {G{a_2}}\\ \vdots \\ {G{a_n}} \end{array}} \right]$

(14)

${{{s}}_0} = \left[ {\begin{array}{*{20}{c}} {{A_1}}\\ {{B_1}}\\ C \end{array}} \right]$

(15)

根据最小二乘拟合方法，使式（12）所述残差平方和 $\varepsilon$ 最小^[15]，则式（12）的最小二乘解为：

${\hat {{s}}_0} = {({{D}}_0^{\rm{T}}{{{D}}_0})^{ - 1}}({{D}}_0^{\rm{T}}{{y}})$

(16)

故正弦曲线的拟合函数为：

$\tilde Ga(i) = {A_1}\;\cos (\tilde \omega i) + {B_1}\;\sin (\tilde \omega i) + C$

(17)

其幅度和相位表达式为：

$\tilde Ga(i) = A\cos (2{\text{π}}f{x_i} + \theta ) + C$

(18)

$\!{\text{其中}} \qquad\qquad\qquad A = \sqrt {A_1^2 + B_1^2}\quad\quad$

(19)

$\theta = \left\{ \begin{array}{l} \arctan \dfrac{{ - {B_1}}}{{{A_1}}},\;\;\;A_1 \geqslant \;0 \\ \arctan \dfrac{{ - {B_1}}}{{{A_1}}} + {\text{π}},\;A_1\; < \;0 \\ \end{array} \right.$

(20)

拟合残差有效值为：

$\rho = \sqrt {\frac{\varepsilon }{n}}$

(21)

上面所推导的是信号频率估计算法和三参数正弦拟合算法的相关公式。将信号频率估计算法和三参数正弦拟合算法进行结合，即为快速正弦曲线拟合方法。利用快速正弦曲线拟合方法获取正弦曲线的幅度、频率、相位和直流分量等4个参数，再利用式（18），就可拟合反演形成扇区伽马成像图。相比于四参数拟合方法，上述拟合过程只涉及四则运算，不需要迭代运算，属于闭合的线性过程，绝对收敛。快速正弦曲线拟合方法具有速度快、过程简洁的特点，适合用于提取井下伽马成像图正弦曲线特征参数。

3. 快速正弦曲线拟合固件算法的实现

近钻头伽马成像测井仪器下短节采用DSP进行各种计算和处理，DSP的结构和指令集可高效地实现乘积、累加和傅里叶变换等功能，数据处理主要包括滤波、伽马计数累加、多扇区伽马测量值拟合、数据压缩和可变增益控制等算法。井下DSP实时多扇区伽马测量值拟合是近钻头伽马成像测井仪器的关键技术之一，为此，设计了近钻头伽马成像测井的快速正弦曲线拟合固件算法，以拟合多扇区伽马成像测量值。快速正弦曲线拟合方法具有绝对收敛、算法简洁、运算速度快和内存占用空间小等特点，可以较好地实现多扇区伽马成像测量值的拟合。

仪器下短节DSP控制及数据处理流程如图4所示。仪器处于地面工作模式时，主要利用SCI接口或485总线接口，通过SCI接口完成工作参数设置，通过485总线接口读取仪器FLASH存储器中的数据。仪器处于延时工作模式时，禁止DSP内部的各种外围设备，并保存关键参数后进入休眠模式；到达设定的时间间隔后，由微控制器通过产生外部中断的方式唤醒DSP，使其进入随钻测井工作模式。随钻测井工作模式下，仪器按照预设的方式工作，启动定时器，采集8扇区或16扇区伽马测量值、动态工具面等参数，并写入缓存中形成一帧数据序列。在微控制器的控制下，按照预设命令同步启动随钻测井工作循环，由微控制器实现大量的控制工作，产生仪器各部分的工作时序，使DSP实时完成每一帧伽马测量值的处理与拟合，并进行数据传输。

图 4 仪器下短节DSP控制及数据处理流程

Figure 4. DSP control and data processing workflow of instrument below the sub

下载: 全尺寸图片幻灯片

井下实时多扇区伽马测量值的拟合采用快速正弦曲线拟合方法，其流程如图5所示。首先，从缓存中读取一帧数据序列中的原始伽马测量值数据，并进行数据预处理操作，剔除奇异值，形成有效伽马测量值数据序列；然后，对伽马测量值数据序列进行频率估计，获取正弦曲线频率的最小二乘估计值；最后，在频率已知的条件下，进行三参数正弦曲线拟合，获得正弦曲线的幅度、频率、相位和直流分量等4个参数。由于钻井液脉冲传输数据的速率较低，大量的原始伽马测量数据和实时计算结果需要保存到井下仪器中的FLASH存储器中，只是将井下实时拟合获得的正弦曲线4个特征参数传输到地面测控系统，再由地面处理程序进行人机交互拟合反演形成8扇区或16扇区伽马成像图，就可以获得比较可信的成像效果。

图 5 伽马测量值快速正弦曲线拟合流程

Figure 5. Rapid sinusoidal fitting process for gamma measurement value

下载: 全尺寸图片幻灯片

4. 近钻头伽马成像测井仪器成像试验

用标准岩样构造分层倾斜地层模拟井眼，进行近钻头伽马成像测井仪器成像试验，以验证快速正弦曲线拟合方法的可行性和正确性。

建立的井眼模型如图6所示。模型外径为800.0 mm，长度为2 200.0 mm；井眼直径为215.9 mm，井眼内充满空气；井眼周围地层由2种不同放射性强度的分层倾斜地层组合而成，高放射性地层A为天然放射性标准花岗岩样混合介质，低放射性地层B为模拟沉积岩混凝土介质；倾斜地层的边界与井眼轴线的夹角为32°。该井眼模型既可以竖直放置，模拟直井分层倾斜地层；也可以水平放置，模拟水平井分层倾斜地层。

图 6 用标准岩样构成的分层倾斜地层井眼模型

Figure 6. Layered dipping strata wellbore model fabricated with standard rock samples

下载: 全尺寸图片幻灯片

井眼模型水平放置时，近钻头伽马成像测井仪器下短节旋转，从井眼的右端移动到井眼的左端进行8扇区伽马成像数据采集，地层放射性强度从250 API变化到50 API，转速30 r/min，水平移动速度0.25 m/min。仪器下短节在井眼内从右端向左端移动时，首先探测到高放射性花岗岩层；然后逐渐向左端移动，从底部靠近倾斜地层边界时，在方位角0°处首先探测到低放射性沉积岩混凝土层；仪器下短节继续向左端移动，其他方位开始探测到低放射性沉积岩混凝土层；最后，仪器下短节完全通过地层边界，并进入到低放射性沉积岩混凝土层中，得到8条方位伽马曲线和8扇区伽马成像图（见图7）。

图 7 8条方位伽马曲线和8扇区伽马成像图

Figure 7. 8 azimuth gamma curves and 8-sector gamma image

下载: 全尺寸图片幻灯片

从图7可以看出，当仪器下短节从高放射性地层穿过倾斜地层边界进入低放射性地层时，伽马测量值从大变小，并且8条方位伽马曲线在地层边界处有明显差异，地层边界在扇区伽马成像图中呈现正弦曲线特征。

采用井下固件程序中的快速正弦曲线拟合算法自动提取正弦曲线4个特征参数，并拟合得到每个深度位置的扇区伽马成像数据曲线（见图8，图中离散点表示方位伽马测量值）。

图 8 8扇区伽马成像数据拟合曲线

Figure 8. Data fitting curve of 8-sector gamma image

下载: 全尺寸图片幻灯片

从图8可以看出，当仪器下短节完全处于高放射性或低放射性地层时，在每个深度位置处测量的扇区伽马成像数据拟合呈直线（见图8中红线和蓝线），各个方位伽马成像测量值基本一致，存在一定的统计涨落误差；当仪器下短节穿过倾斜地层边界时，在每个深度位置处测量的扇区伽马成像数据呈正弦曲线特征，拟合确定系数大于0.75，拟合效果较好（见图8中绿线）。

基于快速正弦曲线拟合算法，利用正弦曲线4个特征参数在地面上反演，得到的8扇区伽马成像和地层界面预测结果如图9所示（图中离散黑点表示伽马成像测量值，紫色曲线表示倾斜地层边界）。

图 9 拟合反演得到的8扇区伽马成像和地层界面预测结果

Figure 9. 8-sector gamma image obtained by fitting inversion and the prediction results of strata interface

下载: 全尺寸图片幻灯片

从图9可以看出，拟合反演得到的8扇区伽马成像图与原始的8扇区伽马成像图基本一致，扇区伽马成像图的正弦曲线特征比较明显，能够准确反映倾斜地面界面信息。依据拟合反演的8扇区伽马成像图，按照式（1）计算出地层倾斜界面与井眼相交的相对倾角为55.6°，即倾斜地层边界与井眼轴线夹角为34.4°，与模拟井眼实际倾斜角的相对误差为7.6%。由此可知，利用快速正弦曲线拟合方法能够提取正弦曲线4个特征参数，进而可以拟合反演形成8扇区伽马成像图，验证了利用快速正弦曲线拟合方法能够准确反演形成扇区伽马成像图。根据倾斜地层界面相对倾角、不同方位伽马曲线中值和深度位置，确定了伽马成像测量值和倾斜地层边界。倾斜地层边界在扇区伽马成像中呈正弦曲线特征，与上文数值模拟获取的8扇区伽马成像图一致，这也验证了利用快速正弦曲线拟合方法反演形成扇区伽马成像图具有可行性。

5. 结论与建议

1）近钻头伽马成像测井仪器穿过倾斜界面进入高放射性泥岩层测量得到的扇区伽马成像图呈现正弦曲线特征，利用扇区伽马成像图可以准确反演井眼相对倾角和围岩厚度，从而确定地层边界，精确控制井眼轨迹。

2）设计了近钻头伽马成像测井中的快速正弦曲线拟合固件算法，该算法将最小二乘估计算法和三参数正弦拟合相结合来获取正弦曲线的幅度、频率、相位和直流分量等4个参数，为拟合反演得到8扇区或16扇区伽马成像图提供数据支持。

3）近钻头伽马成像测井仪器成像试验结果表明，快速正弦曲线拟合方法可以实现扇区伽马成像测量值的拟合，且误差较小。拟合反演得到的8扇区伽马成像图的正弦曲线特征比较明显，能够准确反映倾斜地面界面信息，验证了快速正弦曲线拟合方法的正确性和可行性。

4）在现有随钻测井数据遥测系统传输能力有限的条件下，建议加强快速正弦曲线拟合方法在井下数据压缩技术方面的应用，以提高随钻成像测井技术的实时性。

图 1 泡排剂加入不同助剂后的发泡及稳泡性能

Figure 1. Foaming properties and foaming stability of foamingagent after adding different auxiliary agents

下载: 全尺寸图片幻灯片

图 2 不同温度下的初始泡沫体积

Figure 2. Initial foam volume at different temperatures

下载: 全尺寸图片幻灯片

图 3 不同温度下稳定5 min后的泡沫体积

Figure 3. Foam volume after stabilization for 5 min at differenttemperatures

下载: 全尺寸图片幻灯片

图 4 不同矿化度下稳定5 min后的泡沫体积

Figure 4. Foam volume after stabilization for 5 min at differentsalinity

下载: 全尺寸图片幻灯片

图 5 JY03井生产曲线

Figure 5. Production curves of Well JY03

下载: 全尺寸图片幻灯片

表 1 泡排剂主剂忧选试验结果

Table 1 Experimental results of main agent selection for foaming agent

表面活性剂	V₀/mL	V₅/mL	ω，%
BS-12	450	290	35.56
CAB	480	350	27.08
EHSB	528	305	42.23
OB	425	345	18.82
PKO	460	315	31.52

下载: 导出CSV

表 2 因素水平表

Table 2 Factor level

水平	A	B	C
水平	CAB	OB	Guar
1	1	1	0.05
2	2	2	0.10
3	3	3	0.15
4	4	4	0.20

下载: 导出CSV

表 3 正交试验结果

Table 3 Orthogonal test results

序号		CAB	OB	Guar	V₀/mL	V₅/mL
1		1	1	0.05	450	430
2		1	2	0.10	440	425
3		1	3	0.15	430	420
4		1	4	0.20	425	410
5		2	1	0.10	485	460
6		2	2	0.05	475	455
7		2	3	0.20	470	450
8		2	4	0.15	470	445
9		3	1	0.15	495	465
10		3	2	0.20	480	465
11		3	3	0.05	475	460
12		3	4	0.10	470	455
13		4	1	0.20	490	460
14		4	2	0.15	480	455
15		4	3	0.10	475	450
16		4	4	0.05	475	450
V₀	K₁	634.25	480.00	468.75
	K₂	475.00	468.00	467.50
	K₃	480.00	462.00	468.75
	K₄	480.00	460.00	466.25
	R	43.75	20.00	2.50
V₅	K₁	421.25	453.75	448.75
	K₂	452.50	450.00	447.50
	K₃	461.25	445.00	446.25
	K₄	453.75	440.00	446.25
	R	40.00	13.75	2.50

下载: 导出CSV

表 4 不同COG加量下的发泡性能

Table 4 Foaming properties under different dosage of COG

COG加量，%	V₀/mL	V₅/mL
0.05	440	290
0.10	570	460
0.15	580	465
0.20	600	475
0.25	590	475
0.30	590	475

下载: 导出CSV

表 5 3种泡排剂在页岩气井产出水中的发泡稳泡性能

Table 5 Foaming stability performance of three forming agents in produced-water of shale gas wells

表面活性剂	测试时间	V₀/mL	V₅/mL	ω，%
COG	老化前	575	404	29.74
KD01	老化前	470	274	41.70
XH02	老化前	475	276	41.89
COG	老化后	510	360	29.41
KD01	老化后	420	236	43.81
XH02	老化后	430	245	43.02

下载: 导出CSV

参考文献(21)

[1]	路保平. 中国石化石油工程技术新进展与发展建议[J]. 石油钻探技术,2021,49(1):1–10. doi: 10.11911/syztjs.2021001 LU Baoping. New progress and development proposals of Sinopec’s petroleum engineering technologies[J]. Petroleum Drilling Techni-ques, 2021, 49(1): 1–10. doi: 10.11911/syztjs.2021001
[2]	王友启. 高钙镁油藏烷醇酰胺类表面活性剂及其复配体系驱油性能研究[J]. 石油钻探技术,2018,46(3):98–102. WANG Youqi. Oil displacement performance of alkanolamide surfactants and its compound system for high calcium and magnesium reservoirs[J]. Petroleum Drilling Techniques, 2018, 46(3): 98–102.
[3]	夏海帮,袁航,岑涛. 彭水区块页岩气生产井排采方式研究与应用[J]. 石油钻探技术,2014,42(4):21–26. XIA Haibang, YUAN Hang, CEN Tao. Study and application of drainage methods for shale gas wells in Pengshui Block[J]. Petroleum Drilling Techniques, 2014, 42(4): 21–26.
[4]	王志彬,白慧芳,孙天礼,等. 泡沫排水采气井井下节流压降规律实验及模型修正[J]. 石油钻采工艺,2021,43(3):341–347. WANG Zhibin, BAI Huifang, SUN Tianli, et al. Experimental study and model modification of downhole throttling pressure drop laws in the gas wells with foam drainage gas recovery process[J]. Oil Drilling & Production Technology, 2021, 43(3): 341–347.
[5]	吕伟,刘笑春,白海龙,等. CO₂响应性增强泡沫体系室内试验研究[J]. 石油钻探技术,2021,49(5):88–93. doi: 10.11911/syztjs.2021119 LYU Wei, LIU Xiaochun, BAI Hailong, et al. Laboratory test study of CO₂ responsive enhanced foam system[J]. Petroleum Drilling Techniques, 2021, 49(5): 88–93. doi: 10.11911/syztjs.2021119
[6]	武俊文,熊春明,雷群,等. 阳离子型双子表面活性剂在制备耐高温、高矿化度泡排剂中的应用[J]. 石油钻采工艺,2016,38(2):256–259. WU Junwen, XIONG Chunming, LEI Qun, et al. Application of gmini surfactant in preparing foaming drainage agent with resistance to high temperature and high salinity[J]. Oil Drilling & Production Technology, 2016, 38(2): 256–259.
[7]	赵梓平. 驱油用两性离子型双子表面活性剂的合成及应用[J]. 断块油气田,2019,26(1):119–122. ZHAO Ziping. Synthesis and application of zwitterionic gemini surfactant flooding agent[J]. Fault-Block Oil & Gas Field, 2019, 26(1): 119–122.
[8]	周舰. 低压低产气井井下智能机器人排水采气技术[J]. 石油钻探技术,2020,48(3):85–89. doi: 10.11911/syztjs.2020059 ZHOU Jian. Drainage gas recovery technique with downhole intelligent robots in low pressure and low production gas wells[J]. Petroleum Drilling Techniques, 2020, 48(3): 85–89. doi: 10.11911/syztjs.2020059
[9]	孙玉鹏,吴向阳,张颖,等. 抗凝析油泡排剂SH-1的泡沫性能研究与应用[J]. 应用化工,2021,50(7):1867–1869. doi: 10.3969/j.issn.1671-3206.2021.07.028 SUN Yupeng, WU Xiangyang, ZHANG Ying, et al. Study and application of foam properties of anti-condensate defoaming agent SH-1[J]. Applied Chemical Industry, 2021, 50(7): 1867–1869. doi: 10.3969/j.issn.1671-3206.2021.07.028
[10]	田雨露,王纪伟,李加玉. 气田泡沫排水采气起泡剂研究进展[J]. 应用化工,2021,50(1):183–188. doi: 10.3969/j.issn.1671-3206.2021.01.042 TIAN Yulu, WANG Jiwei, LI Jiayu. Research progress in foaming agent for drainage gas recovery of gas field[J]. Applied Chemical Industry, 2021, 50(1): 183–188. doi: 10.3969/j.issn.1671-3206.2021.01.042
[11]	王伟,崔丹丹,严曦,等. 聚合物表面活性剂与油藏匹配性及液流转向能力研究[J]. 特种油气藏,2021,28(1):111–117. doi: 10.3969/j.issn.1006-6535.2021.01.016 WANG Wei, CUI Dandan, YAN Xi, et al. Study on the compatibility of polymer surfactant with reservoirs and the capability of liquid flow diversion[J]. Special Oil & Gas Reservoirs, 2021, 28(1): 111–117. doi: 10.3969/j.issn.1006-6535.2021.01.016
[12]	李佳欣,张宁波,周成香. 页岩气井泡沫排水采气技术应用研究:以平桥南区为例[J]. 油气藏评价与开发,2020,10(5):91–97. LI Jiaxin, ZHANG Ningbo, ZHOU Chengxiang. Application of foam drainage-gas recovery technology in shale gas wells: A case study of southern Pingqiao[J]. Reservoir Evaluation and Development, 2020, 10(5): 91–97.
[13]	邓强,严娇,白云,等. 表面活性剂起泡动力学方程及温度对其影响研究[J]. 中外能源,2017,22(7):46–51. DENG Qiang, YAN Jiao, BAI Yun, et al. Study on foaming dynamics equation of surfactant and the influence of temperature on it[J]. Sino-Global Energy, 2017, 22(7): 46–51.
[14]	李吉,王江,吴文祥,等. 新型表面活性聚合物驱油剂的研制及应用[J]. 断块油气田,2020,27(6):803–807. LI Ji, WANG Jiang, WU Wenxiang, et al. Development and application of novel surface-active polymer flooding agent[J]. Fault-Block Oil & Gas Field, 2020, 27(6): 803–807.
[15]	郭东红,杨晓鹏,孙建峰,等. 产水气井高效廉价排水采气技术的研究与应用[J]. 现代化工,2018,38(10):137–139. GUO Donghong, YANG Xiaopeng, SUN Jianfeng, et al. Research and application of efficient and economic draining and gas-recovery technique for water-producing natural gas well[J]. Modern Chemical Industry, 2018, 38(10): 137–139.
[16]	秦正山,周建良,谢晶. 致密砂岩油藏注表面活性剂驱油体系提高采收率实验[J]. 断块油气田,2020,27(5):628–632. QIN Zhengshan, ZHOU Jianliang, XIE Jing. Experiment on enhanced oil recovery by surfactant-injection oil displacement system in tight sandstone reservoirs[J]. Fault-Block Oil & Gas Field, 2020, 27(5): 628–632.
[17]	梅绪东,金吉中,王朝强,等. 涪陵页岩气田绿色开发的实践与探索[J]. 西南石油大学学报(社会科学版),2017,19(6):9–14. MEI Xudong, JIN Jizhong, WANG Chaoqiang, et al. Practice and research of green development in Fuling shale gas field[J]. Journal of Southwest Petroleum University(Social Sciences Edition), 2017, 19(6): 9–14.
[18]	陈晓宇. 泡沫排水采气技术在涪陵页岩气田的应用[J]. 天然气技术与经济,2019,13(5):49–53. CHEN Xiaoyu. Technologies of foam drainage gas recovery and their application to Fuling Shale-Gas Field[J]. Natural Gas Technology and Economy, 2019, 13(5): 49–53.
[19]	ADEBAYO A R. Sequential storage and in-situ tracking of gas in geological formations by a systematic and cyclic foam injection: a useful application for mitigating leakage risk during gas injection[J]. Journal of Natural Gas Science and Engineering, 2019, 62: 1–12. doi: 10.1016/j.jngse.2018.11.024
[20]	GHORBANI S, SHARIFI S, DE BRITO J, et al. Using statistical analysis and laboratory testing to evaluate the effect of magnetized water on the stability of foaming agents and foam concrete[J]. Construction and Building Materials, 2019, 207: 28–40. doi: 10.1016/j.conbuildmat.2019.02.098
[21]	WANG Hetang, LI Jia, WANG Zhan, et al. Experimental investigation of the mechanism of foaming agent concentration affecting foam stability[J]. Journal of Surfactants and Detergents, 2017, 20(6): 1443–1451. doi: 10.1007/s11743-017-2004-2

施引文献(11)

期刊类型引用(7)

1.	王忠良，徐文远，王文泽，梁梅，贺杰，钱晨，邹华宝. 结合伽马成像技术的三维地质建模在页岩油地质导向中的应用. 西安石油大学学报(自然科学版). 2024(02): 112-119 . 百度学术
2.	曲博文，谭宝海，张凯，陈雪莲. 自适应声波测井换能器激励电路设计. 石油钻探技术. 2024(06): 141-147 . 本站查看
3.	孙轶伦，潘磊，李海滨，王来智，周建新. 旋转导向近钻头随钻伽马数据成像处理技术. 石油地质与工程. 2022(01): 99-103+108 . 百度学术
4.	刘卫彬，徐兴友，张君峰，陈珊，白静，刘畅，李耀华. 陆相页岩地层地质-工程一体化水平井精确钻探技术——以松辽盆地吉页油1HF井为例. 中国地质. 2022(06): 1808-1822 . 百度学术
5.	马荣华，赵嵩. 基于果蝇算法的空间机械臂最优抑振轨迹规划方法. 机械设计与研究. 2022(06): 49-52 . 百度学术
6.	李继博，钱德儒，郑奕挺，张卫，吴金平. 近钻头伽马高精度实时成像技术研究与应用. 石油钻探技术. 2021(03): 135-141 . 本站查看
7.	李洪强，王瑞和. 近钻头伽马成像仪测量结果环境校正方法研究. 石油钻探技术. 2021(03): 142-150 . 本站查看

其他类型引用(4)

资源附件(0)

图(5) / 表(5)

计量

文章访问数: 273
HTML全文浏览量: 181
PDF下载量: 35
被引次数: 11

1. 倾斜界面近钻头伽马成像特征
2. 快速正弦曲线拟合方法
3. 快速正弦曲线拟合固件算法的实现
4. 近钻头伽马成像测井仪器成像试验
5. 结论与建议

二元复合型泡排剂COG的研制及现场试验

作者简介: 王维恒（1989—），男，江苏南京人，2011年毕业于中国地质大学（北京）材料化学专业，2014年获中国地质大学（北京）应用化学专业工学硕士学位，助理研究员，从事油田化学助剂开发工作。E-mail：beiton@163.com

计量

出版历程