全球约90%以上的天然气水合物（以下简称水合物）都赋存于海底浅层含泥质沉积物中^[1]，由于储层渗透性差，曾一度被认为是不适宜开采的矿体。2017和2020年，我国在南海北部泥质粉砂型水合物储层中成功实施两轮试采，证明了泥质储层中的水合物也具备一定的技术可采性^[2-3]。然而，平均约3×10⁴ m³的日均产气量离工业气流标准还有很大差距，开采效率亟待提升^[4-5]。因此，加强我国南海北部海域水合物增产研究非常必要。

水合物开采产能主要受到储层等地质因素的制约和开采井配置等工程因素的影响^[6–11]。理论研究和实际试采结果均表明，与垂直井相比，水平井因有效扩大了分解面积而可获得更高产气速率^[12–14]。但超长井段水平井在水合物储层中的应用仍受成本和技术难度的限制^[14]。以多分支井为代表的复杂结构井被认为是实现水合物产能提升的关键^[5,14-16] 。然而，降压开采过程中多分支井不同射孔井段的协同降压会干扰储层中的传热传质^[17]。虽然目前已基本明确水平分支射孔程度对产能的影响规律^[18]，但垂直主井同时射孔及其射孔程度对多分支井开采产能的影响仍不清楚。

此外，储层特性对海域水合物产气、产水能力制约明显^[19]。现已初步探明储层孔隙度、渗透率等条件控制下的产能演化规律^{[20-22,12,23–25]}，地层倾角作为水合物储层非均质性的重要构造特征之一^[13]，已逐渐引起学者们的重视^[24,26]。大量地震调查和实地钻探数据均显示海域水合物通常赋存于海底凹凸起伏的地层中^[27-29]。相关研究表明，地层倾角使得水合物储层的初始温压场分布特征产生差异^[27]，进而影响水合物开采过程中的流体运移和生产行为^[24,30]，尤其是水在储层孔隙中的分布及流动性^[31-32]。此外，倾斜水合物储层在开采过程中也更易发生变形，进而对开采井的稳定性和产能造成负面影响^[30]。学者们认为，在日本倾斜水合物储层的生产过程中，位于储层构造高部位的流体因地层存在倾角发生了破坏性流动，进一步损坏了防砂装置，导致地层严重出砂，从实际证明了地层倾角的存在会对水合物开采行为造成显著影响^[33]。我国南海北部神狐海域泥质水合物储层（如SH2、W11、W17和W19站位）主要分布于峡谷地区^[34]。地层倾角也可能会对该区域泥质水合物储层的开采产能造成影响。然而，目前与此相关的研究开展较少，还未能清晰刻画南海北部具有倾角的泥质水合物储层的开采过程。多分支井降压开采条件下，倾斜泥质水合物储层的长期开采响应行为和产气产水规律仍不明确。

基于已形成的复杂水合物储层建模新技术和复杂结构井建模新方法^{[18, 26]}，本研究以南海北部神狐海域X01站位水合物储层的地质参数和地形参数为依据，建立了代表实际情况的倾斜泥质水合物储层模型，通过TOUGH+HYDRATE水合物产能模拟软件进行数值模拟研究。首先，探讨垂直主井是否射孔及射孔程度对多分支井开采产能的影响和增产效果，确定垂直主井和水平分支的最佳打开长度比值，优选多分支井井型；然后，考虑储层倾角作用下重力分异对水合物开采的影响，揭示位于储层不同构造位置（即构造低部位、倾斜部位和构造高部位）的多分支井的降压开采响应规律，确定优势布井位置。上述研究可深化海域倾斜泥质水合物开采传热传质与产气规律认识，为后续试采和开采工程设计提供一定的理论基础和工程依据。

1.   数值模型

1.1   3D倾斜储层模型建立和网格划分

中新世以来，南海北部神狐海域进入构造沉降期，高沉积速率为水合物的赋存创造了良好的地质条件^[35-36]。2015年，广州海洋地质调查局（GMGS）在神狐海域进行了GMGS3科考，获取了高质量的地球物理测井资料和水合物岩心样品^[34]。调查结果显示，神狐海域海底温度大约为3～6 ℃，地热梯度为4.5～6.7 ℃/100 m。水合物层分布深度约150～400 m，厚度为10～80 m。平均水合物饱和度为20%～40%，最高达76%。神狐海域W17、W19等站位的水合物储层主要位于海底峡谷区的构造高部位或海脊边缘。其中，X01站位含水合物沉积物主要为黏土质粉砂或粉砂质黏土，水合物储层厚度为76 m，位于海脊顶部，见图1。根据多波束探测数据可知，X01站位西北区域地层倾角约为10°～15°。Wang等^[37]对X01站位所在水合物层的空间分布研究结果认为，该区域的水合物储层近似平行于海底展布。

