基于NAD算法的声波方程时间四阶差分解法

波动方程数值模拟是研究地震波传播机理的重要工具,有限差分求解波动方程是当前地震波数值模拟的主要方法之一。当地下介质中的地震波速度较低或地震波高频成分丰富时,常规有限差分技术常常产生严重的数值频散误差,这种误差会降低数值模拟的精度,影响对地震波传播机理的分析。为压制地震波数值模拟时产生的数值频散误差,提高波场模拟精度,提出了基于NAD算子的时间四阶精度波动方程差分格式。根据对应的差分格式,分析了该差分格式的数值频散关系。与常规四阶精度差分算法的频散曲线相比,基于NAD时间四阶精度差分方法不但能够实现时间频散的有效压制,同时其基于更多网格点的位移分量和位移梯度分量空间微分求解方法还能够实现空间频散的有效压制。另外在相同模型条件下,基于NAD算法的声波方程时间四阶差分解法可采用大网格对模拟空间进行差分离散,减少网格数,提高计算效率。     

Characterization of Continental Intercalaire aquifer (CI) in the Tinrhert-East Area-Illizi Basin on the Algerian-Libyan Border

This paper addresses the characterization of the Continental Intercalaire aquifer(CI) in the Tinrhert-East area of Illizi Basin on the Algerian-Libyan border, which belongs to the SASS1 system, one of the biggest transboundary aquifers in the world. This study concerns a superficies of 4 300 km~2. On the basis of Mud Logging borehole data conducted in this part of the aquifer, a realistic characterization of the aquifer was done. The thickness of the CI aquifer varies from 300 m in the south to 700 m in the north, and the depth ranges from 180 m to 320 m. The interpretation of the logs showed that the aquifer is characterized by a maximum net thickness in its southwestern part(more than 600 m), the porosity is very high, ranging from 30% in the west to 24% at the Libyan borders, the permeability is low to medium around 10~(-5) m/s, and the maximum transmissivity values of about 8× 10~(-3) m~2/s were recorded at the center of the study area. The depth of water varies from 235 m to 312 m, and the water flows from south to north, in accordance with the general direction observed in the CI aquifer in the Northern Sahara Aquifer System(SASS). The porosity values obtained from the interpretation of the sonic and density logs permit to estimate the water reserves of this aquifer considered fossil, at thresholds much higher than what was considered until now.     

南海北部不同涡旋对内潮的影响

南海北部吕宋海峡是内潮最为活跃的区域之一,且涡旋种类繁多,不同特性的涡旋对内潮的影响不同。基于近岸与区域海洋共同模式(coastal and regional ocean community model,CROCO),模拟探究理想涡旋存在时,涡旋位置、极性、峰值流速和半径对内潮的影响。结果表明:涡旋位置是影响内潮的直接因素,位于涡旋区域内的内潮是主要影响对象,涡旋中心以西内潮方向变化的角度是以东的3倍。气旋涡和反气旋涡分别使潮能通量的方向向南和向北偏转,最大偏转角度超过12°,当涡旋所致背景流与内潮传播方向一致时,内潮群速度增强,反之减弱。涡旋对内潮的影响范围和幅度随着涡旋的半径和峰值流速的增大而变大。当涡旋峰值速度变大时,反气旋涡心以北的潮能通量增长量超过15 kW/m。当涡旋半径增大时,涡旋峰值速度的位置发生变化,涡旋的峰值流速和半径共同影响潮能通量水平分布结构,使其呈现纬向单峰或多峰结构。     

华北地区"12·7"降雪过程的数值模拟研究

