我国化石能源呈现多煤少油少气的特征，煤炭在现在和未来一段时间内都是我国能源消费的主体之一^[1]。随着我国煤炭开发重心逐步向西部干旱半干旱地区转移，部分矿区呈现富煤少水特征，但也有多个矿区矿井水量大、含盐量偏高，因此，探索矿井水的处理、减排方法十分必要。

保水采煤自20世纪末提出后，在保护我国西北干旱半干旱地区水资源完整方面做出突出贡献^[2]。在煤炭开采过程中，煤层顶板产生失稳垮落，垮落岩体上形成导水裂隙；当煤层距离含水层近，裂隙导通上方直接充水含水层时，水沿裂隙向下进入开采工作面、巷道，被矿井排出地表，使含水层失水；当失水量超过含水层自稳定限度后会形成当地水资源不可逆缺失^[3-5]。在干旱半干旱缺水矿区，为保护受煤层开采影响的浅部含水层中地下水资源的可利用性，邵飞燕^[6]提出将浅部含水层地下水向煤层底板以下含水层中进行“转移存储”的技术思路。武强等^[7-8]基于矿山开采过程中矿井水量大但利用程度不高的情况提出“煤-水”双资源型矿井开采技术方法，以梧桐庄矿突水治理为例，对矿井水进行控制、处理、利用、回灌与生态环保“五位一体”的优化结合。我国要求水资源短缺矿区矿井水利用率要达到100%，高盐矿井水减量减排是实现矿井水零排放的重要途径。目前，高盐矿井水处理多集中在矿井水中盐分消除^[9]，达标后作为矿区生产、生活和生态用水。此外，地下水转移存储也是一种有效处理高盐矿井水的手段^[10]。国外较先开展地下水转移存储研究^[11-12]，我国主要集中在近地表地层和采空区^[13-14]，深部转移存储进行过CO₂封存相关的研究与实践，如神华CCS示范工程^[15]，其主要注入层位为刘家沟组、石千峰组、山西组、石盒子组和马家沟组，工程结果表明，鄂尔多斯盆地三叠系刘家沟组是地质封存CO₂的良好储层，至今未发生泄露；说明刘家沟组具有一定的转移存储潜力。

为保障鄂尔多斯盆地X矿安全高效开采煤炭资源，实现高盐矿井水减量减排，开展高盐矿井水异位转移存储储层的水文地质结构特征和介质特征研究，通过对比分析研究区宝塔山砂岩与刘家沟组砂岩的介质特征与水动力学特征，优选目的层，并对矿井水回灌后的水动力影响进行了模拟，以期指导矿井水回灌工程的实际运行。

1   研究区概况

1.1   研究区自然条件

研究区X矿位于鄂尔多斯盆地，属大陆性沙漠气候。含煤地层为侏罗系延安组，含煤5组(2、3、4、5、6煤组)，最低可开采煤层为6^-2中煤层。现开采煤层为3^-1煤，矿井涌水量已达1 900 m³/h，矿化度最高为2 g/L。该矿现有矿井水处理方法有深度处理与中转水池蒸馏结晶法、化工用水、矸石处理用水和离层注浆用水，但以上水处理方式已趋于饱和，难以满足开采需要。

1.2   转移存储目的层特征及选取原则

转移存储目的层选取应满足以下2点准则：① 从含水层介质特征方面，要求目的层必须分布稳定，厚度较大，具有一定程度的孔隙、裂隙发育，以提供储水空间；②从水文地质结构方面，位于煤层底板以下较大深度，以免受开采顶底板破坏带影响，且不能造成底板突水危险性增大，最好为相对封闭的独立水文地质单元。

3^-1煤直接充水含水层为顶板直罗组砂岩孔隙、裂隙含水层和延安组砂岩孔隙、裂隙含水层。延安组下段发育宝塔山中粗砂岩，已有研究表明，宝塔山砂岩孔隙发育，渗透性好，具有相当大的储水空间^[16-17]；三叠系刘家沟组主要由细砂岩、砂质泥岩组成，研究区内刘家沟组距离3^-1煤1 200 m以上；据文献[18]，鄂尔多斯盆地刘家沟组砂岩段裂隙发育，钻井液周期性漏失，堵漏27次，漏失量最大达4 621 m³，这与CCS示范工程模拟结果^[15]相吻合。宝塔山砂岩、刘家沟组砂岩与最低可开采煤层6^-2中煤层距离如图 1所示。

图  1  开采煤层下地层柱状图

Figure  1.  Columnar diagram of strata under mining coal seam

下载: 全尺寸图片 幻灯片

依据转移存储目的层选取准则，分析研究区含水层介质特征和水文地质条件，选择煤层下储水空间用以转移存储矿井水，钻孔揭露延安组底部宝塔山砂岩厚47.1 m、三叠系刘家沟组地层厚416.4 m，储水空间大，提出宝塔山砂岩段和深部刘家沟组砂岩段作为转移存储层备选目标。

