生成对抗网络算法在电成像测井裂缝空白条带填充中的应用

康正明, 武辰升, 杨国栋, 武迪生, 王瑞飞, 杨湘雨, 甘伟

康正明,武辰升,杨国栋,等. 生成对抗网络算法在电成像测井裂缝空白条带填充中的应用[J]. 煤田地质与勘探,2025,53(3):197−207. DOI: 10.12363/issn.1001-1986.24.07.0472
引用本文: 康正明,武辰升,杨国栋,等. 生成对抗网络算法在电成像测井裂缝空白条带填充中的应用[J]. 煤田地质与勘探,2025,53(3):197−207. DOI: 10.12363/issn.1001-1986.24.07.0472
KANG Zhengming,WU Chensheng,YANG Guodong,et al. Application of a generative adversarial network algorithm to filling blank strips of fractures in formation microresistivity imaging images[J]. Coal Geology & Exploration,2025,53(3):197−207. DOI: 10.12363/issn.1001-1986.24.07.0472
Citation: KANG Zhengming,WU Chensheng,YANG Guodong,et al. Application of a generative adversarial network algorithm to filling blank strips of fractures in formation microresistivity imaging images[J]. Coal Geology & Exploration,2025,53(3):197−207. DOI: 10.12363/issn.1001-1986.24.07.0472

 

生成对抗网络算法在电成像测井裂缝空白条带填充中的应用

基金项目: 国家自然科学基金项目(42304138);陕西省教育厅重点科学研究计划项目(22JY053);西安石油大学研究生创新与实践能力培养 计划(YCS23215348)
详细信息
    作者简介:

    康正明,1989年生,男,陕西靖边人,博士,副教授。E-mail:190720@xsyu.edu.cn

  • 中图分类号: P631

