Study on Water Level Decline Control of Yueyaquan Lake in Dunhuang
-
摘要:
敦煌鸣沙山月牙泉以“山泉共处，沙水共生”的奇妙景观著称于世，被誉为“塞外风光之一绝”，是甘肃省著名风景名胜区。但从20世纪70年代中期开始，由于党河水库的修建，河道断流，地下水补给量减少，区域地下水位下降，引起月牙泉湖水位急剧下降，泉湖底多次部分露出水面，造成周围地质环境恶化，严重影响敦煌旅游业的发展，引起社会广泛关注。自1986年开始，有关部门多次对月牙泉湖实施抢救措施。笔者阐述了月牙泉湖水文地质条件、水位下降的原因及治理过程，运用FEFLOW软件模拟分析预测月牙泉湖不同治理方案情况下，月牙泉湖水位及面积变化情况。2018年4月，月牙泉湖治理工程实施后，月牙泉湖水面逐渐上升扩大，效果良好。月牙泉湖水位下降治理的成功经验对于西北干旱地区环境治理与修复具有典型的示范和推广价值。
Abstract:Mingsha mountain and Yueyaquan lake in the city of Dunhuang are famous for their wonderful landscapes of “coexistence of mountains and springs, symbiosis of sand and water”. They are known as “one of the most beautiful scenery outside the Great Wall” and are famous scenic spots in Gansu province. However, since the mid–1970s, due to the construction of the Danghe reservoir, the river was cut off, the groundwater replenishment was reduced, and the regional groundwater level dropped, causing the water level of the Yueyaquan lake to drop sharply and the bottom of the spring lake had repeatedly partially exposed to the surface of the water, causing the deterioration of the surrounding geological environment, seriously affecting the development of Dunhuang’s tourism industry, and causing widespread concern in the society. Since 1986, relevant departments have implemented multiple rescue measures on Yueyaquan lake. This paper describes the hydrogeological conditions of Yueyaquan lake, the causes of water level decline and the treatment process, FEFLOW software is used to simulate, analyze and predict the different treatment scenarios of the Yueyaquan Lake, changes in water level and area of the Lake. After the implementation of the Yueyaquan lake treatment project in April 2018, the water surface of the Yueyaquan lake has gradually risen and expanded, with good results. The successful experience in the treatment of the decline in the water level of the Yueyaquan lake has a typical demonstration and promotion value for environmental treatment and restoration in the arid area of Northwest China.
-
Keywords:
- Yueyaquan lake in Dunhuang /
- water level decline /
- governance research
-
蒸汽辅助重力泄油（SAGD）是重油油藏开发的一项关键技术。该技术通过在水平生产井上方部署注入井持续注入高温蒸汽，形成高温蒸汽腔，在改善地下油藏的流动性的同时其前缘以热传导的形式增加地下岩石局部的孔隙压力，导致岩石发生微破裂，其破裂能量以弹性波的形式进行传播，这种现象也被称为“热凯赛尔效应”（ 李彦平等，2003；董宏等，2017）。SAGD快速预热技术利用短时间、人工精细调节的高压蒸汽注入，在SAGD的生产井和注汽井之间形成一个高孔隙度和高渗透率的扩容区域，从而提早建立生产井和注入井的均匀连通。这项技术能大大的缩短SAGD生产的预热周期，解决注汽量大、阶段能耗大、地面循环返出液处理压力大、建产慢等多种问题。
微地震监测技术可以连续监测热蒸汽注入过程中诱发地震事件在时间和空间域的动态变化，通过反演的震源空间与震级能量的分布，可描述高温蒸汽腔体的几何形态、均匀程度和形成过程，预测蒸汽腔前缘扩展速度与扩展规律，从而为重油油田的注采参数调整、开采方案的决策和优化提供理论依据（宋维琪，2014；杨瑞召，2017；李萌，2020）。震源定位是微地震监测的核心技术，其结果直接影响蒸汽腔体几何形态和扩展规律的预测精度。目前，常规的微地震震源定位方法需要拾取准确的初至，通过射线追踪最小化拾取与计算初至之差来获取震源的位置（Wang et al.，2016；Li et al.，2019a）。然而，实际采集的SAGD微地震数据普遍包含强噪声干扰，难以拾取准确可靠的初至，降低了基于初至信息定位方法的精度和可靠性。
微地震逆时干涉震源定位技术（TRI）基于逆时偏移原理，同时利用微地震波场的运动学和动力学特征，通过有限差分或相移法进行波场的逆时延拓并施加一定的成像条件来获取震源的位置和起震时刻。该方法无需拾取震相走时，避免了人工拾取初至引入的定位误差，是一种适用于低信噪比微地震监测数据的震源定位方法（Gajewski，2005；Artman et al.，2010；Maxwell，2015；李萌等，2016）。近年来，优化改进逆时干涉定位成像条件以获取高精度震源定位结果成为了研究热点。Zhu（2014）在粘弹性介质中对逆时干涉定位法进行了介质固有衰减的校正，但忽略了散射衰减对定位的影响。Douma等（2015）利用地震勘探反褶积思想提高了逆时干涉定位的空间和时间分辨率。Nakata等（2016）对各接收器波场进行独立逆时传播并将这些子波场的乘积作为最终的干涉波场从而提高定位结果的聚焦性。Li等（2019b）采用多尺度双树复小波变换将原始数据分解到多个时频域分别进行逆时干涉，选择有效尺度TRI的乘积压制了稀疏观测系统引起的成像假象。Li等（2021）引入二维VTI介质波场解耦算法，采用解耦qP和SV波场的自相关和互相关压制了耦合波场互相干涉引起的定位成像假象，提高了逆时干涉震源定位的精度。大量研究表明，逆时干涉震源定位在速度模型合理的情况下，对低信噪比微地震监测数据具有较高的定位精度。
笔者优选风城油田八道湾组超稠油油藏的一组SAGD快速预热注入与生产井组进行为期5天的微地震井中监测，采用逆时干涉定位算法成功实现弱诱发地震事件的震源定位，以此估算不同注入时期的蒸汽腔体在空间的几何展布形态，得到其在水平方向和深度方向的扩展规律，分析岩性夹层和非均质性对蒸汽腔体延伸形态的影响，为该地区SAGD重油油田开发方案的布置与优化奠定基础。
1. 观测系统布置方案
研究区块位于准噶尔盆地西部隆起乌夏断裂带、夏红北断裂上盘的风城油田重58井区，其整体构造形态为被断裂切割的南倾单斜，南部发育断裂，走向近东西向，倾向向北。区块目的层为受构造和岩性控制的单斜超稠油油藏，岩性以砾砂岩为主，油质稠，原油粘度较高，储层物性较差（何文军等，2023）。研究区块优选I井作为高温蒸汽注入井，P井为生产井，F井为微地震监测井，布置8级井下三分量检波器，进行为期5天的井中监测，实时记录高温蒸汽注入时的诱发地震信号。项目选取邻井S井布置射孔信号，校正井下检波器水平分量的方位（图1）。图1a中红色轨迹为选取的井中微地震监测的SAGD快速预热井组。图1b中红色三角为三分量检波器，蓝色星形为射孔点。观测系统以监测井井口坐标为坐标原点，检波器阵列中点距离SAGD井组小于500 m，满足井中微地震监测条件。
![]()
图 1 研究区域井位部署图（a）和井中微地震观测系统（b）Figure 1. (a) Well deployment of studied area and (b) downhole microseismic recording geometry2. 井中微地震逆时干涉震源定位方法
传统的微地震震源定位方法需要拾取准确的纵波和横波初至，采用网格搜索法通过三维射线追踪最小化拾取初至与理论初至之差来推断震源的位置。这类方法的定位精度强烈依赖于拾取初至的准确性，在SAGD快速预热中产生的低信噪比微地震监测中难以取得满意的震源定位结果。逆时干涉震源定位方法基于三维弹性波波动方程，能够同时利用微地震波场的运动学和动力学特征。这类方法通过在检波器位置逆时加载实际记录的原始波形，当速度模型合理时，各检波器逆时传播的微地震波场能量会在真实震源位置干涉叠加到达最大值，该成像最大能量所对应的网格即为震源位置。这类方法无需拾取纵波与横波初至，避免了人工拾取初至引入的定位误差，适用于SAGD快速预热中产生的低信噪比微地震事件的定位。
微地震逆时干涉震源定位算法的基本假设为波场互易原理及弹性波动方程在非耗散介质中的时不变性。该方法以检波器作为震源，以实际接收的微地震原始数据作为震源子波，采用时间域或频率域有限差分或有限元方法进行波场的逆时重构，通过施加合适的成像条件获取逆时干涉震源定位成像，拾取定位成像在空间分布的极大值即可估算震源的位置。以二维弹性波方程为例，由下式估算每个时间步长逆时传播的波场：
$$ \rho \frac{{{\partial ^2}{\mathbf{u}}}}{{\partial {t^2}}} - (\lambda + 2\mu )\nabla \nabla \cdot {\mathbf{u}} + \mu \nabla \times \nabla \times {\mathbf{u}} = {\mathbf{d}}(x,z,T - t) $$ (1) 式中：d为接收到的微地震波形数据；t为时间步长；T为逆时延拓的总时长；
$ \lambda 和\mu $ 为拉梅系数；$ \rho $ 为介质密度；u为重构的逆时波场。为了保证逆时传播的波场能量最终能够收敛，T可设置为微地震主事件波形持续时间的2倍。目前，常用的成像条件是通过选定时间窗口逆时延拓波场模的积分作为成像条件，从而消除时间轴，得到只与空间位置相关的逆时干涉定位成像：$$ I(x,z) = \int\limits_{{T_1}}^{{T_2}} {|{\mathbf{u}}|dt} $$ (2) 式中：|u|为逆时波场的模；T1和T2为起始时刻和终止时刻；I为成像函数，逆时干涉定位成像的最大值即为推断的震源位置（x0，z0）。震源的起震时刻为逆时延拓波场在震源位置最大能量所对应的时刻如公式（3）所示。
$$ {t_0} = \max [{\mathbf{u}}({x_0},{z_0},t)] $$ (3) 微地震事件的震级采用公式（4）进行估算矩震级Mw。
$$ {M_w} = \frac{2}{3}{\log _{10}}({M_0}) - 6.06 $$ (4) 式中：M0为地震矩，可由地震波形的振幅谱进行估算，如公式（5）所示。
$$ {M_0} = \frac{{4\pi \rho {v^3}r{W_0}}}{U} $$ (5) 式中：v为波形速度；r为震源到接收器的距离；ρ为密度；U为震源辐射特性；W0为振幅谱的低频截止值。对于以剪切应力为主的震源可近似取Up = 0.52和Us = 0.63。由于单垂直井的微地震监测系统缺乏水平方向上对震源位置的约束，其逆时干涉震源定位成像在x–y平面为以检波器为中心、以震源至检波器距离为半径的圆，该圆上任何一点都是潜在的震源位置。因此，为了避免方位不确定性，采用微地震事件P波的偏振方向作为约束，选择距离热蒸汽注入中心较近的震源为最终确定的震源位置。图2为井中微地震监测不同时刻逆时延拓的波场快照，当波场传播至起震时刻时，逆时传播的波场会在真实的震源位置处干涉叠加至最大值。由于观测系统仅能采集到来自震源一侧的信息，因此震源逆时干涉定位成像为一个椭圆区域而非一个完美的点，该区域的大小为评价定位不确定性提供了理论依据。
![]()
图 2 井中微地震逆时延拓的波场在不同时刻的波场图a. t=25 m；b. t=15 ms；c. t=0 msFigure 2. Snapshots of time reversed wavefiled in various time for borehole microseismic3. SAGD快速预热微地震逆时干涉震源定位结果分析
3.1 SAGD快速预热微地震事件空间分布分析
实际采集的SAGD快速预热期间产生的微地震典型波形显示（图3），蓝、红和黑实线分别为同一道检波器的X、Y和Z分量波形，黑虚线与蓝虚线分别为拾取的纵波和横波初至。由于SAGD快速预热期间流体注入的压力较小，因此产生的微地震有效事件波形幅度较弱且包含强噪声干扰。尤其是部分纵波初至被湮没在背景噪声之中，难以准确地拾取，使得基于初至信息的常规定位方法可靠性较低。因此，笔者采用三维弹性波逆时干涉震源定位算法，对5天的SAGD高温蒸汽快速预热井中微地震监测数据进行处理，成功定位了33个有效的微地震事件。逆时干涉震源定位采用监测井F井的声波测井曲线，对监测目的层段进行速度分层，建立一维P波层状速度模型（图4a）。由于缺失横波测井信息，S波速度模型由P波速度模型与研究区域的纵横波比经验值进行估算。
![]()
图 3 SAGD快速预热期间记录的井中微地震波形图Figure 3. Measured microseismic waveforms of downhole array during fast SAGD start–up enhancement![]()
图 4 速度模型（a）、逆时干涉震源定位成像俯视图（b）与侧视图（c）Figure 4. (a) Velocity model, (b) map view and (c) side view of time reversed image for event location本次研究的检波器安置在近似垂直的监测井中，缺乏在水平面上对震源位置的约束。因此，笔者采用微地震有效事件P波的偏振方向作为震源的方向，以该方向高斯射线束与逆时干涉定位成像x–y切片的乘积作为最终的成像结果（图4b）。由于单垂直井监测观测系统的震源定位结果具有180°的不确定性，来自真实震源位置和其关于检波器中心对称的位置都有可能是潜在的震源。为解决这个问题，以SAGD注入井的注入中心（图4b中的红星）为基准，近似认为与注入中心距离较近的成像椭圆的中心为真实的震源位置（图4b中的红圈）。图4c该事件三维弹性波逆时干涉震源定位结果在震源y坐标的x–z切片结果（图4c）。受制于稀疏的检波器阵列，逆时干涉成像在检波器阵列存在成像假象。由于SAGD井组的整体深度与检波器阵列存在一定的深度偏差，因此该成像假象不足以干扰真实震源位置的识别。
拾取所有事件逆时干涉震源定位成像最大能量所对应的位置作为估算的震源空间坐标，得到SAGD快速预热期间的微地震震源空间分布（图5）。由此可粗略估算此次注入预热产生的裂缝缝长为570.19 m，缝高为10.56 m，裂缝波及宽度为127.85 m，裂缝总体延伸方位为NE 95.21°，扩容蒸汽腔体的体积为7.7×105 m3。根据震源位置的空间分布得到以下几点认识。
![]()
图 5 SAGD井组微地震震源俯视图（a）与侧视图（b）Figure 5. (a) Map view and (b) side view of microseismic events for SAGD well（1）I井水平段垂深之上的事件有10个，占总数的30.3%，即 I井上方的非均质物性夹层（不均匀泥岩夹层）发育较多微地震事件，说明I井上方储层产生了较多的微裂缝，已经达到了储层预热改造的目的。
（2）I井与P井水平段垂深之间的事件有18个，占总数的54.6%，说明I井与P井之间已经连通。
（3）J1b油层顶界面之上的事件有0个，说明扩容区域有部分未延伸至油层顶界面。
（4）J1b油层底界面之下的事件有2个，说明扩容区域有小部分延伸至油层底界面。
（5）以I井水平段在水平面的投影为界，I井水平段以北有19个微地震事件，I井水平段以南有14个微地震事件，说明预热期间水平段以北更容易发生破裂、预热效果较好，而水平段以南由于存在非均匀夹层，使得微地震事件个数降低。
3.2 SAGD快速预热微地震空间分布动态分析
SAGD快速预热高温蒸汽注入每天所产生诱发微地震事件分布显示（图6，表1），扩容之前（9月25日）由于排量和压力较小，未能观测到有效的微地震事件；第一阶段扩容（9月25日～ 9月27日）期间，压力和排量增大，地层孔隙压力局部集中导致地层开始破裂，其对应的微地震事件个数急剧增加，这段时间微地震事件个数为23，占总数的69.7%；第二阶段扩容（9月28日04点～9月28日18点）对应的微地震事件个数为10个，占总数的30.3%；9月28日18点停泵，再未观测到微地震事件。
![]()
图 6 高温蒸汽注入诱发的微地震事件个数日期分布图Figure 6. The variation of microseismic event number with date induced by high temperature steam injection表 1 SAGD蒸汽腔几何参数随时间变化表Table 1. Geometric parameters of SAGD steam chamber with various date日期 蒸汽腔网络长（m） 蒸汽腔网络宽（m） 蒸汽腔网络高（m） 蒸汽腔
网络方位累积事
件个数蒸汽腔体积（m3） 9月25日 215.14 51.64 1.04 NE113.78° 4 0.11×105 9月26日 296.26 100.91 2.31 NE98.10° 11 0.69×105 9月27日 570.32 107.04 5.24 NE95.99° 23 3.20×105 9月28日 570.19 127.85 10.56 NE95.21° 33 7.70×105 为进一步分析高温蒸汽注入期间形成蒸汽腔体的发育情况，分别展示9月25日至9月28日的累积微地震事件空间分布（图7、图8）。图中不同颜色的圆点代表不同日期新增的微地震事件。9月25日开始压力和排量增大，地层开始发生破裂，此时的微地震事件主要集中在I井的水平段中心附近；9月26日，扩容段进一步增加，扩容区域沿着I井水平段进行延伸，此时微地震事件已经向上延伸，但还未延伸至油层顶界面；9月27日，扩容区域在水平方向上已经趋于稳定，在垂直方向上进一步延伸，且I井和P井已经连通；9月28日，随着第二段扩容，P井之下开始出现微地震事件，扩容区域向下延伸至油层底界面。
![]()
图 7 高温蒸汽腔空间分布动态扩展过程俯视图a. 9 月 25 日的微地震事件空间分布；b. 9 月 26 日的微地震事件空间分布；c. 9 月 27 日的微地震事件空间分布；d. 9 月 28 日的微地震事件空间分布Figure 7. The geometric extension of high temperature steam chamber (map view)![]()
图 8 高温蒸汽腔空间分布动态扩展过程侧视图a. 9 月 25 日的微地震事件空间分布；b. 9 月 26 日的微地震事件空间分布；c. 9 月 27 日的微地震事件空间分布；d. 9 月 28 日的微地震事件空间分布Figure 8. The geometric extension of high temperature steam chamber (side view)3.3 SAGD快速预热微地震能量分布动态分析
采用公式（4）和公式（5）计算已定位微地震有效事件的矩震级，对不同日期的震级能量进行三维插值和平滑，得到注入期间蒸汽腔体震源能量的动态变化（图9、图10）。其中红色区域为震级较强的部分，蓝色为震级较弱的部分。由图可知，以I井水平段为界，微地震能量分布在I井南北方向分布并不对称，其中北部事件较多、能量较为集中，这与北部微地震事件多而南部微地震事件少有关。整体上看，扩容区域按照震级大小大致可分为2个部分，如图9和图10d的虚线所示。两根虚线中间的区域震级较小，而两边的区域震级能量较大，这与本次监测的观测系统有关，在距离检波器阵列较远的位置只有震级能量较大的事件才能被接收，而震级较小的微地震事件在高温蒸汽腔中快速衰减，难以被监测。
![]()
图 9 高温蒸汽腔能量分布动态扩展过程俯视图a. 9 月 25 日的微地震事件能量分布；b. 9 月 26 日的微地震事件能量分布；c. 9 月 27 日的微地震事件能量分布；d. 9 月 28 日的微地震事件能量分布Figure 9. The energy extension of high temperature steam chamber (map view)![]()
图 10 高温蒸汽腔能量分布动态扩展过程侧视图a. 9 月 25 日的微地震事件能量分布；b. 9 月 26 日的微地震事件能量分布；c. 9 月 27 日的微地震事件能量分布；d. 9 月 28 日的微地震事件能量分布Figure 10. The energy extension of high temperature steam chamber (side view)4. 结论
（1）井中三维弹性波逆时干涉震源定位算法适用于SAGD快速预热产生的弱微地震有效事件震源定位，能够提供高精度的震源空间分布形态。高温蒸汽注入形成的蒸汽腔体长度为570.19 m，高度为10.56 m，宽度为127.85 m，长宽比为4.5，总体延伸方位为NE 95.21°，由此可粗略估算此次注入预热的扩容区域为7.7×105 m3。
（2）监测结果表明，I井上方储层产生了较多的微裂缝，已击破I井上方的非均质物性夹层，达到I井上方储层改造的目的；I井与P井之间已经连通，且扩容区域已部分延伸至油层底界面，但是未延伸至油层顶界面。
（3）SAGD快速预热期间产生的微地震能量在I井南北方向分布并不对称，以I井水平段为界，北部的微地震能量较为集中且事件个数更多，表明蒸汽腔体在I井水平段北部延伸更广。
致谢：本次研究由陕西省自然科学基础研究计划面上资助项目（2023-JC-YB-220）资助，中石油新疆油田分公司提供了宝贵数据资料。
-
![]()
图 1 敦煌月牙泉地区地下水埋深及等水位线图
Figure 1. Groundwater depth and water level map in the Yueyaquan area of Dunhuang
![]()
图 2 党河河道地下水与月牙泉水位关系剖面图
Figure 2. Section diagram of relationship between water level of the Danghe river and groundwater of the Yueyaquan
![]()
图 3 1947～2020年月牙泉湖水位动态变化趋势图
Figure 3. Active change trend of water level of the Yueyaquan lake from 1947 to 2020
![]()
图 5 补水前后地下水流场对比图
Figure 5. Comparison of groundwater flow field before and after water replenishment
![]()
图 7 补水前后地下水流场对比图
Figure 7. Comparison of groundwater flow field before and after water replenishment
-
安志山, 张克存, 牛清河, 等. 敦煌鸣沙山月牙泉景区高大沙丘短期动态变化特征[J]. 干旱区研究, 2016, 33(5): 981-987 AN Zhishan, ZHANG Kecun, NIU Qinghe, et al. Short-term Dynamic Change of Mega-dunes around the Crescent Spring in Dunhuang[J]. Arid Zone Research, 2016, 33(5): 981-987.
安志山, 张克存, 王晓雷, 等. 敦煌月牙泉风沙环境及沙山动态监测[J]. 干旱区资源与环境, 2013, 27(3): 115-120 AN Zhishan, ZHANG Kecun, WANG Xiaolei, et al. Dynamic monitoring of sand hill and wind-blown sand environment in the scenic spots of Crescent Moon Spring in Dunhuang, China[J]. Journal of Arid Land Resource and Environment, 2013, 27(3): 115-120.
柴娟, 刘莎. 兴平市地下水位动态变化及影响因素分析[J]. 地下水, 2018, 40(5): 46-49 CHAI Juan, LIU Sha. Characteristics and Influencing Factors of Groundwater Level Dynamics in Xingping City[J]. Underground Water, 2018, 40(5): 46-49.
丁宏伟, 龚开诚. 敦煌月牙泉湖水位持续下降原因及对策分析[J]. 水文地质工程地质, 2004, 31(6): 74-77 DING Hongwei, GONG Kaicheng. Analyses of the reasons and countermeasures for the decline in the water level of the Crescent Spring Lake near Dunhuang[J]. Hydrogeology and Engineering Geology, 2004, 31(6): 74-77.
韩积斌, 许建新, 徐凯, 等. 柴达木盆地尕斯库勒盐湖地表水—地下水的转化与铀的补给通量[J]. 湖泊科学, 2019, 31(06): 1738-1748 doi: 10.18307/2019.0609 HAN Jibin, XU Jianxin, XU Kai, et al. The exchange relationship of surface water-groundwater and uranium flux in the Gas Hure Salt Lake of northwest Qaidam Basin, China[J]. Journal of Lake Sciences, 2019, 31(06): 1738-1748. doi: 10.18307/2019.0609
冀钦, 杨建平, 陈虹举. 1961-2015年青藏高原降水量变化综合分析[J]. 冰川冻土, 2018, 40(6): 1090-1099 JI Qin, YANG Jianping, CHEN Hongju. Comprehensive analysis of the precipitation changes over the Tibetan Plateau during 1961~2015[J]. Journal of Glaciology and Geocryology, 2018, 40(6): 1090-1099.
李平平, 王晓丹, 缪云腾, 等. 敦煌月牙泉湖近百年水位变化及其原因分析[J]. 地质论评, 2020, 66(6): 1619-1625. LI Pingping, WANG Xiaodan, MIAO Yunteng, et al. The Change of water level and its causes in Crescent Lake in Dunhuang area in the last 100 years[J]. Geological Review, 2020, 66(6): 1619-1625.
李平平. 甘肃省地下水超采区地面沉降控制区判定方法和结果探讨[J]. 中国农村水利水电, 2019, (6): 74-77 LI Pingping. Determination Method and Discussion of Ground Subsidence Control Area of Groundwater Overmining Area in Gansu Province[J]. China Rural Water and Hydropower, 2019, (6): 74-77.
黎涛, 杨俊仓, 苏春丽, 等. 敦煌月牙泉域沉积环境及泉湖水化学成因分析[J]. 干旱区地理, 2013, 36(5): 812-817 LI Tao, YANG Juncang, SU Chunli, et al. Strata sedimentary environment and hydrochemistry cause of Lake Grescent Spring in Dunhuang[J]. Arid Land Geography, 2013, 36(5): 812-817.
吕晓立, 刘景涛, 周冰, 等. 新疆塔城盆地地下水中铁锰分布特征及人类活动的影响[J]. 中国地质, 2020, 47(6): 1765-1775 LV Xiaoli, LIU Jingtao, ZHOU Bing, et al. Fe and Mn distribution of groundwater and impact of human activities in the Tacheng Basin, Xinjiang[J]. Geology in China, 2020, 47(6): 1765-1775.
祁泽学, 汪生斌, 王万平, 等. 气格尔木河冲洪积平原地下水开采潜力分析[J]. 人民黄河, 2018, 40(06): 66-71 doi: 10.3969/j.issn.1000-1379.2018.06.015 QI Zexue, WANG Shengbin, WANG Wanping, et al. Analysis of Sustainable Groundwater Resources Development Scenarios in the Golmud Alluvial Plain[J]. Yellow River, 2018, 40(06): 66-71. doi: 10.3969/j.issn.1000-1379.2018.06.015
桑学锋, 张明泉, 王浩, 等. 敦煌盆地地下水数值模拟及可视化与管理[J]. 兰州大学学报: 自然科学版, 2007, 43(3): 18-22 SANG Xuefeng, ZHANG Mingquan, WANG Hao, et al. Visual simulation and management of groundwater in Dunhuang Basin[J]. Journal of Lanzhou University: Natural Sciences, 2007, 43(3): 18-22.
施锦, 王建红. 机井开采对月牙泉水位的影响及应对措施[J]. 地下水, 2014, (6): 57-58 SHI Jin, WANG Jianhong. Groundwater Exploitation Influence on Crescent Lake Water Level and Its Countermeasures[J]. Underground Water, 2014, (6): 57-58.
孙显科, 吕亚军, 张大治, 等. 风成沙地地形1/10定律的研究与敦煌鸣沙山成因的猜想[J]. 中国沙漠, 2006, 26(5): 704-710 SUN Xianke, LV Yajun, ZHANG Dazhi, et al. Probe into the "one-tenth" Law of Aeolian Sandy Landform and Guess at Causes of Sand-dinging Mountain in Dunhuang[J]. Journal of Desert Research, 2006, 26(5): 704-710.
许朋柱, 秦伯强. 太湖湖滨带生态系统退化原因以及恢复与重建设想[J]. 水资源保护, 2002, 18(3): 31-36 doi: 10.3969/j.issn.1004-6933.2002.03.010 XU Pengzhu, QIN Boqiang. Degeneration of ecosystem of lakeside zone around Taihu Lake and planning for its rehabilitation[J]. Water Resources Protection, 2002, 18(3): 31-36. doi: 10.3969/j.issn.1004-6933.2002.03.010
杨俊仓, 张川. 月牙泉渗流场数学模型及治理方案研究[J]. 西北水资源与水工程, 2003, 14(3): 25-28 YANG Juncang, ZHANG Chuan. The selection of leakingflow fields mathematical model and harnessing program in Crescent Spring[J]. Northwest Water Resources & Water Engineering, 2003, 14(3): 25-28.
尹念文, 魏玉涛. 月牙泉的成因分析[J]. 地下水, 2010, 32(2): 20-22 YIN Nianwen, WEI Yutao. The analysis about the causes of crescent lake[J]. Ground Water, 2010, 32(2): 20-22.
袁国映, 赵子允. 楼兰古城的兴衰及其与环境变化的关系[J]. 干旱区地理, 1997, 20(3): 7-12 YUAN Guoyang, ZHAO Ziyun. Relations between the rise and the decay of the ancient city, Loulan and its environmental change[J]. Arid Land Geography, 1997, 20(3): 7-12.
岳峰, 董霁红, 温秀琴. 敦煌市月牙泉景观衰退的原因及对策[J]. 水土保持研究, 2007, 14(2): 200-202+206 YUE Feng, DONG Wenhong, WEN Xiuqin. The Decadent Reason and Countermeasure of Crescent Spring View of Dunhuang City[J]. Research of Soil and Water Conservation, 2007, 14(2): 200-202+206.
张晨, 来世玉, 高学平, 等. 气候变化对湖库水环境的潜在影响研究进展[J]. 湖泊科学, 2016, 28(4): 691-700 doi: 10.18307/2016.0401 ZHANG Chen, LAI Shiyu, GAO Xueping, et al. A review of the potential impacts of climate change on water environment in lakes and reservoirs[J]. Journal of Lake Sciences, 2016, 28(4): 691-700. doi: 10.18307/2016.0401
张号, 张克存, 安志山, 等. 敦煌月牙泉景区沙丘粒度特征[J]. 中国沙漠, 2014, 34(3): 645-649. ZHANG Hao, ZHANG Kecun, AN Zhishan, et al. Grain size characteristics of sand dunes around the Crescent Moon Spring in Dunhuang, China[J]. Journal of Desert Research, 2014, 34(3): 645-639.
张克存, 牛清河, 屈建军, 等. 敦煌鸣沙山月牙泉景区风沙环境分析[J]. 中国沙漠, 2012, 32(4): 896-900 ZHANG Kecun, NIU Qinghe, QU Jianjun, et al. Analysis of Wind-blown Sand Environment in the Singing Sand Mountain & Crescent Moon Spring Scenic Spot in Dunhuang, China[J]. Journal of Desert Research, 2012, 32(4): 896-900.
张明泉, 曾正中, 蔡红霞, 等. 敦煌月牙泉水环境退化与防治对策[J]. 兰州大学学报(自然科学版), 2004, 40(03): 99-102 ZHANG Mingquan, ZENG Zhengzhong, CAI Hongxia, et al. Astudies of degradation of water system and some control measures in the Moon Lake of Dunhuang Region[J]. Journal of Lanzhou University (Natural Sciences), 2004, 40(03): 99-102.
张文化, 魏晓妹, 李彦刚. 气候变化与人类活动对石羊河流域地下水动态变化的影响[J]. 水土保持研究, 2009, 16(01): 183-187 ZHANG Wenhua, WEI Xiaomei, LI Yangang. Groundwater Dynamic Evolution Under Climatic Change and Human Activites in Shiyang River Basin[J]. Research of soil and water conservation, 2009, 16(01): 183-187.
朱亮, 刘景涛, 杨明楠, 等. 1998年以来兰州市地下水环境变化及驱动因素[J]. 中国地质, 2020, 47(6): 1677-1687 ZHU Liang, LIU Jingtao, YANG Mingnan, et al. Changes and driving factors of groundwater environment in Lanzhou since 1998[J]. Geology in China, 2020, 47(6): 1677-1687.
Garth V D K, Dwayne K, Marlene S E. Long-Term Water Level Changes in Closed-Basin Lakes of the Canadian Prairies[J]. Canadian Water Resources Journal, 2008, 1(1): 23-38.
Katsifarakis K L. Groundwater Pumping Cost Minimizationan-an Analytical Approach[J]. Water Resources Management, 2008, 8(8): 1089-1099.
Lan Yingying, Jin Menggui, Yan Chun, et al. Schemes of groundwater exploitation for emergency water supply and their environmental impacts on Jiujiang City, China[J]. Environmental Earth Sciences, 2015, 5(5): 2365-2376.
Li Zhi, Liu Wenzhao, Zhang Xunchang, et al. Impacts of land use change and climate variability on hydrology in an agricultural catchment on the Loess Plateau of China[J]. Journal of Hydrology, 2009, (1-2): 35-42.
Phan S V, Leavitt P R, Mcgowan S. Spatial variability of climate and land-use effects on lakes of the northern Great Plains[J]. Limnology and Oceanography, 2008, 2(2): 728-742.
Gaur S, Mimoun D, Graillot D. Advantages of the analytic element method for the solution of groundwater management problems[J]. Hydrological Processes, 2011, 25(22): 3426−3436.
TU J. Combined impact of climate and land use changes on streamflow and water quality in eastern Massachusetts, USA[J]. Journal of Hydrology, 2009, (3-4): 268-283.
Zhu Lin, Gong Huili, Li Xiaojuan, et al. Land subsidence due to groundwater withdrawal in the northern Beijing plain, China[J]. Engineering Geology, 2015, 193(02): 243-255.
下载:
计量
- 文章访问数: 173
- HTML全文浏览量: 34
- PDF下载量: 122
