水源热泵井群布置方案对地下水流场影响分析

以某地下水源热泵系统工程为例,依据工程场地实际抽水-回灌试验数据,建立研究区地下水流数值模拟模型,分析该地区适宜的井群布置方案,并定量研究不同井群布置方案对地下水流场的影响范围与程度,为确定合理的井群布置方案提供技术支撑。研究表明：对于同一含水层回灌,随着抽水-回灌时间的延长,就有利于抽水-回灌能力及尽量减小对周边地下水开采的影响而言,抽水井、回灌井交叉布置,部分回灌井集中布置的方案较为合理。     

地形地貌对地震波放大效应数值模拟研究

采用FLAC3D 软件进行大量数值模拟,研究了地形地貌对地震波放大效应规律。结果表明:对于单面坡模型,斜坡较低时,坡体对天然地震波和简谐波的峰值加速度在竖直向和顺坡面向产生线性放大效应。斜坡较高时,若坡角低,两种波形的加速度等值线近似平行于坡面,天然地震波坡顶正下方极值性不明显,但随坡角增大而明显,坡面正下方可见若干条与坡面相平行具节律性的极小值条带,随坡角增大向一条极值带转化,简谐波在高低坡脚时坡顶面正下方区域内都出现等间距相间分布具节律性的极值条带,坡面正下方的极值带状性不如天然地震波明显,随坡角增大而消失; 若坡角高,两波形加速度放大系数顺坡面呈相间的极值圈分布。双面坡坡顶宽度越小,山脊越单薄,加速度放大效应越强烈,随坡顶宽度增大,坡肩加速度放大系数趋于稳定,双面坡对加速度放大效应强于单面坡。对加速度放大效应,多峰形双峰形单峰形,坡顶越不规则,加速度放大效应越强烈。V形河谷坡面上半部分加速度放大效应要强于U形河谷,U形河谷的谷底中部对地震波动力响应要强于V形谷。坡面局部不规则性只影响到坡面局部范围内的峰值加速度分布,坡高、坡角、坡顶几何特征为整个峰值加速度分布的控制因素。     

含节理黄土渗流数值模型研究

本文建立了在考虑黄土中有节理存在时,水分场数值计算的有限元方法。采用质量守恒的观点推导了质点元中饱和度的变化与流速的关系。进而利用达西定律得到以水头为变量的渗流基本方程。针对黄土垂直节理的渗流特点,确定了节理渗流基本方程中的参数。采用四边形等参元,利用Galerkin加权余量法建立考虑节理影响的黄土非饱和渗流的有限元形式。对局部水头边界条件下的黄土节理二维水分入渗问题进行了数值分析。结果表明,节理对黄土场地湿润峰的迁移有很大影响。     

芦山地震的发生对周围断层影响的数值模拟

芦山地震发生后,地震的发生造成周围断层应力变化值得关注。本文基于川滇地区的三维非线性有限元模型,利用芦山地震同震静态滑移量结果,分析地震的发生对川滇地区主要断裂的同震加卸载效应。初步结果表明,芦山地震的发生造成龙门山断裂中南段、岷江断裂、马尔康断裂、鲜水河断裂北西段、大凉山断裂南段、小江断裂南段不同程度的应力增加。其中龙门山断裂中南段增加最为显著,最大库仑应力增加量达0.035 MPa;岷江断裂次之,最大达0.0075 MPa;马尔康断裂增加量达0.0031 MPa;鲜水河断裂北西段达0.0008 MPa。而从断裂带同震应变积累与释放方面分析的结果同样表明上述四条断裂地震危险性增强。该结果可以为川滇地区地震危险性分析提供一定的参考依据。     

Effect of time-evolving age and convergence rate of the subducting plate on the Cenozoic adakites and boninites

Partial melting of subducting oceanic crust expressed as high-Mg volcanic rocks such as adakites and boninites has been actively studied for decades, and Lee and King (2010) reported that time-evolving subduction parameters such as the age and the subduction rate of the converging oceanic plate play important roles in transient partial melting of the subducting oceanic crust (e.g., Aleutians). However, few subduction model experiments have considered time-evolving subduction parameters, posing problems for studies of transient partial melting of subducting oceanic crust in many subduction zones. Therefore, we constructed two-dimensional kinematic–dynamic subduction models for the Izu–Bonin, Mariana, Northeast Japan, Kuril, Tonga, Java–Sunda, and Aleutian subduction zones that account for the last 50 Myr of their evolution. The models include the time-evolving age and convergence rate of the incoming oceanic plate, so the effect of time-evolving subduction parameters on transient partial melting of oceanic crust can be evaluated. Our model calculations revealed that adakites and boninites in the Izu–Bonin and Aleutian subduction zones resulted from transient partial melting of oceanic crust. However, the steady-state subduction model using current subduction parameters did not produce any partial melting of oceanic crust in the aforementioned subduction zones, indicating that time-evolving subduction parameters are crucial for modeling transient eruption of adakites and boninites. Our model calculations confirm that other geological processes such as forearc extension, back-arc opening, mantle plumes and ridge subduction are required for partial melting of the oceanic crust in the Mariana, Northeast Japan, Tonga, and southeastern Java–Sunda subduction zones.     