图  1  中国南海珠江口盆地和X01站位所在水合物储层的厚度示意图(A—B)，多分支井结构示意图(C)，基于实际地形参数和地质参数建立的倾斜3D储层模型及多分支井在储层不同构造部位的布设示意图(D)

图中A—B修改自Li等^[2] 和Wang等^[37]，红色虚线框架表示本研究区域。

Figure  1.  Map of Pearl River Mouth basin, South China Sea and thickness of gas hydrate-bearing sediments at site X01(A-B) ; the schematic diagram of a multilateral horizontal well configuration(C) and a inclined 3D hydrate-bearing reservoir model and schematic diagram of the layout of multilateral horizontal wells in different structural positions (D)

A-B are modified from Li et al ^[2] and Wang et al ^[37], red dotted frame is the research area of this study.

下载: 全尺寸图片 幻灯片

本研究采用国际上最认可的TOUGH+HYDRATE水合物产能预测软件进行水合物开采产能模拟^[38-39]。3D储层建模对网格数量和计算时间要求较高，基于前期确立的优势网格划分方案^[18]，本研究采用并行版TOUGH+HYDRATE软件（pT+H）^[40]，建立了与实际海底情况一致的倾角为10°～15°的3D水合物储层模型（1000 m(x)× 1000 m(y)× 312.5 m(z)），见图1D。考虑本研究主要探究实际水合物储层的倾角对开采产能的影响，建模时将水合物层设置为与海底地层平行分布且水合物饱和度一致。除地形参数外，3D倾斜储层的其他参数均参考X01站位的主要储层物性参数（表1）。此外，鉴于X01站位与神狐海域SH2站位的水合物储层物性相似，相对渗透率模型和孔隙水压力模型中经验参数的设定主要参考SH2站位水合物储层建模时选取的经验参数^[41]。

表  1  数值模拟采用的主要储层参数^[34,41]

Table  1.  Main reservoir parameters employed in numerical simulation

参数	数值	参数	数值
储层基础参数
上覆盖层厚度/m	116.5	颗粒密度/(kg∙m−3)	2650
水合物储层厚度/m	76	地温梯度/(℃/100 m)	5.46
下伏地层厚度/m	120	热导率(湿)/(W∙m−1∙℃)	2.917
储层孔隙度/%	34.50	热导率(干)/(W∙m−1∙℃)	1.00
储层绝对渗透率/(10−3μ㎡)	kx = ky = kz = 0.22	压缩系数/Pa−1	1.00×10−8
水合物饱和度/%	34.00	比热容/(J∙kg−1∙℃−1)	1000
X01站位处水深/m	1309.75	孔隙水盐度/‰	3.05
相对渗透率模型(Relative permeability model) [39]
KrA =（SA*）n		KrG =（SG*）nG
SA*=（SA−SirA）/（1−SirA）		SG* =（SG −SirG）/（1−SirA）
SirA	0.50	nG	3.00
SirG	0.05	n	5.00
孔隙水压力模型(Capillary pressure model) [42]
Pcap = −P0[（S*）−1/λ−1]1−λ		S* =（SA−SirA）/（SmxA−SirA）
P0/Pa	105	SmxA	1.00
λ	0.15	SirA	0.50

下载: 导出CSV 

| 显示表格

根据以往经验，本研究在建模时，将上覆盖层和下伏地层厚度分别设定为 116.5 m和120 m，以满足热量和质量的交换^[43-44]。多分支井的共面水平分支垂直于主井布设，位于垂直方向上的储层中部，套管射孔完井方式，各分支打开长度为20 m（图1C）。在x轴和y轴方向上，模型采用不等间距的划分方式。井眼附近传热传质比较剧烈，故在区域离散过程中，对井周网格进行了较为精细的划分：网格尺寸以0.1 m为主。其他储层的网格尺寸随着与开采井距离的增大而增大。在z轴方向上，水合物储层的网格尺寸设定为1.0 m；上覆盖层和下伏地层的网格划分主要遵循Moridis和Reagan^[45]提及的划分规律，即储层的网格尺寸随着与储层距离的增大而规律性增大。储层上、下边界处的网格在z方向的厚度均为0.001 m。

1.2   初始和边界条件

