Numerical simulation for effect of discrete fracture network on shale gas productivity of horizontal well
-
摘要: 页岩作为典型的非常规储层,基质孔隙小,渗透率极低,水平井多级水力压裂为其商业开发的主要手段。准确模拟页岩气产能,应同时考虑水力裂隙和天然裂隙的渗流。基于离散裂隙模型和等效连续模型建立页岩气渗流数学模型,利用有限元分析方法进行数值求解,研究不同走向裂隙组对页岩气井产能的影响。研究认为,页岩基质为气体的生产提供了主要气源,天然裂隙作为渗流的主要通道,将气体输送到水力裂缝,进而到达井筒。模拟结果表征,离散裂隙的渗流特征对于页岩气井的产能有重要影响。根据页岩储层的天然裂隙走向,可以优化相应的水平井方位。对于二维离散裂隙网络模型,水平井沿着2个裂隙组夹角的平分线更有利于生产。Abstract: As a typical unconventional reservoir, shale gas reservoirs have well-developed natural micro-fractures, in which the pores of matrix are small and the permeability is extremely low. Multi-staged hydraulic fracturing for horizontal well has become a key technology for commercial development. Therefore, in order to accurately simulate the shale gas production, hydraulic fractures and natural fractures should be considered simultaneously. In this paper, the mathematic model is established based on discrete fracture model (DFM) and equivalent continuum model (ECM). The model is solved by Finite Element Method (FEM). According to the model, the effect of different fracture sets with different fracture orientation is investigated. The results show that matrix provides main gas resource for gas production and natural fracture as main passage undertakes a task of transporting gas into hydraulic fracture, and then to wellbore. The characteristic of gas flow in discrete fracture network has an important impact on gas production. According to orientation of natural fracture, corresponding horizontal wellbore azimuth can be optimized. For 2D discrete fracture model, the horizontal well along the angle bisector of two fracture groups is more advantageous to gas production.
-
习近平总书记在第75届联合国大会上提出“中国二氧化碳排放力争于2030年前达到峰值,争取在2060年前实现碳中和”[1]。据国际能源机构(IEA)统计,2021年,与能源相关CO2排放量达到363亿t,煤炭使用产生的碳排放占全球CO2排放增量的40%以上,达到153亿t,均创历史新高。中国CO2排放量达到119亿t,约占全球碳排放总量的33%[2]。2021年10月,国务院印发《2030年前碳达峰行动方案》,提出“构建新能源占比逐渐提高的新型电力系统,推动清洁电力资源大范围优化配置”。当前我国正构建以光、风、水、核为主体的新能源电力系统,并推动电力行业从高碳向低碳、从以化石能源为主向以清洁能源为主的发展方向[3-4]。在新能源电力系统发展过程中,风、光等可再生能源存在强间歇性、随机性与波动性等诸多问题[5],迫切需要建设储能电站,保障电力系统的稳定性与安全性。在诸多储能方式中,抽水蓄能电站是促进新能源大规模利用、实现多能互补、打造安全可靠与经济灵活电力系统的最可靠方法[6]。
抽水蓄能电站具有技术成熟、容量大、能源效率高、使用寿命长等优点,然而其建设受到地形、土地可利用性和较高的资金成本的限制,同时易对周边生态环境产生负面影响,因此,抽水蓄能电站选址成为亟待解决的重要科学问题。针对矿产资源开发过程中形成的大小、形态不一的巷道、硐室及采空区等地下空间[7],袁亮、谢和平、顾大钊等院士和专家提出了废弃矿井抽水蓄能电站利用的战略构想[8]。废弃矿井巷道具有空间大、结构简单等特点,进行简易加固密封后即可使用,是建设抽水蓄能电站地下水库的理想场所[9],为抽水蓄能电站选址提供了新的方向。利用废弃矿井建设抽水蓄能电站可以减少土地占用和景观破坏,降低电站建设投资[10],缩短建库周期[11],具有良好的经济、社会和环境效益。
伴随煤炭供给侧结构性改革的不断深化,黄河流域九省区均推行煤炭行业化解过剩产能政策,关闭了数量众多的资源枯竭型或落后产能煤矿,形成了大量废弃矿井地下空间。为实现空间资源高效利用,笔者开展废弃矿井抽水蓄能利用潜力评估,提出一种废弃矿井抽水蓄能选址评价指标体系,并进行典型废弃矿井抽水蓄能电站实例分析,以2016—2020年黄河流域九省区关闭煤矿资料为基础,对黄河流域废弃矿井抽水蓄能电站利用潜力进行动静态效益定量计算分析,以期为黄河流域九省区废弃矿井抽水蓄能电站建设提供参考。
1 废弃矿井抽水蓄能电站发展现状与工作原理
20世纪初,美国麻省注册了一个结合风车的半地下、半地表蓄能装置[12],但该前瞻性想法在之后几十年并未实现。20世纪60年代起,以美国、德国为主的许多国家开始进行地下抽水蓄能电站的基础研究[13-15]。1960年,R. D. Harza[16]提出可使用地下或露天矿山建设抽水蓄能电站。美国提出建设Mount Hope抽水蓄能电站项目[17],其上水库是在地面开挖的人工水池,下水库则是利用地下约760 m深处废弃的矿井空间,但该项目并未成功利用废弃矿井作为下水库。随着理论与技术的不断完善,废弃矿井抽水蓄能电站逐步由理论走向工程实践。奥地利在阿尔卑斯山的Nassfeld成功建造首个真正意义的半地下抽水蓄能电站,上水库为2 km的人工隧道群,下水库为天然湖泊[18],为废弃矿井抽水蓄能电站工程建设奠定了基础。澳大利亚计划用露天金矿建设抽水蓄能电站,储能规模预计为250 MW[19],并于2021年举行Kidston抽水蓄能发电站破土动工典礼。南非规划利用Fast West Rand区废弃的深井金矿建设一个大型全地下抽水蓄能电站[20],目前该项目正在做大量详细的数据分析及可行性分析。德国正在利用废弃的Prosper-Haniel煤矿改建为半地下抽水蓄能电站,若项目实施成功将成为世界第一个废弃煤矿抽水蓄能电站,同时下萨克森州能源研究中心也在规划将废弃金属矿巷道改建为全地下的抽水蓄能电站[21-22],当前进行试验型抽水蓄能电站建设。国外废弃矿井抽水蓄能电站工程设计参数见表1。
表 1 国外部分废弃矿井抽水蓄能电站工程设计参数[23]Table 1. Engineering design parameters of some pumped storage power plants using abandoned mines abroad[23]电站名称 类型 水库巷道长/m 水头/m 库容/103 m3 功率/MW 储能容量/(MW·h) 西班牙Asturian发电站 半地下 5 700 300~600 170 23.52 141 南非Fast West Rand发电站 全地下 67 000 1 200/1 500 1 000 抽955/发1 230 6 800 德国Prosper-Haniel发电站 半地下 15 500 560 600 200 820 德国Grund ore mine发电站 全地下 25 000 700 240~260 100 400 当前国内学者在废弃矿井建设抽水蓄能电站方面研究较少。2010年,神东大柳塔矿建设煤矿地下分布式水库示范工程[24],顾大钊院士团队建立了煤矿地下水库理论框架和技术体系[25],为废弃矿井抽水蓄能地下水库建设提供了借鉴。2015年,谢和平院士团队全面阐述了煤矿井下抽水蓄能发电的概念与技术[26]。2021年中煤能源大屯煤电公司对2022年底即将闭坑的江苏沛县龙东煤矿抽水储能电站建设项目开展前期可行性研究。2022年山东省地矿局八〇一队与中国华电集团山东公司、山东省科学院合作,在淄博市启动废弃矿井抽水蓄能发电多能互补能源综合体项目,计划利用闭坑煤矿建设分布式矿坑水抽水蓄能电站及供暖能源站,开启我国废弃矿井抽水蓄能电站工程应用新篇章。
废弃矿井抽水蓄能电站工作原理与常规抽水蓄能原理一致,其不同之处在于水库建设采用矿井原有的塌陷区或巷道群。当电网在负荷低谷期时,利用富余电能将下水库水抽至上水库,当电网处于负荷高峰期时,将上水库水放至下水库发电。依据国内外现有废弃矿井抽水蓄能电站工程实践规划与工作原理,黄河流域九省区废弃矿井建设抽水蓄能电站构建原理如图1所示。本文提出黄河流域九省区废弃矿井抽水蓄能电站建设的2种模式:第一种为半地下式,上水库利用塌陷区,下水库利用巷道群(图1a);第二种为全地下式,上下水库为处于不同水平的回风巷道群与运输巷道群(图1b)。废弃矿井抽水蓄能电站模型主要由上水库(矿山巷道群/塌陷区)、下水库(矿山下水平巷道群/上水平巷道群),及在矿井主/副井中布置的引水系统与地下厂房等组成。
不同模式的废弃矿井抽水蓄能电站各有其优缺点,半地下式抽水蓄能电站优点在于地下空间需求较少和改建要求相对较低,缺点是上水库的建造需利用沉陷区或工程开掘,且存在水的蒸发损失;全地下抽水蓄能电站的优点是减少了上水库的土地占用、减轻了水资源的蒸发损失,缺点是对地下水库要求较高,既要保持上下水库的有效库容大致相当,又需避免上下水库串通。黄河流域矿区生态以绿色、稳定和可持续发展为导向,许多矿区采用绿色采矿方法防止形成塌陷区,同时加强了矿区土地复垦与生态恢复工作,该区蒸发强烈。据此,全地下式模式更适合黄河流域九省区废弃矿井抽水蓄能电站的建设。
2 黄河流域九省区电力与矿产资源概况
2.1 电力情况
2016—2020年,黄河流域九省区总发电装机容量与发电量全国占比均在36%以上,发电装机容量增速与发电量变化如图2所示。从图2a来看,黄河流域呈现用电量持续增长的特点,但火电依旧占据主体地位。由图2b可知,5年间黄河流域九省区发电装机容量增速均高于全国平均水平,反映出黄河流域九省区在全国能源发展中的重要战略地位。其中,火、电装机容量增速处于全国领先水平,风、光等可再生能源装机容量增速较快,但水电装机容量增速仅在2017年高于全国平均值,而最近2年增速不到1%,水电发展陷入桎梏阶段。在碳达峰碳中和(“双碳”)目标下,火电企业纷纷转型为新能源电力运营商,黄河流域九省区电力供需结构将会被打破,必将构建以风、光、水电等清洁能源为主体的电力系统。黄河流域九省区加快风光等可再生能源开发和抽水蓄能电站建设势在必行。
2.2 矿产资源概况
目前我国主要的能源矿产是煤炭、石油与天然气,其中,黄河流域九省区的煤炭储量全国占比70.9%,石油占比34.6%,天然气64.2%,是我国能源生产的重要战略地带(表2)。
表 2 2020年黄河流域九省区主要能源矿产储量情况Table 2. Main energy mineral reserves in the nine provinces of Yellow River Basin in 2020省/自治区 煤炭/亿t 石油/万t 天然气/亿m3 青海 2.26 8 251.85 1 055.32 四川 26.66 555.40 15 274.98 甘肃 15.31 39 560.97 588.00 宁夏 35.01 4 670.51 280.67 内蒙古 194.47 6 676.91 10 123.53 陕西 293.90 36 812.83 11 096.45 山西 507.25 0 1 402.04 河南 33.65 3 022.24 62.82 山东 41.32 25 493.92 343.52 合计 1 149.83 125 044.63 40 227.33 全国 1 622.88 361 885.75 62 665.78 煤炭作为黄河流域九省区最重要的一次能源,在区域经济发展中起到了关键作用。随着资源枯竭和国家政策等原因,大量矿井被关闭或废弃[27]。据统计,2016—2020年黄河流域九省区关闭煤矿数量达到1 314个,其中,小型矿井数量占79%,中型矿井占17%,大型矿井占4%,总规模达38 896万t/a(表3)。黄河流域九省区废弃煤矿具有体量大、分布广泛等特征,为建造废弃矿井抽水蓄能电站创造了良好的基础条件。
表 3 2016—2020年黄河流域九省区关闭煤矿情况Table 3. Closed coal mines in the nine provinces of Yellow River Basin from 2016 to 2020省/自治区 不同产能规模/(万t·a−1) 合计 (0,30] (30,120) [120,500] 数量 规模 数量 规模 数量 规模 数量 规模 青海 19 168 0 0 1 120 20 288 四川 447 4 925 11 840 1 150 459 5 915 甘肃 92 965 8 435 0 0 100 1 400 宁夏 31 581 13 823 3 617 47 2 011 内蒙古 41 1 230 22 1 225 5 660 68 3 085 陕西 116 2 457 32 1 890 7 1 250 155 5 597 山西 33 894 80 4 965 25 5 030 138 10 889 河南 208 3 841 29 1 561 5 900 242 6 302 山东 51 969 29 1 535 5 905 85 3 409 总计 1 038 16 030 224 13 274 52 9 632 1 314 38 896 注:数据统计来源,2016—2020年黄河流域九省区煤炭行业化解过剩产能关闭煤矿名单,关停煤矿不计入。 3 废弃矿井抽水蓄能选址评价指标体系
3.1 选址评价指标体系的建立
废弃矿井抽水蓄能电站目前处于探索发展阶段,并非所有废弃矿井都适宜于建设抽水储能电站,在工程建设时首先需要综合考虑各种因素,从废弃矿井中选择适宜建设的矿井。以废弃矿井抽水蓄能电站建设过程中地下水库库容、地下水循环、围岩稳定性与水质因素这4个关键科学问题[8]为基础,考虑废弃矿井区域的经济、社会与资源条件因素,综合大量研究文献[28-37],确定用于选址问题的一级指标(B)6个,二级指标(C)24个,具体如图3所示。
3.2 基于AHP的选址步骤
由于废弃矿井抽水蓄能选址影响因素的模糊性,选用AHP(层次分析法)确定各指标的影响权重,该方法是将复杂问题分解为多个组成因素,并将这些因素按支配关系分组形成递阶层次,进行定性、定量分析评价[38]。具体评估步骤如下:
第一步:确定标度和构造判断矩阵。邀请5位废弃矿井抽水蓄能相关领域专家对选址因素进行重要性评判,采用1—9标度打分法对各影响指标两两比较,构造判断矩阵。
第二步:计算各层次指标所占的权重。计算各判断矩阵的最大特征值及对应的特征向量,利用一致性指标、随机一致性指标和一致性比率进行一致性检验。
一致性指标CI为:
$$ {{\rm{CI}}} = \frac{{{\lambda _{\max }} - n}}{{n - 1}} $$ (1) 式中:λmax为矩阵最大特征根;
$ n $ 为判断矩阵阶数。随机一致性比率CR为:
$$ {{\rm{CR}}} = \frac{{{{\rm{CI}}}}}{{{{\rm{RI}}}}} $$ (2) 式中:RI为平均随机一致性指标。
第三步:一致性检验分析。当CR<0.1,表明通过一致性检验,反之则没有通过一致性检验。当一致性检验未通过时,需返还给专家检查是否存在逻辑问题,重新判断,再次录入判断矩阵进行分析。
3.3 选址因素权重计算
通过5位专家对各项指标的重要性进行评判,建立准则层对目标层(一级指标对目标层)、分析层对准则层(二级指标对一级指标)各自的判断矩阵,采用加权平均方法并进行一致性检验,求得废弃矿井抽水蓄能选址体系中各影响因素权重。所得一致性检验与权重结果见表4—表5。
表 4 一致性检验结果Table 4. Consistency test results专家 一致性检验 判断矩阵 A-B B1-C B2-C B3-C B4-C B5-C B6-C 专家1 λmax 6.015 9 4.061 8 4.068 3 4.166 7 6.211 8 3.009 2 3.094 8 CR 0.002 5 0.023 1 0.023 9 0.062 5 0.033 6 0.008 9 0.091 1 专家2 λmax 6.093 9 4.089 0 4.045 8 4.191 6 6.273 6 3.018 3 3.005 5 CR 0.014 9 0.033 3 0.017 1 0.071 8 0.043 4 0.017 6 0.005 3 专家3 λmax 6.087 2 4.083 9 4.166 2 4.045 4 6.412 7 3.009 2 3.018 3 CR 0.013 8 0.031 4 0.062 2 0.017 0 0.065 5 0.008 9 0.017 6 专家4 λmax 6.059 2 4.051 9 4.109 1 4.159 3 6.577 9 3.044 4 3.053 9 CR 0.009 4 0.019 4 0.040 8 0.059 7 0.091 7 0.042 7 0.051 9 专家5 λmax 6.032 8 4.186 5 4.144 8 4.131 6 6.338 6 3.000 0 3.018 3 CR 0.005 2 0.069 8 0.054 2 0.049 3 0.053 7 0 0.017 6 表 5 废弃矿井抽水蓄能选址各指标权重Table 5. Weights of each index for site selection of pumped storage plant using abandoned mines准则层 权重 指标层 综合权重 排序 空间条件 B1 0.304 4 巷道空间体积C1 0.123 9 3 巷道地形条件C2 0.020 4 12 上下水库水位差C3 0.142 7 1 巷道瓦斯含量C4 0.017 4 15 地质条件 B2 0.299 2 巷道围岩稳定性C5 0.133 2 2 巷道围岩渗透率C6 0.123 3 4 地质灾害类别C7 0.013 9 18 地质灾害频率C8 0.028 7 9 水文条件 B3 0.151 6 水源水质C9 0.014 2 17 矿井涌水量C10 0.030 8 8 水库注水水源C11 0.015 5 16 地下水循环特征C12 0.091 1 5 社会条件 B4 0.073 7 当地用电需求量C13 0.022 2 11 矿区交通便捷度C14 0.005 5 22 与电网距离C15 0.008 3 21 与城市距离C16 0.008 8 20 当地政策支持度C17 0.023 9 10 就业率C18 0.005 0 23 经济条件 B5 0.129 1 单位储能成本C19 0.055 9 7 投资回收年限C20 0.012 8 19 工程收益率C21 0.060 5 6 资源条件 B6 0.042 0 当地有效风速C22 0.019 4 13 平均降水量C23 0.018 7 14 太阳辐射时间C24 0.004 0 24 综合以上计算结果可知:C1(巷道空间体积)、C3(上下水库水位差)、C5(巷道围岩稳定性)、C6(巷道围岩渗透率)、C12(地下水循环特征)是选址的重要指标。其中,上下水库水位差影响着抽水蓄能电站的经济效益和发电效率、装机容量等;巷道空间体积则决定着上下水库容量;巷道围岩稳定性和巷道围岩渗透率影响废弃矿井抽水蓄能电站地下空间重构;地下水循环特征影响水库环保性与机组设备效率。这些指标对于抽水蓄能电站建设的可行性、经济性、安全性与环保性等有着举足轻重的影响,是工程选址过程中需要考虑的关键因素。
4 黄河流域典型废弃矿井抽水蓄能电站设计实例
利用废弃矿井抽水蓄能选址评价方法,依据专家对黄河流域九省区5座废弃矿井综合评价结果,选取得分最高的某矿建立全地下式抽水蓄能电站(图4)。该矿地下空间充足,地表土地并未形成塌陷区,采取全地下式抽水蓄能电站设计。电站上水库选用一采区回风巷与运输大巷联合构建,并选取2条联络巷保证2条巷道之间较好的水体流动性。下水库由三采区东西两翼运输巷道组成,同样采用车场与联络巷作为水体流动性保障。电站采用轨道上下山与地表相连接,采用运输斜巷、二采区运输上山、二采区车场、三采区运输上山等组成输水通道,通过井巷参数确定该抽水蓄能电站主要参数(表6)。假设废弃矿井抽水蓄能电站所产生的电力均并入电网并出售,抽水时间和发电时间为6 h,巷道所建地下水库不考虑死库容,水库运行期间全部排空。
表 6 废弃矿井抽水蓄能电站基本参数Table 6. Basic parameters of pumped storage power plant using abandoned mine参数 数值 上水库平均水位/m 1 210 上水库总容积/m3 34 600 下水库正常储水高度/m 4.5 下水库正常储水位高程/m 830.5 下水库总容积/m3 41 200 上下水库自然高程差/m 338 上下水库容积之比 0.84 输水巷道长度/m 2 016 4.1 电站耗电量
(1) 平均水头(即为上下水库的平均高程差):
$$ \overline H = {H_{{\rm{d}}}} + \left( {{H_{{{\rm{sd}}}}} + \frac{1}{2}{H_{{{\rm{sg}}}}}} \right) - \left( {{H_{{{\rm{ud}}}}} + \frac{1}{2}{H_{{{\rm{ug}}}}}} \right) $$ (3) 式中:Hd为上水库和下水库的自然高程差,338 m;Hsd、Hud分别为上下水库死水位,0 m;Hsg、Hug分别为上下水库储水深度,分别取值92、4.5 m。计算可得平均水头
$ \overline H $ 为381.75 m。(2) 抽水流量:
$$ Q = \frac{V}{t} $$ (4) 式中:V为上水库库容,34 600 m3;t为抽水时间,6 h;计算抽水流量
$ Q $ 为1.60 m3/s。(3) 输水路管道直径:
$$ d = \sqrt {\frac{{4{Q_{\max }}}}{{{{\text{π}} }{{v}}}}} $$ (5) 式中:Qmax为管道最大水流量,取抽水流量1.60 m3/s[39];v为输水管中的水流速度,取经济流速2.4 m/s[40]。计算得出输水路管道直径d约为0.9 m。
(4) 抽水总水头损失量Hp,主要包括沿程水头损失量(Hf)与局部水头损失量(Hj):
$$ {H_{{\rm{f}}}} = \beta \frac{{l{v^2}}}{{2dg}} $$ (6) $$ {H_{{\rm{j}}}} = \xi \frac{{{v^2}}}{{2g}} $$ (7) $$ {H_{{\rm{p}}}} = {H_{{\rm{f}}}} + {H_{{\rm{j}}}} $$ (8) 式中:
$ \beta $ 为输水管糙率系数,取值0.018[41];l为输水巷道长度,2 016 m;g为重力加速度,9.8 m/s2;$\xi $ 为局部水头损失系数,取值1[40]。通过计算得出沿程水头损失量为11.85 m,局部水头损失量为0.29 m,抽水总水头损失量为12.14 m。
(5) 水泵最大耗电量:
$$ {W_{{{\rm{pc}}}}} = \frac{{\left( {\overline H + {H_{{\rm{p}}}}} \right)V{\rho _{{\rm{w}}}}g}}{{{\eta _1}}} $$ (9) 式中:Wpc为水泵最大耗电量,J;ρw为水的密度,1.0×103 kg/m3;η1为电力系统的整体机械效率,取值87%。
水泵的最大耗电量为1.535×1011 J,相当于42 644 kW·h,在抽水时长为6 h的模式下,电机的平均功率为7 107 kW。
4.2 电站发电量
(1) 最大发电量:
$$ {W_{{\rm{pg}}}} = \left( {\overline H + {H_{\rm{t}}}} \right)V{\rho _{\rm{w}}}g{\eta _2} $$ (10) 式中:Wpg为最大发电量;Ht为发电总水头损失,假设其近似与抽水总水头损失相同[40];η2为发电系统的整体机械效率,取值87%。
电机最大发电量为1.162×1011 J,相当于32 277 kW·h,在发电时长为6 h情况下,平均发电机功率为5 380 kW。
(2) 发电效率:
$$ {\eta _0} = \frac{{{W_{{{\rm{pg}}}}}}}{{{W_{{{\rm{pc}}}}}}} \times 100{\text{%}} $$ (11) 式中:η0为废弃矿井抽水蓄能电站发电效率,%。
计算得到该抽水蓄能电站发电效率为75.7%,可媲美常规抽水蓄能电站。
综上可知,该废弃矿井抽水蓄能电站装机容量可以选择为5 000 kW(考虑死库容情况和预留库容情况水库库容会减少),按照每天运行6 h,每月3 d停运检修计算,发电效率约为75.7%,则年发电量可达7.47×106 kW·h。如果建设成本在可控制范围内,则该小型废弃矿井抽水蓄能电站方案具有一定的可行性。
5 废弃矿井抽水蓄能电站运行效益
5.1 静态效益
对于废弃矿井抽水蓄能电站静态效益的评价可以从容量效益、调峰填谷效益、排放效益等[42]方面进行分析。
1) 容量效益
容量效益,是指废弃矿井抽水蓄能电站投入后,减少同规模火电电站建设投入节约的效益。其计算公式[43]如下:
$$ {B_{{\rm{c}}}} = \left[ {{C_0} + {I_0}\left(\frac{A}{P},i,{y_0}\right)} \right] - \left[ {{C_1} + {I_1}\left(\frac{A}{P},i,{y_1}\right)} \right] $$ (12) 式中:Bc为抽水蓄能电站每年的容量效益,万元;C0为火电机组固定运行费用,万元;I0为火电机组的建设费用,万元;A为年终支付金额,万元;P为净值,万元;i为年利率;y0为火电机组成本回收周期,年;C1为抽水蓄能机组运行费用,万元;I1为抽水蓄能机组建设成本,万元;y1为抽水蓄能机组成本回收周期,年。
2) 调峰填谷效益
抽水蓄能重要作用之一就是调峰填谷,其中,调峰效益指抽水蓄能机组运行时减少常规调峰机组燃料所获得的效益;填谷效益是指用电低谷期抽水蓄能电站为保持功率平衡,消耗富余电量进行抽水,避免火电机组降负荷运行,提升火电机组燃料利用率所带来的效益。综合计算公式[42]如下:
$$ {B_{{\rm{p}}}} = \sum\limits_{i = 1}^{365} {({B_0} + B_0^{'} - {B_1} - B_1^{'})} \cdot \omega $$ (13) 式中:Bp为废弃矿井抽水蓄能电站每年调峰填谷效益,万元;
$B_0$ 、$B_0^{'} $ 分别为调峰时火电机组在调峰时段和非调峰时段燃料消耗量;$B_1$ 、$B_1^{'} $ 分别为填谷时火电机组在非填谷时段和填谷时段的燃料消耗量;ω为市场实时煤炭价格。3) 排放效益
排放效益是指废弃矿井抽水蓄能电站发电量替代火电站同等发电量时所减少的有害气体排放量,主要通过减少的二氧化碳、二氧化硫与氮氧化物等排放量进行量化。
4) 建设成本节约效益
建设成本节约效益是指利用废弃矿井进行改建所节约的建设成本,当前常规抽水蓄能电站投资标准为5 000~6 000元/kW。目前关于废弃矿井抽水蓄能电站成本有2种预测结果:第一种利用煤矿建设抽水蓄能电站工程投资为17 959元/ kW[44];第二种基于废弃矿井改造的抽水蓄能电站的投资为2 982元/ kW[39]。这2种建设投资成本差距在于废弃矿井地下空间重构。对于巷道加固密闭所需成本存在争议,因此,目前对于废弃矿井抽水蓄能电站建设成本预估方案,需待废弃矿井抽水蓄能地下空间稳定性与密闭性问题的进一步解决。
5.2 动态效益
抽水蓄能电站具有适应负荷快速变化的特性,机组运行快速灵活,负荷调整范围大,调频、调相性能好,可满足系统运行需要,提高电网安全可靠性,由此产生的效益为动态效益。废弃矿井抽水蓄能动态效益主要体现在调频效益、调相效益、备用效益与黑启动效益。
1) 调频效益
调频效益是抽水蓄能电站代替小型火力机组进行增减荷,以保持频率并提高电网的稳定性,由此带来的减少能源损耗的效益。实际工程中存在一些不确定参数,计算过程中结合经验值及实际情况进行取值,其公式[42]如下:
$$\begin{aligned} &{B_{\rm{f}}} = \sum\limits_{i = 1}^{365} {\left[ {\sum\limits_{i = 1}^n {({{\text{c}}_{1i}} \cdot {x_i} + {{\text{c}}_{2i}} \cdot {T_i} \cdot {D_i}) \cdot {f_i}} } \right]} - \\ &\qquad\sum\limits_{i = 1}^{365} {\left(\frac{{{V_0} \cdot H \cdot \eta }}{{365}} \cdot {\textit{z}_i} \cdot {m_i} \cdot {j_{\rm{s}}}\right)}\end{aligned} $$ (14) 式中:Bf为废弃矿井抽水蓄能电站的调频效益;c1i为第i台调峰火电机组启停时煤炭的消耗量,g/(kW·h);xi为第i台火电机组每日启动次数;c2i为第i台机组调频耗煤量;Ti为第i台机组参与调频的总时长,h;Di为第i台机组的机组容量,万kW;V0为机组启动时消耗水量;H为机组的平均工作水头,m;η为水轮机效率;zi为废弃矿井抽水蓄能电站日启动次数;
$ {m_i} $ 为废弃矿井抽水蓄能电站机组数量;js为当前电价,元/(kW·h)。2) 调相效益
调相效益是废弃矿井抽水蓄能机组替代火电机组,为系统提供调相时节约的火电动态成本[43]。其计算公式为:
$$ \begin{gathered} {B_{{\rm{t}}}} = P{}_0 \cdot U \cdot \cos \theta \cdot \sin \theta \cdot \left(\frac{A}{P},{i_{{\rm{s}}}},{e_{{\rm{s}}}}\right) + \\ \qquad P{}_0 \cdot U \cdot \cos \theta \cdot \sin \theta \cdot 2{\text{%}} \\ \end{gathered} $$ (15) 式中:Bt为废弃矿井抽水蓄能电站年调相效益;P0为不考虑视在容量时机组的单位有功容量投资;U为机组无功容量;cosθ为额定功率因数;is为社会折现率;es为废弃矿井抽水蓄能电站设计使用期限。
3) 备用效益
备用效益是指废弃矿井抽水蓄能机组替代火电机组,为系统提供旋转备用所节约的火电动态成本[43]。其计算公式如下:
$$ {B_{{\rm{s}}}} = E \times 30{\text{%}} \times 2\sigma $$ (16) 式中:Bs为废弃矿井抽水蓄能电站备用效益;E为抽水电站总装机容量;σ为废弃矿井抽水蓄能电站综合效率系数。
4) 黑启动效益
发电厂在失去电源时自行启动,称“黑启动”服务,确保在电力系统局部或整体瓦解时系统重新运行,其对电网产生的经济效益称为黑启动效益[45]。黑启动效益与电网全部停电概率、电网平均负荷、黑启动时间等条件相关,其经济效益属于广义范畴,暂无具体的计算公式。
6 废弃矿井抽水蓄能电站利用效益
由表2可知,我国主要煤炭资源储量大多集中在黄河流域九省,由于煤炭资源的不可再生性,未来废弃矿井数量将不断增加。文献[46]表明2018年我国共有各类废弃矿山约9.9万座,黄河流域九省作为国家矿业生产的主要地区,废弃矿井数量十分巨大。本次以2016—2020年黄河流域各省区废弃煤矿为研究对象,计算其可利用地下空间,采用静态效益和动态效益2个维度,分析其应用潜力。
将表3数据与文献[47]中的全国煤矿规模与其井巷可利用的地下空间量的比例系数相结合,估算得出2016—2020年黄河流域九省区废弃矿井井巷地下空间7 846.47万m3。按照煤矿关停后60%的井巷空间可以有效利用估算[47],黄河流域九省区5年来关闭矿井的可有效利用井巷空间为4.7×107 m3,如图5所示。
因废弃矿井抽水蓄能电站对两水库间距有一定要求,所以图5中关闭煤矿的井巷可利用空间并非都可以作为废弃矿井抽水蓄能电站地下水库载体,需根据矿井开采情况,对废弃煤矿抽水蓄能地下水库可利用巷道进一步估算。根据现有研究表明:当上下水库高程差在200~400 m时,发电机组与设备成熟,但存在因水头高度引起的电站发电效率与经济效益的下降;在400~600 m时,废弃矿井抽水蓄能电站效率与效益均较好;在600~800 m时,电站所需发电机组及设备难度提升,但其发电效率和经济效益会增加;当小于200 m,电站效率较低且建设成本较高;当大于800 m 时,现有的高水头水泵水轮机不能满足废弃矿井抽水蓄能电站需求[48-49]。通过调研收集黄河流域九省区132座关闭煤矿基本情况,其中小型矿井占比70%,中型矿井占比23%,大型矿井占比7%;关闭矿井采深小于200 m的占比19%,200~400 m占比30%,400~600 m占比34%,600~800 m占比17%(图6)。依据黄河流域九省区关闭矿井调研信息,按照矿井各煤层采深初步划分为不适宜建设、半地下式和全地下式3种类型,并按照上下水库高差进一步划分,其中,采深小于200 m且规模小于3万t/a不适宜建设。半地下式与全地下式模式,根据各煤层埋深进行上下水库高差分析。统计结果表明:其中25座不适宜废弃矿井抽水蓄能利用,95座适宜于半地下式抽水蓄能电站建设模式,12座适宜于全地下式抽水蓄能电站建设模式。具体各省市废弃矿井建设模式上下水库高差情况见表7。
表 7 黄河流域九省区部分废弃矿井抽水蓄能电站建设模式Table 7. Construction mode of some pumped storage power plants using abandoned mines in nine provinces of Yellow River Basin省/自治区 半地下式建设矿井数 全地下式建设矿井数 200~400 m 400~600 m 600~800 m 200~400 m 400~600 m 青海 1 0 0 0 1 四川 9 22 4 4 1 甘肃 3 3 2 0 1 宁夏 0 2 0 1 0 内蒙古 3 0 0 0 0 陕西 5 4 1 0 1 山西 5 4 0 0 0 河南 10 7 2 1 1 山东 3 3 2 0 1 本文利用统计学估算的方法,提出各省区废弃煤矿抽水蓄能电站各模式井巷资源量估算公式:
$$ {v_{ik}} = {V_i} \cdot {a_{{{ik}}}} \cdot \tau $$ (17) 式中:vik为i省区废弃煤矿抽水蓄能电站建设k模式可利用井巷地下空间体积;Vi为i省区井巷可利用地下空间总体积;aik为i省区废弃矿井抽水蓄能建设k模式所占比;τ为废弃矿井抽水蓄能电站地下水库的井巷空间比例,取0.27[9]。
各省区废弃煤矿抽水蓄能电站最大可发电量估算公式如下:
$$ {W_i} = \sum\limits_{k=1}^3 {{{\bar H}_{{\rm{s}}k}}{v_{{\rm{s}}k}}{\rho _{\rm{w}}}g{\eta _{{\rm{s}}k}}} + \sum\limits_{k=1}^2 {\frac{{{{\bar H}_{{\rm{a}}k}}{v_{{\rm{a}}k}}{\rho _{\rm{w}}}g{\eta _{{\rm{a}}k}}}}{2}} $$ (18) 式中:Wi为i省区废弃煤矿抽水蓄能电站发电总量;
${\bar H_{{\rm{s}}k}}$ 、${\bar H_{{\rm{a}}k}}$ 分别为i省区半地下式、全地下式抽水蓄能电站建设k模式平均水头高;vsk、vak为i省区半地下式、全地下式抽水蓄能电站建设k模式可利用井巷地下空间总体积;ηsk、ηnk为i省区半地下式、全地下式抽水蓄能电站建设k模式发电效率。通过井巷资源量估算公式计算出,黄河流域九省区5年废弃矿井抽水蓄能电站各建设模式可利用井巷空间为1.34×107 m3(表8)。假设上下水库高差200~400 m的废弃矿井抽水蓄能电站平均净水头高度为300 m,抽水蓄能效率为70%;400~600 m的平均净水头高度为500 m,抽水蓄能效率为75%;600~800 m的平均净水头高度为700 m,抽水蓄能效率为80%。按照抽水蓄能电站1年运行329 d,通过式(18)计算得出2016—2020年黄河流域九省区废弃煤矿抽水蓄能电站每年发电总量约为3.78×109 kW·h(图7),根据抽水蓄能电站电能消耗与发电量之比约4∶3,则每年可消纳电能5.04×109 kW·h。
表 8 2016—2020黄河流域九省区废弃矿井抽水蓄能电站可利用井巷资源量Table 8. Available shaft and drift resources for pumped storage power plants in abandoned mines in nine provinces of Yellow River Basin in 2016-2020单位:万m3 省/自治区 半地下式 全地下式 [200, 400) m [400, 600) m [600, 800) m [200, 400) m [400, 600) m 青海 1.74 0.00 0.00 0 1.7424 四川 72.12 176.46 32.08 32.08 8.02 甘肃 18.12 18.12 12.08 0 6.04 宁夏 0 20.75 0 10.36 0 内蒙古 27.44 0 0 0 0 陕西 72.60 58.08 14.52 0 14.52 山西 103.43 82.74 0 0 0 河南 181.72 127.20 36.34 18.17 18.17 山东 58.41 58.41 38.94 0 19.47 按照2021年全国平均标准煤耗为302.5 g/(kW·h)计算,2016—2020黄河流域九省区废弃矿井抽水蓄能电站全年发电量相当于节约了114.3万t标准煤。2021年,黄河流域九省区累计弃风电量为123.5亿kW·h,弃光电量为54亿kW·h,总计177.5亿kW·h,如果利用这些废弃煤矿建设抽水蓄能电站可满足流域内28.4%的弃风、弃光电量的消纳需求。统计2021年黄河流域九省区峰谷电价差数据(表9),选用各地区峰谷电价差平均值与废弃矿井抽水蓄能发电量计算,求得2016—2020黄河流域废弃煤矿抽水蓄能电站的直接经济效益约达21.2亿元/年。
表 9 2021年黄河流域九省区峰谷电价差统计数据Table 9. Statistical data of peak-to-valley electricity price difference in nine provinces of Yellow River Basin in 2021地区 一般工商业峰谷价差/(元·kW−1·h−1) 大工业峰谷价差/(元·kW−1·h−1) <1 kV 1~10(20) kV ≥35 kV 1~10(20) kV 35 kV 110 kV 220 kV 青海 0.444 6 0.438 4 0.398 8 0.424 2 0.411 6 0.399 0 0.390 4 四川 0.615 2 0.602 2 0.589 2 0.540 2 0.520 2 0.498 2 0.475 2 甘肃 0.581 4 0.571 4 0.561 4 0.438 3 0.428 3 0.418 3 0.409 3 宁夏 0.407 6 0.391 6 0.375 6 0.318 8 0.294 8 0.278 8 0.262 8 内蒙古东部 0.717 1 0.680 0 0.594 3 0.492 1 0.485 1 0.445 7 0.422 7 内蒙古西部 0.217 9 0.217 9 0.217 9 0.217 9 0.217 9 0.217 9 0.217 9 陕西榆林 0.564 7 0.544 7 0.524 7 0.530 7 0.511 8 0.486 6 0.486 6 陕西其他 0.577 4 0.557 4 0.537 4 0.633 4 0.608 2 0.583 0 0.576 8 山西 0.455 1 0.435 8 0.421 3 0.438 7 0.409 8 0.385 6 0.366 3 河南 0.624 5 0.596 0 0.566 4 0.622 4 0.606 3 0.590 3 0.581 7 山东 0.744 5 0.728 1 0.711 6 0.735 7 0.717 7 0.699 7 0.681 7 7 结 论
a. 提出黄河流域九省区废弃矿井抽水蓄能电站建设的半地下式、全地下式2种模式,并建立抽水蓄能电站选址决策指标体系,包含6个一级指标、24个二级指标,其中,巷道空间体积、上下水库水位差、巷道围岩稳定性、巷道围岩渗透率、地下水循环特征是影响选址的重要因素。
b. 开展了黄河流域某矿改建抽水蓄能电站的实例分析,计算出矿山可选择机组装机容量为5 000 kW,年发电量可达7.47×106 kW·h,发电效率可达到75.7%,对废弃矿井抽水蓄能电站的可行性具有一定的参考价值与指导意义。
c. 统计出黄河流域九省区2016—2020年能源消耗与关闭煤矿情况,据此计算出有效利用井巷空间约4.7×107 m3,通过统计学分析计算出黄河流域九省区5年来废弃矿井抽水蓄能电站可利用井巷空间为1.34×107 m3,进而预估利用其建设抽水蓄能电站每年可发电3.78×109 kW·h,可满足黄河流域九省区28.4%的弃风、弃光电量的消纳需求,每年直接经济效益约达21.2亿元。
-
[1] LIU J, WANG J G, GAO F, et al. Impact of micro-and macro-scale consistent flows on well performance in fractured shale gas reservoirs[J]. Journal of Natural Gas Science and Engineering, 2016, 36(Part B):1239-1252.
[2] 杨峰, 宁正福, 胡昌蓬, 等. 页岩储层微观孔隙结构特征[J]. 石油学报, 2013, 34(2):301-311. YANG Feng, NING Zhengfu, HU Changpeng, et al. Charaterization of microscopic pore structres in shale reservoirs[J]. Acta Petrolei Sinica, 2013, 34(2):301-311.
[3] 谢和平,高峰,鞠杨,等. 页岩气储层改造的体破裂理论与技术构想[J]. 科学通报, 2016, 61(1):36-46. XIE Heping, GAO Feng, JU Yang, et al. Novel idea of the theory and application of 3D volume fracturing for stimulation of shale gas reservoirs(in Chinese)[J]. Chinese Science Bulletin, 2016, 61(1):36-46.
[4] 庄茁,柳占立,王永亮. 页岩油气高效开发中的基础理论与关键力学问题[J]. 力学季刊, 2015,36(1):11-25. ZHUANG Zhuo, LIU Zhanli, WANG Yongliang. Fundamental theory and key mechanics problems of shale oil gas effective extraction[J]. Chinese Quarterly of Mechanics, 2015,36(1):11-25.
[5] 梁冰,代媛媛,陈天宇,等. 裂缝参数对考虑滑脱效应页岩气水平井产能的影响[J]. 煤田地质与勘探, 2015,43(2):44-47. LIANG Bing, DAI Yuanyuan, CHEN Tianyu, et al. Influence of fissure parameters on production capacity of shale gas horizontal wells considering slippage effect[J]. Coal Geology & Exploration, 2015,43(2):44-47.
[6] 邹才能,朱如凯,白斌,等. 中国油气储层中纳米孔首次发现及其科学价值[J]. 岩石学报, 2011,27(6):1857-1864. ZOU Caineng, ZHU Rukai, BAI Bin, et al. First discovery of nano-pore throat in oil and gas reservoir in China and its scientific value[J]. Acta Petrologica Sinica, 2011,27(6):1857-1864.
[7] LIU J, WANG J G, GAO F, et al. Flow consistency between non-darcy flow in fracture network and nonlinear diffusion in matrix to gas production rate in fractured shale gas reservoirs[J]. Transport in Porous Media, 2016, 111(1):97-121.
[8] 糜利栋,姜汉桥,李涛,等. 基于离散裂缝模型的页岩气动态特征分析[J]. 中国石油大学学报(自然科学版), 2015,39(3):126-131. MI Lidong, JIANG Hanqiao, LI Tao, et al. Characterization and dynamic analysis of shale gas production based on discrete fracture model[J]. Journal of China University of Petroleum(Edition of Natural Science), 2015,39(3):126-131.
[9] PENG Yan, LIU Jishan, PAN Zhejun, et al. A sequential model of shale gas transport under the influence of fully coupled multiple processes[J]. Journal of Natural Gas Science and Engineering, 2015, 27:808-821.
[10] CAO Peng, LIU Jishan, LEONG Y-K. A fully coupled multiscale shale deformation-gas transport model for the evaluation of shale gas extraction[J]. Fuel, 2016, 178:103-117.
[11] MI Lidong, JIANG Hanqiao, LI Junjian, et al. The investigation of fracture aperture effect on shale gas transport using discrete fracture model[J]. Journal of Natural Gas Science and Engineering, 2014, 21:631-635.
[12] 陈必光,宋二祥,程晓辉. 二维裂隙岩体渗流传热的离散裂隙网络模型数值计算方法[J]. 岩石力学与工程学报, 2014, 33(1):43-51. CHEN Biguang, SONG Erxiang, CHENG Xiaohui. A numerical method for discrete fracture network model for flow and heat transfer in two-dimensional fractured rocks[J]. Chinese Journal of Rock Mechanics and Engineering, 2014,33(1):43-51.
[13] ODA M. An equivalent continuum model for coupled stress and fluid flow analysis in jointed rock masses[J]. Water Resources Research, 1986, 22(13):1845-1856.
[14] GARCIA-TEIJEIRO X, RODRIGUEZ-HERRERA A, FISCHER K. The interplay between natural fractures and stress as controls to hydraulic fracture geometry in depleted reservoirs[J]. Journal of Natural Gas Science and Engineering, 2016, 34:318-330.
-
期刊类型引用(17)
1. 刘德旭,杨迎,黄宏旭,钟海旺. 新型电力系统大规模抽水蓄能调度运行与控制综述及展望. 中国电机工程学报. 2025(01): 80-98 . 百度学术
2. 王双明,宋世杰,魏江波,刘浪,牛瑞琳,李源红,谈庆钰,田雨. 陕北黄土覆盖区煤层上覆岩系层状结构对开采区土壤损伤效应研究. 煤炭学报. 2025(01): 23-42 . 百度学术
3. 李大成,朱天生,吴迪,马黎,项华伟. 基于滑动叠加组合的抽水蓄能电站洪水调节计算模型及应用. 中国农村水利水电. 2024(01): 244-250+256 . 百度学术
4. 叶润成,李茂斌,梁毅. 基于EOD模式的“矿山抽蓄+生态修复”开发前景思考. 能源与环保. 2024(03): 40-46 . 百度学术
5. 修雅馨,刘钦节,付强,杨卿干,包兴芮,邓卓越. 废弃矿井地下空间物理储能方式对比与优选. 绿色矿冶. 2024(02): 6-13 . 百度学术
6. 文志杰,姜鹏飞,宋振骐,李利平,蒋宇静,梁桐凯. 关闭/废弃矿井抽水蓄能开发利用现状与进展. 煤炭学报. 2024(03): 1358-1374 . 百度学术
7. 王方田,张辰凯,张村,贾胜. “双碳”目标驱动下煤系伴生资源协同开发技术框架. 科技导报. 2024(07): 40-55 . 百度学术
8. 孙中博,赵毅鑫,任建东. 废弃矿井地下空间再利用根源因素识别及评价模型构建. 中国矿业大学学报. 2024(03): 585-599+612 . 百度学术
9. 刘飞跃,陈澫赟,杨科,段敏克,刘钦节,杨凌越. 安徽省关闭煤矿资源调查及其综合利用研究. 煤田地质与勘探. 2024(06): 1-11 . 本站查看
10. 宋世杰,张家杰,杨帅,刘志坚,江宁,陈平,刘汉斌,魏江波,刘露,陈宝灯,李源红. 黄河上中游采煤沉陷区水土流失效应的探索与思考. 绿色矿山. 2024(02): 169-182 . 百度学术
11. 霍超,王蕾,谢志清,潘海洋,方惠明,赵岳,王丹丹,王丹凤. 新时期我国煤矿地下空间综合利用现状及展望. 地质论评. 2024(04): 1455-1468 . 百度学术
12. 文志杰,姜鹏飞,宋振骐,蒋宇静. 废弃矿井抽水蓄能面临的关键问题与对策思考. 山东科技大学学报(自然科学版). 2023(01): 28-37 . 百度学术
13. 顾清华,李学现,卢才武,阮顺领,江松. “双碳”背景下露天矿智能化建设新模式的技术路径. 金属矿山. 2023(05): 1-13 . 百度学术
14. 王双明,申艳军,宋世杰,刘浪,顾霖骏,魏江波. “双碳”目标下煤炭能源地位变化与绿色低碳开发. 煤炭学报. 2023(07): 2599-2612 . 百度学术
15. 于凯,黄馨阅. 基于知识图谱的抽水蓄能专利信息挖掘与可视化研究. 科学技术与工程. 2023(26): 11364-11374 . 百度学术
16. 高学平,袁泽雨. 抽水蓄能电站地下水库建设进展及关键水力学问题. 水利学报. 2023(09): 1058-1069 . 百度学术
17. 刘具,秦坤. 我国煤炭绿色开采技术进展. 矿业安全与环保. 2023(06): 7-15 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 63
- HTML全文浏览量: 12
- PDF下载量: 5
- 被引次数: 19