对2001年12月7日一次引发北京交通堵塞的降雪过程成因作了模拟研究.模拟结果显示PSU/NCAR的MM5有可能模拟出此次北方较弱的降雪过程,模拟的降雪量、落区以及持续时间与观测较一致.在成功模拟的基础上,利用模式输出的时空分辨率较高的资料,对此次降雪的发生、发展和水汽输送过程等进行了分析,结果表明:(1)降雪发生前对流层中层先出现饱和,而低层并未饱和,这种弱降雪的产生似乎并不需要中低层有深厚的湿层存在;(2)此次降雪是由对流层中层快速移动的短波槽和近地面出海高压后部的回流共同影响的结果,近地面的高压回流主要对增加低层的湿度有贡献,槽前的西南气流将水汽由南向北输送到华北地区,辐合引起的上升运动又将水汽输送到对流层中上层,槽前的水汽输送和辐合上升是此次降雪过程的触发机制之一;(3)华北地区大气中可降水量达到7 mm以上时,就可能有弱降雪发生,并有可能根据可降水量判断降雪的维持时间;(4)冰相云物理过程对成功地模拟降雪是不可忽视的.     

Fully coupled hydraulic fracture simulation using the improved element partition method

The improved element partition method (IEPM) is a newly developed fracture simulation approach. IEPM allows a fracture to run across an element without introducing extra degrees of freedom. It can also simulate any number of fractures in a prescribed mesh without remeshing. In this study, the IEPM is extended to hydraulic fracture simulation. First, the seepage and volumetric storage matrix of a cracked element are derived using virtual nodes (the intersection points of a crack with element edges). Subsequently, the fully coupled hydromechanical equation is derived for this cracked element. To eliminate the extra degrees of freedom (virtual nodal quantities), the water pressure and displacement of the virtual nodes are associated with their adjacent nodes through least squares interpolation. Finally, the fully coupled equation in terms of nodal quantities is obtained. The verification cases validate the method. By using this method, the field-scale hydraulic fracturing process is well simulated. The proposed approach is simple and efficient for field-scale hydraulic fracture simulation.     

Numerical study for vegetation effects on coastal wave propagation by using nonlinear Boussinesq model

The vegetation has important impacts on coastal wave propagation. In the paper, the sensitivities of coastal wave attenuation due to vegetation to incident wave height, wave period and water depth, as well as vegetation configurations are numerically studied by using the fully nonlinear Boussinesq model. The model is based on the implementation of drag resistances due to vegetation in the fully nonlinear Boussinesq equation where the drag resistance is provided by the Morison’s formulation for rigid structure induced drag stresses. The model is firstly validated by comparing with the experimental results for wave propagation in vegetation zones. Subsequently, the model is used to simulate waves with different height, period propagating on vegetation zones with different water depth and vegetation configurations. The sensitivities of wave attenuation to incident wave height, wave period, water depth, as well as vegetation configurations are investigated based on the numerical results. The numerical results indicate that wave height attenuation due to vegetation is sensitive to incident wave height, wave period, water depth, as well as vegetation configurations, and attenuation ratio of wave height is increased monotonically with increases of incident wave height and decreases of water depth, while it is complex for wave period. Moreover, more vegetation segments can strengthen the interaction of vegetation and wave in a certain range.     

How soil texture,channel shape and cross-sectional area affect moisture dynamics and water loss in irrigation channels

To enhance the utilization efficiency of farmland irrigation water and reduce the leakage of water conveyance channels, the leakage process of channels was simulated dynamically. The simulated results were compared with data measured in laboratory experiments, and the performance of the model was evaluated. The results indicated that the simulated values of the model were consistent with the observation values, and the R² values varied between 0.91 and 0.99. In addition, based on the laboratory experiments, a water supply system (Mariotte bottles) and soil box were built using plexiglass. Three influencing factors, namely, the channel form, soil texture and channel cross-sectional area, were varied to observe and calculate the resulting cumulative infiltration amount, infiltration rate and wetting front migration distance. HYDRUS-3D software was used to solve the three-dimensional soil water movement equation under different initial conditions. The results demonstrated that the U-shaped channel was more effective than the trapezoidal channel in increasing the utilization efficiency of the water resources. A U-shaped channel with a small channel cross-sectional area should be adopted and the soil particle size should be prioritized in the construction of water conveyance channels for farmlands. The simulation results were in agreement with the observed results, which indicates that HYDRUS-3D is a reliable tool that can accurately simulate the soil moisture movement in water conveyance channels. The research results can provide a reference for the design and operation of farmland irrigation systems.     

A two-layer model for studying 2D dissolved pollutant runoff over impermeable surfaces