X01站位水合物储层的初始孔隙压力随地层深度的增加而增大。在5.46 ℃/100 m的地温梯度下^[34]，储层初始温度从顶边界到底边界呈线性升高。海水密度平均为1023 kg/m³。3D储层模型的外围边界为绝热无渗流边界，上下边界均为恒温、恒压边界^[46]。本研究涉及的水合物储层初始化细节参考前人研究^[37,47-48]。完成模型初始化后，水合物储层底部温度为15.35 ℃；因地层倾角的原因，水合物层底部压力分布不一，最低压力为15.96 MPa（图1D为初始化后地层温度分布情况）。根据水合物的压力-温度平衡曲线认为水合物可稳定存在。

根据Moridis等^[45]学者的早期研究，可将开采井单元作为一种“伪多孔介质”，并假定井内流动满足达西定律。该“伪多孔介质”的孔隙度值为1，渗透率设定为k_r = k_z = 5.0×10⁻⁹ m²，且毛细压力设定为P_c = 0。在模拟时认为开采井眼稳定、防砂可靠、出水已妥善处理。

1.3   井型结构与布设位置设计

先前研究证实，与一分支井、二分支井和三分支井相比，四分支井对泥质水合物储层的开采能力相对最高^[18]。因此，本研究采用的多分支井为四分支井。如图2所示，为了探讨多分支井垂直主井是否射孔及打开长度对开采产能的影响和增产效果，共设计5个案例。Case 1-1代表多分支井的垂直主井不射孔，Case 1-2— Case 1-5的多分支井的垂直主井打开长度分别为5、10、20 和28 m，即射孔垂直主井长度与分支长度比值（L_v/L_h）分别为0.25、0.5、1和1.4。模拟时，5个案例的多分支井均布设于图1D中的较水平储层，以最大程度控制变量，减少构造条件对井型优选的影响。根据海域实际开采现场试验结果^[44]，设定多分支井采用4.5 MPa进行恒压开采。

图  2  具有不同射孔垂直主井的多分支井结构示意图

Figure  2.  Schematic of multilateral horizontal well configuration with different perforated vertical main wells

下载: 全尺寸图片 幻灯片

确定最佳L_v/L_h后，将优选的多分支井部署在储层的不同结构位置进行降压开采，即将多分支井布设于储层的构造低部位（Case 2-1）、构造倾斜部位（Case 2-2）和构造高部位（Case 2-3）（图1D）。在考虑重力分异对水合物开采的影响下，对比开采产能差异，揭示倾斜水合物储层多分支井降压开采的响应特性和储层渗流规律，确定优势多分支井布设位置。

2.   模拟结果与分析

本文以产气速率、产水速率、气水产出比（即累计气体开采量与累计产水量的比值，水合物开采经济性评价相对指标）、储层水合物分解累计产气量和储层内剩余游离气量来综合对比不同开采方案下的产气产水规律及水合物储层响应特征。

2.1   垂直主井射孔程度对开采产能的影响

图3描述了垂直主井不射孔的多分支井（Case 1-1）和垂直主井射孔的多分支井（Case 1-2 — Case 1-5）降压开采条件下，储层产气量及井中产气速率、产水速率和气水产出比。

图  3  多分支井开采时，产气和产水的情况

A. 储层水合物分解累计产气量和储层内剩余游离气量，B. 井中产气速率和产水速率，C. 井中气水产出比，D. 第720天，井中气水产出比。

Figure  3.  Gas and water production performance when hydrate-bearing reservoirs are produced by multilateral horizontal wells

A. Cumulative gas production volume from hydrate deposits and amount of free gas remaining in the reservoirs, B. gas production rate and water production rate in the production well, C. gas-to-water ratio in the production well, D. gas-to-water ratio in the production well at t = 720 days.

下载: 全尺寸图片 幻灯片

从图3A中可以看出，无论采用何种多分支井开采泥质水合物储层，水合物分解累计产气量随着开采的进行而逐渐增大，但储层经过长期开采后，水合物分解累计产气量的增加幅度逐渐降低。这主要是由于：①随着水合物分解和吸热反应的发生，水平分支周围储层的温度明显降低，在一定程度上抑制了水合物的分解；②随着压降的传递和开采的进行，水合物分解前缘逐渐远离开采井，但储层渗透率极低，水合物分解范围主要集中于井周（图4A）。因此，水合物分解累计产气量的增幅逐渐减弱。从图3A还可发现，无论采用何种多分支井开采泥质水合物储层，储层内剩余游离气量随着开采的进行逐渐增多。但当开采时间超过约2 880 d时，储层的游离气量逐渐降低。此现象表明，多分支井的长期气体产出能力比短期气体产出能力更强。这可能是由于在开采后期，分支井周围储层的渗流能力因水合物分解而得到了提升，进而促进储层内游离气体的产出。如图4B所示，随着开采的进行，垂直主井井周和水平分支下部储层的游离气量逐渐增大，促进更多分解气从多分支井中产出。