2   转移存储层介质特性分析

转移存储层的储水空间直接决定了其原始状态下的可转移存储量，转移存储层砂岩的介质特征直接影响转移存储层的储水空间大小^[19]。研究区宝塔山砂岩为粗砂岩和含砾粗砂岩，距离6^-2中煤近。刘家沟组成岩时代较煤层早，主要为砂岩。收集宝塔山砂岩和刘家沟组岩心样本13组，测试孔隙率及孔径小于100 nm、100~10 000 nm和大于10 000 nm孔隙占比(表 1)，获取研究区转移存储层的含水介质条件。

表  1  目标转移存储层孔隙分布及黏土含量

Table  1.  Pore distribution and clay content of the target transfer storage layer

存储层	孔隙率/%	< 100 nm孔隙占比/%	100~10 000 nm孔隙占比/%	>10 000 nm孔隙占比/%	黏土质量分数/%
宝塔山砂岩	6.57~19.89/14.53	5.75~17.55/11.31	66.85~83.80/74.87	10.45~15.60/13.82	12.70~25.90/17.13
刘家沟组砂岩	4.18~7.49/5.50	14.40~42.54/25.78	20.48~65.02/36.89	20.98~55.86/37.33	11.20~23.50/17.45
注：数据表示最小~最大值/平均值。

下载: 导出CSV 

| 显示表格

已有研究表明，含水介质孔隙多少和大小都会对地下水运动造成显著影响(表 2)^[20]。测试结果表明，宝塔山砂岩孔隙率均值为14.53%，重力水可以运动的孔隙占比平均为88.69%，最低为82.45%，其中100~10 000 nm孔隙占比平均为74.87%，表明宝塔山砂岩孔隙多位于这一级别，重力水可在较高水头下运动，作为储层转移存储矿井水潜力巨大。刘家沟组砂岩孔隙率均值为5.50%，重力水可以运动的孔隙占比平均为74.22%，较宝塔山砂岩总孔隙率小，重力水可以运动的孔径占比小，相当一部分孔隙充满结合水，表明在原始状态下，宝塔山砂岩在单位空间内转移存储矿井水潜力优于刘家沟组地层。

表  2  微孔隙结构特征^[20]

Table  2.  Micropore structure characteristics^[20]

渗透性等级	孔径/nm	孔隙分布
孔隙充满结合水	< 100	粒间和溶蚀孔隙
重力水在较高水头下运动	100~10 000	粒间、晶间孔隙
重力水一定水头下运动，毛细上升快且不高	10 000~106	颗粒内和粒间孔隙
重力水可以自由运动	>106	晶粒间孔隙

下载: 导出CSV 

| 显示表格

在长期转移存储过程中，转移存储层的抗压能力也会影响其储水潜力，这是由于转移存储过程中转移存储层水位上升导致回注井孔周围岩层受压产生不同程度的破裂，部分新裂隙产生、部分原有微裂隙扩张形成裂隙网络，进一步增大储水潜力。对宝塔山砂岩4组样本和刘家沟组砂岩13组样本进行单轴抗压强度测试，结果如图 2所示；黏土含量等物性参数测试结果见表 1。

图  2  单轴抗压强度测试结果

Figure  2.  Uniaxial compressive strength test results

下载: 全尺寸图片 幻灯片

结果表明，原始状态下宝塔山砂岩饱和抗压强度低于刘家沟组砂岩；在长期转移存储饱和含水条件下，转移存储层抗压强度均降低，宝塔山砂岩抗压强度从32.1~60.2 MPa降低到22.6~44.3 MPa；刘家沟组岩层抗压强度从59.7~114.0 MPa降低到30.0~66.1 MPa，刘家沟组降低显著。对刘家沟组样本进行X射线衍射得到黏土矿物平均质量分数17.45%，最高为23.50%。已有研究表明黏土含量的存在会弱化砂岩强度，其中含少量或大量泥质的砂岩，其单轴抗压强度较不含泥质砂岩分别折减了22%和4%^[21]。刘家沟组砂岩为泥质胶结，抗压强度较低，在饱和含水与受压双重耦合下易发生开裂，增大原岩渗透性。

3   转移存储层盖层及水化学特征

根据钻孔取心及物探分析，延安组中含2、3、4、5、6煤组，宝塔山砂岩距6^-2中煤仅8~10 m，受采动影响，煤层开采扰动裂隙现已导通直罗组含水层，导致宝塔山砂岩无稳定盖层，某些开采侏罗纪煤田下组煤的煤矿，宝塔山砂岩甚至作为充水含水层考虑并进行疏放水可行性研究，证明煤层开采存在垂直裂隙导通含水层，宝塔山砂岩作为转移存储层时，矿井水会产生沿裂隙扩散逃逸的风险。