Application of a generative adversarial network algorithm to filling blank strips of fractures in formation microresistivity imaging images

  • 摘要:
    目的 

    电成像测井仪器极板之间存在间隔,测得的井壁电阻率图像存在空白条带,对井壁附近裂缝参数评价影响较大。针对裂缝空白条带填充没有全井眼图像,填充质量评价难的问题,提出一种模拟数据与真实数据相结合构造数据集,基于生成对抗网络的电成像测井裂缝空白条带填充方法。

    方法 

    首先,基于三维有限元法对裂缝地层的电阻率测井响应进行模拟。采用过渡边界条件代替实体裂缝的方法提升了多尺度模型的计算效率,满足了深度学习对样本的需求。其次,将模拟得到的裂缝与实际测井图像相结合构建了全井眼图像,用电成像测井图像制作的掩码作为空白条带。通过图像填充指标评价裂缝空白条带填充能力的方式优选超参数,获得最佳空白条带填充效果。最后,采用不同占比的掩码进行裂缝空白条带填充,并评价裂缝填充质量。

    结果和结论 

    结果表明:模型能够在不同占比的空白条带情况下进行填充,对粗糙裂缝有较好的填充效果且裂缝形态平滑,能更准确地恢复边缘轮廓与细节。在真实测井数据上的填充进一步验证了其适用性,填充后的图像自然且能有效恢复裂缝特征,有助于裂缝提取及定量计算,为裂缝储层的精准评价提供基础,支撑油气产能的准确预测。

    Abstract:
    Objective 

    Gaps between the electrodes of formation microresistivity imaging (FMI) imagers lead to blank strips in the resistivity images of borehole walls, significantly influencing the parameter assessment for fractures near borehole walls. The absence of full-borehole images renders it challenging to assess the blank strip filling quality for fractures in FMI images. Using a dataset constructed utilizing both simulated and actual data, this study proposed a generative adversarial network (GAN)-based method for filling the blank strips of fractures in FMI images.

    Methods 

    First, the resistivity logging responses of a fractured formation were simulated using the 3D finite element method. Actual fractures were simplified using the transition boundary condition, improving the computational efficiency by avoiding multi-scale models and meeting the requirements of deep learning for substantial samples. Second, full-borehole images were attained by integrating the simulated fractures with actual FMI images, and masks generated from FMI images were used as blank strips. To achieve the optimal filling effects, hyperparameters were optimized by assessing the capacity of blank strip filling for fractures using image filling indicators. Finally, the blank strips of fractures were filled using masks with varying proportions, followed by a filling quality assessment.

    Results and Conclusions 

    The results indicate that the proposed method allows for the filling of blank strips with varying proportions. Notably, it yields encouraging filling effects and smooth morphology for coarse fractures and can restore their margins and details more accurately. Its applicability was further verified using actual log data, revealing that it produced natural-looking filled images that can effectively restore fracture features. The proposed method facilitates the extraction and quantitative analysis of fractures, laying a foundation for the precise assessment of fractured reservoirs and supporting the accurate prediction of hydrocarbon productivity.

  • 井筒气水两相管流节流前后动态特性分析与控制,对保障致密气等含水气井高效稳定安全开采和解决因水合物淤积造成管阻问题具有重要意义[1-3]。由于气井井下节流喷嘴的降压作用,井口压力得到有效调节,井口设备与管道承压下降[4-6];此外,节流后的含水气流一般处于临界流动状态,井口设备的压力波动不会对井底压力产生影响,减少了压力波动和井底出砂。

    对于气液两相流态节流过程,目前主要通过研究节流中流体质量、动量、能量与相态的变化,提出了经验和机理相结合的算法。Shao Huaishuang等[7]在Perkins和Al-safran模型基础上,提出一种节流阀多相流态和质量流量预测方法,该方法考虑了流体经喷嘴节流时的交替气液绝热条件和气液两相间的滑移现象,但模型中两相流为均相流,且仅适用于临界流与亚临界流。C. M. Carstensen等[8]在Gilbert和Ros模型基础上,修正非平衡系数并解释节流阀快速降压的中间状态,提出运用非平衡流动的概念预测冻结流动和绝热膨胀的节流过程,但却忽略了组分、温压、体积与流体温度等因素的影响。王振嘉等[9]综合产量对比分析法和油套环空压差预测法,并结合井下节流喷嘴尺寸、下井深度等参数,提出一种适用于“低渗透率、低孔隙率、低丰度、强非均质性”特征储层井下节流气井的积液量预测方法,其适用范围受限。安永生等[10]建立了井下节流气井渗流、喷嘴气流数学模型,并依此推导数值求解源汇项方程,完成气井渗流和喷嘴气流的耦合处理,仅适用于遵守达西定律的单相气体流动,且假设节流前后为等温过程。曾焱等[11]建立了节流气井储层和井筒耦合生产动态模型,拟合试井解释或生产动态数据,并依此预测节流气井的流量和压力变化情况,但仅适用于单相气体流动工况。

    S. Naseri等[12]在GA-RBF模型基础上,提出一种亚音速气流经喷嘴和孔板式节流阀时人工智能预测流量系数方法,该方法采用径向基函数神经网络结合遗传算法,并考虑了雷诺数以及节流喷嘴直径与管径比值等因素。王荧光等[13]模拟井下节流过程压力和温度变化情况,依据NETOPT网络模型优化节流喷嘴内径等结构参数。这些研究可归纳为网络优化模型,但由于算法的实际应用效果欠佳,目前工程上还没有得到普遍应用。

    目前,大部分研究成果不适用于高压含水气井或存在工程上应用受限等问题,笔者综合等熵绝热膨胀、焓与动能非平衡转化等过程,提出一种适用于致密气、湿气井等高压含水气井开采工况下的井筒气水两相管流节流前后动态特性分析与控制方法,并通过数值模拟和在鄂尔多斯盆地东缘大宁−吉县区块井场验证,分析高压含水气井节流过程随喷嘴尺寸、含水率以及节流前压力和温度动态变化状况等,以期达到保证高压含水气井高效稳定开采和防治水合物的目的。

    气水两相流经节流喷嘴时,井下节流过程的能量变化基本方程[14-15]为:

    $$ Q = ({H_{{\text{out}}}} - {H_{{\text{in}}}}) + 0.051J(v_{{\text{out}}}^{\text{2}} - v_{{\text{in}}}^2) + JW $$ (1)

    气水两相流体流经节流喷嘴的时间较短,节流前后气水两相流体几乎不与井筒及其近井围岩发生热功和机械功的交换,则节流过程气水两相流体的能量组成为:

    $$ {H_{{\text{out}}}} - {H_{{\text{in}}}} = 0.051J(v_{{\text{in}}}^2 - v_{{\text{out}}}^{\text{2}}) $$ (2)

    由此推导出喷嘴出口处两相流的速度为:

    $$ {v_{{\text{out}}}} = \sqrt {\frac{{19.62}}{J}\left( {{H_{{\text{in}}}} - {H_{{\text{out}}}}} \right) + v_{{\text{in}}}^2} $$ (3)

    井下节流前后,流经喷嘴的气水两相流体的焓转化为动能,节流后两相流速得以提升,等熵绝热过程中,节流后的两相流体压力下降,而比热容增大;与此同时,节流后两相流体的熵不变,但两相流体的温度下降。

    不考虑井下节流过程的位能差、系统与外界的热交换,且节流喷嘴内部流体处于等熵绝热状态,依据节流过程能量守恒方程[16-18],气水两相管流井下节流过程的能量变化(图1)为:

    图  1  高压含水气井筒流体能量随井深变化情况
    Figure  1.  Depth-varying fluid energy in wellbores of high-pressure gas wells producing water
    $$ {p}_{\text{in}}{V}_{\text{in}}+\frac{{v}_{\text{in}}^{2}}{2}+{C}_{\mathrm{V}}({T}_{\text{in}}-{T}_{\text{out}})={p}_{\text{out}}{V}_{\text{out}}+\frac{{v}_{\text{out}}^{2}}{2} $$ (4)

    单位质量气水两相流体的能量组成为:

    $$ f_{\text{λ}} ({p_{{\text{in}}}}{V_{{\text{in}}}} - {p_{{\text{out}}}}{V_{{\text{out}}}}) + \frac{{{f_{{{\mathrm{l}}}}}}}{{{\rho _{{{\mathrm{l}}}}}}}({p_{{\text{in}}}} - {p_{{\text{out}}}}) + \frac{{v_{{\text{in}}}^{\text{2}} - v_{{\text{out}}}^{\text{2}}}}{2} = 0 $$ (5)
    $$ f_{\text{λ}} ={f}_{\mathrm{g}}+\frac{({f}_{\mathrm{g}}{C}_{\mathrm{V}\mathrm{g}}+{f}_{{{\mathrm{l}}}}{C}_{\mathrm{Vl}{}}){M}_{\mathrm{g}}}{ZR} $$ (6)

    节流过程中系统热量Q保持不变,则有:

    $$ {\text{d}}Q = {f_{\text{g}}}({C_{{\mathrm{V}}{{{\mathrm{g}}}}}}{\text{d}}T + p{\text{d}}V) + {f_{{{\mathrm{l}}}}}{C_{{\mathrm{V}}{{{\mathrm{l}}}}}}{\text{d}}T = 0 $$ (7)

    基于熵的流体热力学微分方程为:

    $$ p{\text{d}}{V_{\text{M}}} + {V_{\text{M}}}{\text{d}}p = \frac{{ZR}}{{{M_{\text{g}}}}}{\text{d}}T $$ (8)

    由此推导出节流过程系统热量微分方程为:

    $$\begin{aligned} &\left[{f}_{{\mathrm{g}}}\left({C}_{\mathrm{V}\mathrm{g}}{\mathrm{d}}T+\frac{ZR}{{M}_{\mathrm{g}}}\right)+{f}_{{\mathrm{l}}}{C}_{{\mathrm{Vl}}}\right]p{\mathrm{d}}V+\\ &\qquad ({f}_{\mathrm{g}}{C}_{\mathrm{V}\mathrm{g}}+{f}_{\mathrm{l}}{C}_{\mathrm{V}\mathrm{l}}){V}_{\mathrm{M}}{\mathrm{d}}p=0\end{aligned}$$ (9)

    井下节流过程中,气相等压热容(Cp)[19-20]与气相等容热容(CVg)之间的关系式为:

    $$ {C_{\text{p}}} = {C_{{\text{Vg}}}} + \frac{{ZR}}{{{M_{\text{g}}}}} $$ (10)

    考虑等压热容、等容热容等因素后的节流系统热量为:

    $$ {p_{{\text{in}}}}V_{{\text{in}}}^n = {p_{{\text{out}}}}V_{{\text{out}}}^n $$ (11)

    其中,

    $$ n = \frac{{{f_{\text{g}}}{C_{\text{p}}}C_{_{{\text{Vg}}}}^{ - 1} + {f_{\text{l}}}{C_{{\text{Vl}}}}}}{{{f_{\text{g}}}{C_{{\text{Vg}}}} + {f_{\text{l}}}{C_{{\text{Vl}}}}}} $$ (12)

    引入节流后出口与节流前入口的压力比β,减小压力比β,节流过程气水两相流体质量流量有增大的趋势,且节流后出口压力降至一定程度后,流体速度存在音速界限并保持不变,即两相流体达到音速的临界流动状态[21-24],临界流动状态单位质量两相流体能量为:

    $$ f_{\text{λ}} {p_{{\text{in}}}}{V_{{\text{in}}}}\left[1 - {\beta ^{\tfrac{{n - 1}}{n}}}\right] + \frac{{{f_{\text{l}}}}}{{{\rho _{\text{l}}}}}{p_{{\text{in}}}}(1 - \beta ) + \frac{{v_{{\text{out}}}^{\text{2}}}}{2}\left[ {{{\left( {\frac{1}{\beta }} \right)}^2} - 1} \right] = 0 $$ (13)
    $$ \beta = \frac{1}{V}{\left( {\frac{{\partial V}}{{\partial T}}} \right)_{\text{p}}} $$ (14)

    井下节流过程中物质的平衡方程[25-28]为:

    $$ \frac{{{A_{{\text{in}}}}{v_{{\text{in}}}}}}{{{f_{\text{g}}}{V_{{\text{in}}}} + {f_{\text{l}}}/{\rho _{\text{l}}}}} = \frac{{{A_{{\text{out}}}}{v_{{\text{out}}}}}}{{{f_{\text{g}}}{V_{{\text{out}}}} + {f_{\text{l}}}/{\rho _{\text{l}}}}} $$ (15)

    井下节流过程流体的质量流量为:

    $$ w = \frac{{Av}}{{{f_{\text{g}}}{V_{{\text{out}}}} + {f_{\text{l}}}/{\rho _{\text{l}}}}} $$ (16)

    井下节流入口压力与出口压力之间的压力比低于临界压力比[29-31]时,气水两相流体处于临界流动状态,质量流量达到最大并保持不变,即:

    $$ \frac{{\text{d}}}{{{\text{d}}\beta }}{\left( {\frac{w}{{{A_{{\text{out}}}}\sqrt {{p_{{\text{in}}}}/{V_{{\text{in}}}}} }}} \right)^2} = 0 $$ (17)

    选取鄂尔多斯盆地东缘大宁−吉县区块致密气高压含水气井开展井筒气水两相流态井下节流前后动态特性分析。依据井下技术作业公司的检测报告,大宁−吉县区块井场测试的井筒各部分物性参数与生产参数为:节流前压力14、15、16、17、18 MPa,节流前温度313、318、323、328、333 K,气体密度0.68 kg/m3,含水率0、0.1、0.2、0.3、0.4 mol/mol,喷嘴内径3.0、3.5、4.0、4.5、5.0 mm,油管内径63 mm,油套内径139.7 mm。

    高压含水气井节流前后两相流的动态特性会受到喷嘴尺寸等结构参数的影响。节流前压力18 MPa、温度333 K和含水率0.1 mol/mol时,气水两相流体节流过程质量流量随喷嘴内径和压力比的动态变化关系如图2所示。减小压力比β,节流过程两相流体质量流量先是呈指数级增大而后基本保持恒定,压力比的阈值为0.55左右。

    图  2  高压含水气井节流过程随喷嘴尺寸动态变化状况
    Figure  2.  Dynamic variations of the throttling process in high-pressure gas wells producing water under different inner diameters and setting depths of chokeling nozzles

    喷嘴内径等结构参数对高压含水气井两相流节流过程质量流量的影响相对较大,且给定压力比β时,节流过程质量流量随喷嘴内径的不断增大而显著提升,这与安永生等[10]所得出的“井下节流喷嘴直径越大,气井产量越高”的结论一致。图2a中,喷嘴内径由3.0 mm扩至5.0 mm时,气水两相流体的临界质量流量由1.45×104 kg/d增至4.05×104 kg/d,临界质量流量提升幅度达到179.3%。图2b给出不同压力比β工况节流后温度变化情况,当β=1时,节流后温度等于节流前温度;而后随压力比不断减小,节流后温度呈现类“抛物线”式下降的趋势。其主要原因在于节流前入口压力一定的情况下,节流后出口压力的下降会造成压力比β的减小,节流过程气水两相流体质量流量有增大的趋势,综合等压热容、等容热容等因素作用,节流之后出现两相流体温度降低的现象。此外,节流前入口压力一定时,喷嘴结构参数设计不合理易使得节流后出口压力发生“低压”工况(图2b中压力比小于0.1),由此会造成节流之后的两相流体温度降至零度以下;而当节流后出口压力一定时,不断增加典型“高压”含水气井的节流前入口压力(14~18 MPa),也易出现节流后两相流体温度降至零度以下的现象。

    图2c中,给定压力比β时,喷嘴下入深度越深,节流后温度越高,这与曾焱等[11]所得出的“井下节流温度随喷嘴深度的不断增大而逐渐提升”的结论一致。此外,图2c中气水两相流体节流后温度随喷嘴下入深度和压力比的不断减小而呈现类“抛物线”式逐渐下降的趋势;且喷嘴下入深度较浅且压力较小时,节流之后的两相流体温度可降至−293 K~−333 K而出现低温工况,其主要原因在于高压含水气井节流前后两相流动态特性分析充分考虑了等熵绝热、等压热容、等容热容、不同井深节流能量等因素的综合作用,而且低温高压工况的出现易导致水相态发生变化并伴随水合物生成的情况,从而不利于水合物生成的防治;图2c中喷嘴下入深度400 m且压力比0.01的工况下,气水两相流体的节流后温度最低可降至−325.7 K。为此,需合理设计高压含水气井喷嘴的下入深度,避免出现井深超1.6 km而节流作业中喷嘴下入深度不足1.0 km的情况;并且需合理制定节流前后的压力等生产制度,防止发生因压力比过小而造成节流后温度过低等问题;同时,充分利用地层温度对两相流体进行升温,以保证气水两相流体节流过程中的温度和节流后的温度高于水合物的临界生成温度。

    高压含水气井节流前后的动态特性也受到压力等生产参数的影响,喷嘴内径5.0 mm、节流前温度333 K和含水率0.1 mol/mol工况下,气水两相流体节流过程质量流量以及节流后温度随节流前压力和压力比的动态变化关系如图3所示。

    图  3  高压含水气井节流过程随压力动态变化状况
    Figure  3.  Dynamic variations of the chokeling process in high-pressure gas wells producing water under different pressures

    高压含水气井气水两相流体节流过程质量流量随节流前压力的增大而缓慢提升,且较喷嘴内径的影响较低,图3a中,节流前压力由14 MPa增至18 MPa时,气水两相流体的临界质量流量由3.17×104 kg/d升至4.05×104 kg/d,临界质量流量提升的幅度仅为27.8%。图3b中,节流过程给定压力比β时,节流前压力越高,节流后温度越低,且其影响程度随压力比下降而显著增大。

    高压含水气井节流前后的动态特性还受到温度等生产参数的影响,喷嘴内径5.0 mm、节流前压力18 MPa和含水率0.1 mol/mol工况下,气水两相流体节流过程质量流量以及节流温差随节流前温度和压力比的动态变化关系如图4所示。

    图  4  高压含水气井节流过程随温度动态变化状况
    Figure  4.  Dynamic variations of the chokeling process in high-pressure gas wells producing water under different temperatures

    节流前温度等生产参数对高压含水气井气水两相流体节流过程质量流量以及节流温差的影响均较小,节流过程质量流量随节流前温度的升高而略有降低,图4a中,节流前温度由313 K升至333 K时,气水两相流体的临界质量流量由4.27×104 kg/d降至4.05×104 kg/d,下降幅度仅为5.15%。图4b中,节流过程给定压力比β时,节流前温度越高,节流温差ΔT=ToutTin越低,且其影响程度随压力比下降而显著增大。

    高压含水气井节流前后的动态特性同时会受到含水率等物性参数的影响,喷嘴内径5.0 mm、节流前压力18 MPa和温度333 K工况下,气水两相流体节流过程质量流量以及节流后温度随含水率和压力比的动态变化关系如图5所示。

    图  5  高压含水气井节流过程随含水率动态变化状况
    Figure  5.  Dynamic variations of the chokeling process in high-pressure gas wells producing water under different water contents

    含水率等物性参数对高压含水气井气水两相流体节流过程质量流量的影响相对较小,与节流前压力的影响接近,图5a中,含水率由0增至0.4 mol/mol时,气水两相流体的临界质量流量由4.04×104 kg/d升至5.01×104 kg/d,增幅为23.76%。图5b中,节流过程给定压力比时,含水率越高,节流后温度也越高。

    气水两相流体井筒运移节流过程临界质量流量受喷嘴内径的影响最大,这与安永生等[10]所提出的“定气嘴尺寸生产”的概念一致,含水率和节流前压力的影响次之,而节流前温度和喷嘴下入深度的影响最小。

    相较于安永生、曾焱等前人研究[9-11],一方面笔者团队综合等熵绝热膨胀、焓与动能非平衡转化等过程,提出了一种适用于致密气等高压含水气井开采工况和气水两相流体节流前后动态特性分析方法,而不是仅适用于单相气体流动工况、特定特征储层节流积液量预测或假设节流前后为等温过程;另一方面,本文不仅量化分析了喷嘴内径等结构参数对节流过程动态特性的影响,而且揭示出含水率和喷嘴深度等物性参数以及压力和温度等生产参数对高压含水气井节流特性的具体影响程度,初步解决了现有算法实际应用效果欠佳和工程上应用受限的问题,并为优化井下节流器及其节流喷嘴的结构参数和保障气井流动安全提供了理论依据。

    a. 综合等熵绝热、等压热容、等容热容、不同井深节流能量等因素,推导出气水两相流体经喷嘴的能量、动能和温度动态变化等热力学微分方程组,提出了一种高压气水两相节流动态特性分析与节流控制方法,为优化井下节流器及其节流喷嘴的结构参数和保障气井流动安全提供了理论依据,并可避免致密气等高压含水气井因水合物淤积而引起的管阻问题。

    b. 增大喷嘴内径及其下入深度和节流前压力同时降低节流前温度,有利于提高两相流体节流过程中临界质量流量,并提升高压含水气井的产气量;而增加含水率虽然会小幅增加临界质量流量,但产气量也随之下降;且节流临界质量流量受喷嘴内径的影响最大,含水率和节流前压力次之,而节流前温度和喷嘴深度的影响最小。

    c. 在高压含水气井的实际生产中,应合理设计喷嘴下入深度和节流前后压力,并控制含水率和节流温度,同时充分利用地层温度对两相流体进行升温,达到保障高压含水气井高效稳定开采和防治水合物的目的。

    d. 致密气、湿气井等高压含水气井开采工况受多因素影响,且高压气水两相流体节流前后流场及运移规律复杂,本次研究成果仅用于高压含水气井节流前后两相流动态特性分析,后续可继续开展高压含水气井不同开采工况喷嘴动态节流场仿真分析和节流喷嘴优化设计工作。

    符号注释:

    A为连接节流器的油管流道横截面积,m2AinAout分别为节流前喷嘴入口处和节流后喷嘴出口处的节流器流道横截面积,m2CV为流体的等容热容,J/(mol∙K);CpCVgCVl分别为气相等压热容、等容热容和水相等容热容,J/(mol∙K);fgfl分别为气相与水相的滞留率;fλ为滞留率;HinHout分别为喷嘴入口和出口处流体焓,J/kg;J为热功当量,J/kg;Mg为气相的相对分子质量,J/(kg∙K);n为节流系统热量系数;p为流体压力,MPa;pinpout分别为节流前后流体压力,MPa;Q为总交换热,J/kg;R为气体通用常数,J/(mol∙K);T为流体温度,K;TinTout分别为节流前后流体温度,K;v为节流过程中流体速度,m/s;vinvout分别为喷嘴入口和出口处流体的速度,m/s;V为流体的摩尔体积,L/mol;VM为井筒壁面流体的摩尔体积,L/mol;VinVout分别为节流前后流体的摩尔体积,L/mol;w为节流过程流体的质量流量,kg/d;W为总交换功,J/kg;Z为压缩因子;ρl为水相密度,kg/m3

  • 图  1   裂缝地层数值计算模型

    Fig.  1   Numerical calculation model of a fractured formation

    图  2   裂缝成像与数据扩充

    Fig.  2   Fracture imaging and data augmentation

    图  3   数据集制备方法

    Fig.  3   Method of the sample data preparation

    图  4   空白带填充模型工作流程

    Fig.  4   Workflow of the blank strip filling model

    图  5   空白带填充模型网络架构

    Fig.  5   Network architecture of the blank strip filling model

    图  6   不同学习率和批尺寸下测试集均方误差曲线对比

    Fig.  6   Comparison of MSE curves under different learning rates and batch sizes

    图  7   网络训练多轮迭代输出图像

    Fig.  7   Images output after multiple rounds of training

    图  8   粗糙裂缝填充结果

    Fig.  8   Filling results of rough fractures

    图  9   测井图像填充结果

    Fig.  9   The filling result of the Logging image

    表  1   不同批尺寸和学习率的均方误差

    Table  1   MSEs under varying batch sizes and learning rates

    $\eta $ 训练集 测试集
    4 8 16 32 4 8 16 32
    0.000 1 0.000 584 0.000 725 0.000 837 0.001 028 0.000 908 0.001 040 0.001 183 0.001 379
    0.000 3 0.000 580 0.000 631 0.000 706 0.001 013 0.000 841 0.000 926 0.001 177 0.001 291
    0.000 5 0.000 594 0.000 587 0.000 624 0.001 011 0.000 865 0.001 196 0.001 063 0.001 277
    0.001 0 0.000 648 0.000 704 0.000 948 0.001 290 0.000 882 0.001 199 0.001 272 0.001 441
    0.003 0 0.000 876 0.001 041 0.001 278 0.001 438 0.001 044 0.001 399 0.001 416 0.001 894
    0.005 0 0.001 004 0.001 038 0.001 380 0.002 179 0.001 200 0.001 677 0.001 472 0.002 636
    下载: 导出CSV

    表  2   典型仪器掩码示例

    Table  2   Examples of masks of typical imagers

    掩码类型 FMI STARⅡ XRMI
    掩码
    下载: 导出CSV

    表  3   不同掩码情况下空白带网络填充性能

    Table  3   Performance of GAN-based blank strip filling under varying mask types

    掩码类型峰值信噪比/dB结构相似性
    FMI36.307 60.984 4
    STARⅡ32.421 30.956 9
    XRMI30.98310.9421
    下载: 导出CSV
  • [1]

    EKSTROM M P,DAHAN C A,CHEN Minyi,et al. Formation imaging with microelectrical scanning arrays[C]//The SPWLA 27th Annual Logging Symposium. Houston,Texas,June 9−13,1986:SPWLA-1986-BB.

    [2] 侯振学,于雪娴,李东旭,等. 电成像测井处理新技术在储层评价方面的应用[J]. 地球物理学进展,2020,35(2):573−578.

    HOU Zhenxue,YU Xuexian,LI Dongxu,et al. Application of new processing technology of electrical imaging logging in reservoir evaluation[J]. Progress in Geophysics,2020,35(2):573−578.

    [3] 邹才能,朱如凯,吴松涛,等. 常规与非常规油气聚集类型、特征、机理及展望:以中国致密油和致密气为例[J]. 石油学报,2012,33(2):173−187. DOI: 10.1038/aps.2011.203

    ZOU Caineng,ZHU Rukai,WU Songtao,et al. Types,characteristics,genesis and prospects of conventional and unconventional hydrocarbon accumulations:Taking tight oil and tight gas in China as an instance[J]. Acta Petrolei Sinica,2012,33(2):173−187. DOI: 10.1038/aps.2011.203

    [4] 黄振华,程礼军,刘俊峰,等. 微电阻率成像测井在识别页岩岩相与裂缝中的应用[J]. 煤田地质与勘探,2015,43(6):121−123. DOI: 10.3969/j.issn.1001-1986.2015.06.025

    HUANG Zhenhua,CHENG Lijun,LIU Junfeng,et al. Application of micro–resistivity image logging in identifying shale facies and fractures[J]. Coal Geology & Exploration,2015,43(6):121−123. DOI: 10.3969/j.issn.1001-1986.2015.06.025

    [5]

    SIBBIT A M,FAIVRE O. The dual laterolog response in fractured rocks[C]//The SPWLA 26th Annual Logging Symposium. Dallas,Texas,June 17−20,1985:SPWLA-1985-T.

    [6]

    RITTER R N,CHEMALI R,LOFTS J,et al. High resolution visualization of near wellbore geology using while–drilling electrical images[C]//The SPWLA 45th Annual Logging Symposium. Noordwijk,Netherlands,June 6−9,2004:SPWLA-2004-PP.

    [7] 蔡亚琳,柯式镇,康正明,等. 随钻电阻率成像测井在裂缝地层中的响应模拟[J]. 石油科学通报,2020,5(3):327−336. DOI: 10.3969/j.issn.2096-1693.2020.03.028

    CAI Yalin,KE Shizhen,KANG Zhengming,et al. Logging response simulation of a LWD resistivity imaging tool in fractured formations[J]. Petroleum Science Bulletin,2020,5(3):327−336. DOI: 10.3969/j.issn.2096-1693.2020.03.028

    [8]

    KANG Zhengming,LI Xin,NI Weining,et al. Using logging while drilling resistivity imaging data to quantitatively evaluate fracture aperture based on numerical simulation[J]. Journal of Geophysics and Engineering,2021,18(3):317−327. DOI: 10.1093/jge/gxab016

    [9]

    ZHANG Tuanfeng,SWITZER P,JOURNEL A. Filter–based classification of training image patterns for spatial simulation[J]. Mathematical Geology,2006,38:63−80. DOI: 10.1007/s11004-005-9004-x

    [10]

    HURLEY N F,ZHANG Tuanfeng. Method to generate full–bore images using borehole images and multipoint statistics[J]. SPE Reservoir Evaluation & Engineering,2011,14(2):204−214.

    [11] 孙建孟,赵建鹏,赖富强,等. 电测井图像空白条带填充方法[J]. 测井技术,2011,35(6):532−537. DOI: 10.3969/j.issn.1004-1338.2011.06.008

    SUN Jianmeng,ZHAO Jianpeng,LAI Fuqiang,et al. Methods to fill in the gaps between pads of electrical logging images[J]. Well Logging Technology,2011,35(6):532−537. DOI: 10.3969/j.issn.1004-1338.2011.06.008

    [12] 李振苓,沈金松,李思,等. 成像测井电导率图像空白带奇异谱插值和缝洞孔隙度分离方法[J]. 测井技术,2017,41(1):33−40.

    LI Zhenling,SHEN Jinsong,LI Si,et al. Singular spectral interpolation of blank strips in formation micro–scanner conductivity image and separation of porosity between fracture and karst cave[J]. Well Logging Technology,2017,41(1):33−40.

    [13] 张翔,张猛,肖小玲,等. 复杂地层情况下全井周电成像图像修复方法[J]. 石油物探,2018,57(1):148−153. DOI: 10.3969/j.issn.1000-1441.2018.01.019

    ZHANG Xiang,ZHANG Meng,XIAO Xiaoling,et al. Image inpainting for fullbore electrical imaging logging in complex formations[J]. Geophysical Prospecting for Petroleum,2018,57(1):148−153. DOI: 10.3969/j.issn.1000-1441.2018.01.019

    [14]

    KRIZHEVSKY A,SUTSKEVER I,HINTON G E. ImageNet classification with deep convolutional neural networks[J]. Advances in Neural Information Processing Systems,2012,25:1−9.

    [15]

    SCHROFF F,KALENICHENKO D,PHILBIN J. FaceNet:A unified embedding for face recognition and clustering[C]//2015 IEEE Conference on Computer Vision and Pattern Recognition (CVPR). Boston,MA,USA,June 07−12,2015. doi:10.1109/CVPR.2015.7298682.

    [16]

    DONG Chao,LOY C C,HE Kaiming,et al. Image super–resolution using deep convolutional networks[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence,2016,38(2):295−307. DOI: 10.1109/TPAMI.2015.2439281

    [17] 王哲峰,高娜,曾蕊,等. 基于深度学习模型的测井电成像空白条带充填方法[J]. 测井技术,2019,43(6):578−582.

    WANG Zhefeng,GAO Na,ZENG Rui,et al. A gaps filling method for electrical logging images based on a deep learning model[J]. Well Logging Technology,2019,43(6):578−582.

    [18] 杜春雨,邢强,张晋言,等. 基于注意力约束深度生成网络的测井电成像空白条带填充[J]. 地球物理学进展,2022,37(4):1548−1558. DOI: 10.6038/pg2022FF0320

    DU Chunyu,XING Qiang,ZHANG Jinyan,et al. Blank strips filling for electrical logging images based on attention–constrained deep generative network[J]. Progress in Geophysics,2022,37(4):1548−1558. DOI: 10.6038/pg2022FF0320

    [19] 张浩,司马立强,王亮,等. 基于卷积神经网络的电成像图像空白条带填充方法[J]. 地球物理学进展,2021,36(5):2136−2142.

    ZHANG Hao,SIMA Liqiang,WANG Liang,et al. Blank strip filling method for resistivity imaging image based on convolution neural network[J]. Progress in Geophysics,2021,36(5):2136−2142.

    [20] 陈建华,杨丽丽,赵延静,等. 深度神经网络模型在测井电成像图像处理中的应用[J]. 电子测量技术,2021,44(4):138−143.

    CHEN Jianhua,YANG Lili,ZHAO Yanjing,et al. Application of deep neural network model in processing of electrical logging images[J]. Electronic Measurement Technology,2021,44(4):138−143.

    [21] 袁晓涛,马旭成,肖仕军,等. 融合多层语义特征的测井电成像空白条带填充深度神经网络方法[J]. 西安石油大学学报(自然科学版),2022,37(6):133−139.

    YUAN Xiaotao,MA Xucheng,XIAO Shijun,et al. A deep neural network method for filling blank strips in electrical logging images based on fusion of multi–level semantic features[J]. Journal of Xi’an Shiyou University (Natural Science Edition),2022,37(6):133−139.

    [22]

    WU Yuyan,DENG Rui,LINGHU Song,et al. Method of image restoration of the blank strips of electric imaging logs[J]. Arabian Journal of Geosciences,2022,15(13):1189. DOI: 10.1007/s12517-022-10434-6

    [23] 苏乾潇,乔德新,任义丽,等. 基于傅里叶卷积的电成像测井图像修复[J/OL]. 北京航空航天大学学报,2024:1–11 [2025-02-17]. https://doi.org/10.13700/j.bh.1001-5965.2023.0754

    SU Qianxiao,QIAO Dexin,REN Yili,et al. Inpainting of blank strips in imaging logging images based on Fourier convolution[J/OL]. Journal of Beijing University of Aeronautics and Astronautics,2024:1–11 [2025-02-17]. https://doi.org/10.13700/j.bh.1001-5965.2023.0754.

    [24]

    CRESWELL A,WHITE T,DUMOULIN V,et al. Generative adversarial networks:An overview[J]. IEEE Signal Processing Magazine,2018,35(1):53−65. DOI: 10.1109/MSP.2017.2765202

    [25]

    RADFORD A,METZ L,CHINTALA S. Unsupervised representation learning with deep convolutional generative adversarial networks[C]//4th International Conference on Learning Representations. San Juan,2016.

    [26]

    YEH R A,CHEN Chen,LIM T Y,et al. Semantic image inpainting with deep generative models[C]//2017 IEEE Conference on Computer Vision and Pattern Recognition (CVPR). Honolulu,HI,USA,July 21−26,2017. doi:10.1109/CVPR.2017.728.

    [27]

    ISOLA P,ZHU Junyan,ZHOU Tinghui,et al. Image–to–image translation with conditional adversarial networks[C]//2017 IEEE Conference on Computer Vision and Pattern Recognition (CVPR). Honolulu,HI,USA,July 21−26,2017. doi:10.1109/CVPR.2017.632.

    [28]

    GLOROT X,BORDES A,BENGIO Y. Deep sparse rectifier neural networks[C]//Proceedings of the 14th International Conference on Artificial Intelligence and Statistics. Fort Lauderdale,2010.

    [29]

    MISRA D. Mish:A self regularized non–monotonic activation function[C]//Mish Activation Function. Bhubaneswar,India 2019.

    [30]

    JOHNSON J,ALAHI A,LI Feifei. Perceptual losses for real–time style transfer and super–resolution[C]//Computer Vision–ECCV 2016:14th European Conference. Amsterdam,2016:694−711.

    [31]

    GATYS L,ECKER A,BETHGE M. A neural algorithm of artistic style[J]. Journal of Vision September,2016,16:326. DOI: 10.1167/16.12.326

    [32]

    GUPTA P,SRIVASTAVA P,BHARDWAJ S,et al. A modified PSNR metric based on HVS for quality assessment of color images[C]//2011 International Conference on Communication and Industrial Application. Kolkata,India,December 26−28,2011. doi:10.1109/ICCIndA.2011.6146669.

    [33]

    HORÉA,ZIOU D. Image quality metrics:PSNR vs. SSIM[C]//20th International Conference on Pattern Recognition. Istanbul,Turkey,August 23−26,2010. doi:10.1109/ICPR.2010.579.

图(9)  /  表(3)
计量
  • 文章访问数:  79
  • HTML全文浏览量:  5
  • PDF下载量:  15
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-07-20
  • 修回日期:  2025-02-18

目录

/

返回文章
返回