图  4  降压开采第1、5和10年，Case 1-5案例储层中水合物饱和度（A）和气体饱和度（B）的分布图

为了更好地呈现结果，主要呈现多分支井附近的物理场图，下同。

Figure  4.  Spatial distributions of hydrate saturation (A) and gas saturation (B) at t = 1 year, t = 5 year and t = 10 year in Case 1-5

for clearly presenting situations, only the spatial distributions around the multilateral horizontal well are shown, same as below .

下载: 全尺寸图片 幻灯片

相比而言，当采用多分支井降压开采泥质水合物储层时，随着多分支井垂直主井打开长度的增加，储层水合物分解累计产气量和储层内剩余游离气量均增加。但对于X01站位这样具有极低绝对渗透率的海域泥质水合物储层而言，开采潜能并不十分理想，若要实施开采作用，还须结合其他辅助增产措施。本研究主要考虑多分支井垂直主井射孔程度和储层构造条件对多分支井开采产能的影响，暂不考虑将多分支井布设在低渗海域泥质水合物储层中的产量是否达到商业开采目的。

不同案例的井中产气产水规律随开采时间的变化较为规则，不同案例间的产能差异变化较小。从图3B中可以看出，由于地层渗透率极低，储层的降压产气效果较弱，即使是开采早期，多分支井的产气速率也比较低，主要为70～100 ST·m³/d。在130 d的模拟周期内，产气速率降低较快。而后，随着降压开采的进行，产气速率逐渐缓慢降低。与产气速率演变规律相反，不同案例的产水速率随着长期开采的进行而逐渐增加。但无论采用何种多分支井进行降压开采，产水速率的增幅并不十分明显，产水速率基本为3.5～4.5 ST·m³/d。图3C呈现了10 年开采时间内多分支井的气水产出比。从图中可以看出，开采早期气水产出比较高，尤其是130 d内的气水产出比（＞20），但随着开采的进行，5个案例的气水产出比均呈现下降的趋势。

与Case 1-2—Case 1-5的产气产水结果相比，Case 1-1井中产气速率、产水速率和气水产出比均最低。随着垂直主井打开长度的增加，井中产气速率和产水速率相应增大。虽然水平分支下部储层的温压分布和水合物分解情况基本不受垂直主井射孔程度的影响（图5A—C），但相同开采时间内，垂直主井井周储层的压降范围和水合物分解范围随着垂直主井打开长度的增加而增大。因此，气、水产出速率随着垂直主井打开长度的增加相应增大（图3B）。但当垂直主井打开长度增大到一定程度后，增强射孔程度带来的产气增幅不如产水增幅明显，导致气水产出比反而随着垂直主井打开长度的增加而降低（图3C）。这可能是由于垂直主井射孔程度增强后，垂直主井与上覆盖层的距离更近，上覆地层中的水在开采过程中更易流入开采井（图5D），水流的增多导致垂直主井产气速率的增幅随其打开长度的增加而降低。

图  5  开采10 年后，Case 1-2，Case 1-4和Case 1-5储层中压力（A）、温度（B）、水合物饱和度（C）和液体饱和度分布图（D）

Figure  5.  Spatial distributions of the pressure (A), temperature (B), hydrate saturation (C), and water saturation (D) at t = 10 year in Cases 1-2, 1-4 and 1-5

下载: 全尺寸图片 幻灯片

通过进一步对比相同开采时间（以开采第720天为例，图3D）内垂直主井射孔程度不同的多分支井的气水产出比可发现，当L_v/L_h＜0.5时，气水产出比随L_v/L_h的增大而增大；当0.5＜L_v/L_h＜1.0时，L_v/L_h对气水产出比影响不明显，气水产出比值均较高，当L_v/L_h = 0.75时，值最大；当L_v/L_h＞1.0时，随L_v/L_h的增大，气水产出比反而降低。因此，开采低渗泥质水合物储层时，L_v/L_h比值为0.5～1.0时，最利于提高多分支井的气水产出。

2.2   布设于不同构造位置的多分支井的产能差异