刘家沟组上段与和尚沟组主要以泥岩、砂质泥岩和粉砂岩为主，厚度达180 m。鄂尔多斯盆地属于我国典型的克拉通盆地，地壳稳定，地震活动微弱，研究区内和尚沟组盖层和刘家沟组上段不发育垂向裂隙，具有十分良好的封盖作用。

刘家沟组砂岩水赋存于埋藏深、温度高的低孔隙率砂岩中，补径排条件差，TDS达到65.11 g/L，其中主要离子Cl^-质量浓度为39.74 g/L，SO₄^2-为2.08 g/L，Ca²⁺为14.51 g/L，Na⁺为7.82 g/L，表现出极高矿化度、极高盐分的特征，pH呈弱酸性。整体上，矿井水水质优于刘家沟组原生水水质，整体可呈现良水转移存储(表 3)。

表  3  转移存储层矿井水水质

Table  3.  Mine water quality of transfer storage layer

水类型	TDS/(mg·L-1)	PH值	水化学类型
宝塔山砂岩水	653~1 452	11.40~11.70	CO3-Na
刘家沟组砂岩水	65 111	5.35	Cl-Ca· Na
矿井水	1 134~2 411	6.95~8.33	SO4-Na

下载: 导出CSV 

| 显示表格

4   转移存储层水动力特征分析

4.1   宝塔山砂岩转移存储可行性分析

根据区域内宝塔山砂岩放水试验显示，宝塔山砂岩含水层水位1 180.10~1 200.33 m，渗透系数K为0.105 7~2.024 7 m/d。根据《煤矿防治水细则》中突水系数法安全水头值计算公式，对宝塔山砂岩段进行允许最大水位计算：

式中：p_s为底板隔水层安全水头值，MPa；M为底板隔水层厚度，取8~10 m；T_s为临界突水系数，取0.1 MPa/m。

计算得到最低可开采煤层6^-2中煤层可承受宝塔山砂岩含水层最大安全水压为0.8~1.0 MPa，研究区宝塔山砂岩顶界高程440~470 m，因此，允许水位为520~570 m，在不影响未来开采6^-2中煤条件下，宝塔山砂岩含水层还需疏降水位，因此，宝塔山砂岩段可存储量十分有限。在现开采阶段，宝塔山砂岩段不考虑作为转移存储层存储高盐矿井水，可作为该矿绿色闭坑后矿井水的转移存储层，形成分布式地下水库等^[22]。

4.2   刘家沟组水文地质试验

针对刘家沟组进行自然水位恢复试验、水力压裂和多期次注水试验，得到研究区域内刘家沟组的渗透性和富水性。2019年12月21日至12月29日对MC-1试验回注孔进行自然水位恢复观测(图 3)，洗井阶段，刘家沟组无水反出；为准确检测水位变化，在洗井结束后注水至水位-213.10 m开始观测，得到水位恢复速率为5 cm/h，推测刘家沟组砂岩含水层静水位约-100 m。

图  3  MC-1孔水位自然恢复曲线

Figure  3.  Natural recovery curve of water level in hole MC-1

下载: 全尺寸图片 幻灯片

利用Aquifer Test软件微水试验中Bouwer & Rice方法计算原岩状态下刘家沟组的渗透系数，具体计算公式如下：

式中：r_c为井管半径；y₀=y(t₀)为初始时刻水位；y_t=y(t)为t时刻后水位；r_w为过滤段钻孔有效半径；L为筛管长度；R/r_w为水流长度；R为影响半径；t为时间。

利用MC-1孔近似自然水位恢复数据，代入式(2)计算得到刘家沟组的渗透系数K为5.31×10^-6 m/d，表明刘家沟组原始状态下为弱渗贫水含水层，总体呈现弱地下水循环交替；而刘家沟组地下水TDS极高，表明在原岩条件下刘家沟组地下水基本处于停滞状态。

为满足矿区转移存储水量的需要，对刘家沟组开展水力压裂试验进行局部增透，共进行6次注水试验。根据注水试验观测数据，刘家沟组在后期稳定阶段满足地下水向承压水完整井的稳定运动条件，即符合圆岛模型假设(广义)：无限含水层，产状水平、等厚、均质、各向同性；长时间抽水后会出现似稳定状态，取R为影响半径。则有：

式中：M_h为含水层厚度，m，取309 m；s_w为降深，m；Q为流量，m³/h；r_w取53.5 mm。

除第2次注水试验水压波动外，依据5次注水试验稳定阶段进行水文地质参数计算(表 4)，得到刘家沟组经水力压裂后渗透系数为0.008 14~0.015 27 m/d。相较刘家沟组原始渗透系数增大几个数量级，表明水力压裂局部增透可明显改善刘家沟组水文地质条件。