降雨入渗补给地下水数值模拟研究

针对在地下水资源评价中入渗补给量采用经验估算存在较大误差的问题,以神东矿区为例,在广泛调查矿区包气带岩性结构并结合野外取样、室内参数测定的基础上,采用数值模拟的方法建立包气带水分运移数值模拟模型,定量模拟矿区内不同地段降雨入渗强度,探讨影响降雨入渗强度的主要因素,计算得出研究区降雨入渗系数大致在0.18~0.27,分析认为影响降雨入渗强度的因素有降雨量、潜水埋深、包气带岩性等。其中在研究区广泛分布的风积沙对地下水资源起到了一定的保护作用。     

王坡井田煤层气井部署方法研究

王坡煤矿是高瓦斯矿井,已有的井下瓦斯抽排不能满足矿井安全生产的要求,迫切需要通过地面抽采来缓解。在已有地质成果基础上,采用煤层气藏数值模拟方法,对王坡井田地面煤层气抽采工程部署方案进行优化。结果表明:A区块是井田煤层气开发试验首选区,以梅花形井网作为试验井网,井间距300m×250m,井网长轴方向60°。这种布井方式在工程实践中取得了较好效果,具有较大的推广价值。     

Improvement of groundwater resources potential by artificial recharge technique: A case study of Charf El Akab aquifer in the Tangier region,Morocco

In arid and semi-arid zones, water is the most vulnerable resource to climate change. In fact, various techniques such as artificial recharge are adopted to restore aquifers and to ensure aquifer sustainability in relation to the accelerated pace of exploitation. Morocco is a Mediterranean country highly vulnerable to climate change, many of its main aquifers are subjected to excessive drawdowns. This technique is practiced to increase potentiality of these aquifers. In the Northwestern area of Morocco, the significant development experienced by Tangier City in the industrial, tourism, and commercial sectors will lead to increased water requirements-up to 5 067 L/s (159.8 mm3) by 2030. However, the Charf El Akab aquifer system, subject to artificial recharge, is the only groundwater resource of Tangier region; hence, a rational management context is needed to ensure aquifer sustainability, and optimized exploitation under the background of differing constraints, such as increased water requirements, and climate change impacts. This work aims to respond, for the first time, to the Charf El Akab aquifer overexploitation problem, and to evaluate the future scenarios of its exploitation in the event of failure of one of the superficial resources. This work also presents a synthesized hydrodynamic modeling based on the results of the numerical simulations carried out using Feflow software for 2004 (date of cessation of injections) and 2011 (date of resumption of these facilities), making it possible to evaluate the impact of the artificial recharge on the piezometric level of the aquifer on a spatiotemporal scale. Finally, the exploitation scenarios have shown that the aquifer of Charf El Akab will not adequatly provide for the region's water requirements on the future horizon, entailing an optimal management of water resources in the region and an intentionally increased recharge rate.     

Proposed nonlinear 3-D analytical method for piled raft foundations

The load distribution and deformation of piled raft foundations subjected to axial and lateral loads were investigated by a numerical analysis and field case studies. Special attention is given to the improved analytical method (YSPR) proposed by considering raft flexibility and soil nonlinearity. A load transfer approach using p–y, t–z and q–z curves is used for the analysis of piles. An analytical method of the soil–structure interaction is developed by taking into account the soil spring coupling effects based on the Filonenko-Borodich model. The proposed method has been verified by comparing the results with other numerical methods and field case studies on piled raft. Through comparative studies, it is found that the proposed method in the present study is in good agreement with general trend observed by field measurements and, thus, represents a significant improvement in the prediction of piled raft load sharing and settlement behavior.     

柴达木盆地北缘东段石炭纪至白垩纪构造演化 模式及地层缺失原因探讨

为了明确柴达木盆地北缘东段石炭纪至白垩纪盆地演化历史,分析石炭系—侏罗系缺失的原因,本文重新确定了石炭系和侏罗系残余地层分布范围,并根据平衡剖面恢复和断裂落差计算,分析了主要断裂的活动特征,提出了柴北缘东段石炭纪至白垩纪的构造演化模式。应用有限单元法模拟了柴北缘东段印支期(三叠纪)、燕山早期(早—中侏罗世)、燕山晚期(晚侏罗世—白垩纪)的应力场,对构造演化模式加以验证。研究结果表明:柴北缘东段自石炭纪至白垩纪经历了石炭纪—二叠纪伸展、三叠纪挤压褶皱、早—中侏罗世断陷、晚侏罗世—白垩纪早期挤压坳陷和白垩纪末挤压反转五个构造演化阶段。三叠纪,柴北缘东段在印支期发育两排近东西走向的背斜凸起,造成石炭系—二叠系在各地区遭受不同程度的剥蚀;侏罗纪—白垩纪早期,欧南地区为继承性隆起区,未完全接受沉积;白垩纪末,受燕山晚期旋回影响,构造反转,逆冲断裂复活,绿梁山、锡铁山、埃姆尼克山、欧隆布鲁克山等主要山体隆升,遭受剥蚀。