如上所述，随着垂直主井打开长度的增加，多分支井的产气速率逐渐增大。但当L_v/L_h为0.75时，多分支井的气水产出比最高。因此在本小节研究中，以垂直主井打开长度为15 m、水平分支打开长度为20 m的多分支井为例，对比布设于不同构造位置的多分支井的产能差异。

图6呈现了上述多分支井布设在泥质水合物储层不同构造部位时的产气产水情况。从图中可见，经过相同开采时间，不同案例中开采井的产气速率、产水速率和气水产出比随开采时间的变化较为规则。相比而言，Case 2-1和Case 2-2中的产气速率和产水速率均明显高于Case 2-3中的产气速率和产水速率。其主要原因是，在Case 2-3中水合物储层的埋藏深度较浅，水合物储层的初始温度和压力相对较低。在相同的开采条件下，Case 2-3中的水合物受到更弱的压降作用而发生分解反应，水合物分解范围更小（图7）。此外，虽然在开采过程中Case 2-2储层的压降比Case 2-3储层的压降更明显，但经过相同开采时间，布设于储层构造高部位的多分支井的井周储层温压明显更低（图7），在一定程度上抑制了水合物的进一步分解，故Case 2-3在开采时呈现出相对较低的气、水产出速率。在Case 2-1和Case 2-2两个案例中，尽管Case 2-2的储层初始温压较Case 2-1的储层初始温压条件低，但在开采过程中流体流动随时都受到重力分异作用的影响，导致Case 2-1和Case 2-2中的产气速率和产水速率相差不显著（图6）。

图  6  多分支井布设在不同构造位置的产气产水情况

Figure  6.  Gas and water production performance of multilateral horizontal wells deployed in different structural positions

下载: 全尺寸图片 幻灯片

图  7  开采5 年后，Case 2-2和Case 2-3的储层物理场分布图

Figure  7.  Spatial distribution of physical properties in Cases 2-2 and 2-3 at t = 5 year

下载: 全尺寸图片 幻灯片

虽然多分支井布设于储层倾斜处和构造低部位时的产气、产水效果基本一致，但多分支井位于泥质水合物储层不同构造位置时，气水产出比存在差异。相比而言，多分支井布设于泥质水合物储层的构造高部位（Case 2-3）时，气水产出比最大（图6）。随着多分支井布设储层埋深的增加，气水产出比值逐渐降低。对比Case 2-1和Case 2-2的气水产出比可发现：①对于短期开采而言，将多分支井布设于储层倾斜处（Case 2-2），气水产出比较高；②对于长期开采而言，Case 2-1和Case 2-2的气水产出比逐渐趋于一致。水平井开采倾斜泥质水合物储层的研究指出，水合物分解易引起开采井发生变形、井周沉积物发生位移，进而影响开采安全和采收率^[30]。采用多分支井长期开采具有倾角的泥质水合物储层时，采气安全也极可能受到倾角的影响。相比于倾斜地层，水平地层对实际水合物开采过程中井壁稳定性和储层稳定更有利，而产出水也可以依赖设备进行适当处理，故认为将多分支井部署在位于构造低部位的较水平泥质水合物储层可能更利于提高采收率。值得注意的是，对于渗透率极低的泥质水合物储层而言，仅采用降压法进行开采，多分支井的产气产水情况仍不理想，导致地层倾角对多分支井降压开采造成的产能差异的影响有限。若在海域低渗泥质储层中开展实际开采，势必会结合其他辅助增产措施，如储层改造技术。在联合开采技术的应用下，储层倾角的存在对开采产能造成的影响会更明显，尤其对于渗透率较大的泥质水合物储层，更应将多分支井部署在位于构造低部位的水平水合物储层中。此外，本研究仅考虑了多分支井射孔程度和布设位置对倾斜泥质水合物储层开采产能的影响，对于储层变形和非均质性对开采产能的影响在后续的研究中需进一步量化考虑。

3.   结论

（1）相比于仅水平分支射孔的多分支井，水平分支和垂直主井同时射孔的多分支井更利于水合物分解产气。但垂直主井的打开长度不宜过长，垂直主井与水平分支打开长度比值介于0.5～1.0时，最利于提高多分支井的气水产出比。

（2）海域泥质水合物储层的地层倾角影响多分支井的开采产能。将多分支井布设在储层构造低部位的水平位置更利于水合物的长期开采。但对于X01站位此类具有极低绝对渗透率的泥质水合物储层而言，开采潜能并不十分理想；若要实施开采作用，在将多分支井布设在优势构造储层进行开采的基础上，还须结合其他辅助增产措施，如井网开采模式和储层改造技术。