表  4  注水试验求参成果

Table  4.  Results of water injection experiments

序号	孔径/mm	水位差/m	渗透系数/(m·d-1)
1	107	930	0.011 93
3		700	0.010 62
4		710	0.008 14
5		750	0.010 51
6		760	0.015 27

下载: 导出CSV 

| 显示表格

4.3   刘家沟组长期转移存储水动力分析

前期注水试验初步掌握刘家沟组的水文地质参数。在长期转移存储矿井水过程中，刘家沟组的水文地质条件可能会发生变化。为获取高强度、持续性注水条件下刘家沟组的流场演化及水压变化规律，利用2020年3月8日至4月30日矿井水回注数据(图 4)，分析刘家沟组地层变化情况。

图  4  注水井井口水压和注水量变化曲线

Figure  4.  Change curves of wellhead water pressure and water injection volume of injection well

下载: 全尺寸图片 幻灯片

在水力压裂试验后进行试回灌，并获得相关数据，利用式(3)再次进行参数计算，得到试回灌阶段渗透系数为0.009 11~0.016 93 m/d，与压裂阶段渗透系数差别较小，证明回注地层压裂后保持稳定，刘家沟组压裂后转移存储矿井水的潜力较大。

分析结果表示，在为期2个月的回灌过程中，转移存储层刘家沟组的渗透系数基本稳定，说明介质特征也基本稳定，压裂后刘家沟组不再产生二次裂隙扩展，介质条件稳定，压裂后也有较好的转移存储空间，顶底板盖层发育较好，未发生泄露。

基于以上分析，利用MODFLOW模拟长期转移存储下刘家沟组的流场变化情况，用于定量分析和预测长期转移存储条件下储层的水动力场变化情况。研究区刘家沟组在鄂尔多斯盆地区域发育稳定连续，构造不发育，在实际数值模拟过程中，取100R(R为影响半径)作为模拟边界(零流量边界)。根据现有资料，对刘家沟组进行有限单元网格剖分，并对影响半径内进行网格加密，共计剖分网格数10 147个，模拟步长为1 d。因回注孔MC-1无观测孔，采用石油井田在周边地区的水力压裂影响范围经验值700 m进行计算。利用2个月实际资料进行参数率定，设置在压裂影响范围内渗透系数为0.011 m/d，压裂影响范围外为5.31×10^-6 m/d，初始水位-100 m，回灌量2 400 m³/d，模拟时间60 d。得到回注60 d后流场变化情况(图 5)，在设置回灌量下，井口水位达到773 m，图 4中实际观测压力7.5 MPa，表明模型参数设置较为合理。

图  5  模拟60、180和360 d后刘家沟组流场变化

Figure  5.  Flow field changes in Liujiagou Formation after 60, 180 and 360 days

下载: 全尺寸图片 幻灯片

模拟结果表明，模拟时间为60 d时，转移存储水体扩散半径达到274 m；模拟时间为180 d时，转移存储水体扩散半径达到472 m；模拟时间360 d时，转移存储水体扩散半径达到665 m。在单口井转移存储1 a情况下，转移存储水体扩散半径接近压裂半径，表明刘家沟组的转移存储矿井水前景较好，可作为该矿处理矿井水，实现矿井水减量减排，转移存储水资源的一种新途径。

5   结论

a. 转移存储目的层选取的两点准则为：目的层必须分布稳定，厚度较大，具有一定的孔隙、裂隙发育程度，以提供储水空间；位于煤层底板以下较大深度，以免受开采扰动，最好为相对封闭的独立水文地质单元。基于以上准则，在鄂尔多斯盆地X矿内选择延安组宝塔山砂岩和刘家沟组砂岩作为备选转移存储层。

b. 宝塔山砂岩段厚度47.1 m，平均孔隙率为14.53%，其中重力水可移动的孔隙占比88.69%，但距离煤层近，无稳定盖层，转移存储高盐矿井水会产生垂向扩散逃逸的风险；刘家沟组总厚416.4 m，砂岩段厚度309 m，平均孔隙率为5.50%，其中重力水可运动的孔隙占比74.22%，原生条件下孔隙发育较差，但存在水力压裂的可能，对煤层开采无影响，更适合用于矿井水的转移存储。

c. 刘家沟组水力压裂局部增透后渗透系数增大到0.008 14~0.015 27 m/d，在长期注水试验条件下渗透系数不发生明显变化，保证压裂后介质特征基本稳定。长期转移存储数值模拟结果表明，刘家沟组具有较为良好的转移存储矿井水前景。

d. 宝塔山砂岩可考虑作为该矿闭坑后的矿井水转移存储层，与已采煤工作面形成采空区+宝塔山砂岩含水层联合分布式地下水库，合理调蓄当地水资源。