Dissolved pollutants in stormwater are a main contributor to water pollution in urban environments. However, many existing transport models are semi-empirical and only consider one-dimensional flows, which limit their predictive capacity. Combining the shallow water and the advection–diffusion equations, a two-dimensional physically based model is developed for dissolved pollutant transport by adopting the concept of a ‘control layer’. A series of laboratory experiments has been conducted to validate the proposed model, taking into account the effects of buildings and intermittent rainfalls. The predictions are found to be in good agreement with experimental observations, which supports the assumption that the depth of the control layer is constant. Based on the validated model, a parametric study is conducted, focusing on the characteristics of the pollutant distribution and transport rate over the depth. The hyetograph, including the intensity, duration and intermittency, of rainfall event has a significant influence on the pollutant transport rates. The depth of the control layer, rainfall intensity, surface roughness and area length are dominant factors that affect the dissolved pollutant transport. Finally, several perspectives of the new pollutant transport model are discussed. This study contributes to an in-depth understanding of the dissolved pollutant transport processes on impermeable surfaces and urban stormwater management.     

Assessment of 2DH and pseudo‐3D modelling platforms in a large saline aquatic system: Lake Urmia,Iran

The main objective of this paper is to provide comparative quantitative examinations on the capabilities of two‐dimensional horizontal and pseudo‐three‐dimensional (3D) modelling approaches for simulating spatial and temporal variability of the flow and salinity in Lake Urmia, Iran. The water quality in the lake has been an environmentally important subject partly because this shallow hypersaline aquatic ecosystem is considered to be one of the largest natural habitats of a unique multicellular organism, Artemia urmiana. This brine shrimp is the major food source for many of the protected and rare shorebirds that visit the lake. A. urmiana can grow and survive in certain ranges of salinity, and their disappearance could lead to an alteration of existing equilibria. The lake has also experienced considerable man‐made changes during the past three decades. A newly built crossing embankment almost divided the lake into two northern and southern halves. A relatively small opening of 1.25 km in the new embankment provides water connections between the two halves. As a result, the flow and salinity regimes have been significantly changed. This might have had adverse serious impacts on the lake ecosystem. In the current study, the two‐dimensional horizontal hydrodynamic model has been found to provide reasonable predictions for the flow regime in the lake, whereas its salinity predictions have not been consistent with the field observations. The pseudo‐3D model has produced results fairly close to the salinity measurements and its temporal and spatial variations. The pseudo‐3D model has been used for evaluating the embankment effects on the lake hydrodynamics and on the salinity conditions. The effectiveness of introducing a different number or length of openings in the embankment for restoring the pre‐embankment conditions has also been examined. These remedy options have been found not to offer substantial improvements to the lake's existing ecosystem. Copyright © 2013 John Wiley & Sons, Ltd.     

Simulation of drilling‐induced compaction bands using discrete element method

Rock failure is observed around boreholes often with certain types of failure zones, which are called breakouts. Laboratory‐scale drilling tests in some high‐porosity quartz‐rich sandstone have shown breakouts in the form of narrow localized compacted zones in the minimum horizontal stress direction. They are called fracture‐like breakouts. Such compaction bands may affect hydrocarbon extraction by forming barriers that inhibit fluid flow and may also be a source of sand production. This paper presents the results of numerical simulations of borehole breakouts using 3D discrete element method to investigate the mechanism of the fracture‐like breakouts and to identify the role of far‐field stresses on the breakout dimensions. The numerical tool was first verified against analytical solutions. It was then utilized to investigate the failure mechanism and breakout geometry for drilled cubic rock samples of Castlegate sandstone subjected to different pre‐existing far‐field stresses. Results show that failure occurs in the zones of the highest concentration of tangential stress around the borehole. It is concluded that fracture‐like breakout develops as a result of a nondilatant failure mechanism consisting of localized grain debonding and repacking and grain crushing that lead to the formation of a compaction band in the minimum horizontal stress direction. In addition, it is found that the length of fracture‐like breakouts depends on both the mean stress and stress anisotropy. However, the width of the breakout is not significantly changed by the far‐field stresses. Copyright © 2013 John Wiley & Sons, Ltd.