Numerical investigation of hydraulic characteristics and prediction of cavitation number in Shahid Madani Dam's Spillway

The prediction of the probability of cavitation occurrence to prevent serious damages in the spillways is the major concern for hydraulic engineers. In this research, the three-dimensional model of Shahid Madani Dam’s spillway was simulated with the Flow 3D software and by the comparison of numerical model results with the experi-mental data, the probability of occurrence of the cavitation phenomenon has been investigated. The flow parameters including pressure, velocity, and water depth were calculated for three different flow rates of 495 m3/s, 705 m3/s and 2 205 m3/s respectively. The Renormalization Group (RNG) turbulence model was used to simulate current turbulence. Comparison of simulation results for pressure, velocity and water depth with the results of the experimental model with two statistical indices Root Mean Square Error (RMSE) and Coefficient of Determination (R2) showed that the numerical simulation results are in good consistency with experimental model. However, simulation results indicated that at any flow rates with a return period of 1 000 years, probable maximum flood and designed flow rates, the cavitation number is not lower than the critical cavitation number; Therefore, it is predicted that the cavitation phenomenon in Shahid Madani Dam’s spillway will not happen.     

Discharge prediction using hydraulic characteristics of mean velocity equation

Discharge is an important factor in river design for water utilization, water control and hydraulic structures; therefore, an accurate estimation of the discharge is required. At present, a rating curve depicting the relationship between a stage and discharge is used to calculate the discharge from river systems. Although the rating curve has an advantage in that it can predict and use the discharge during the flood season in which the measurement is difficult, there is room for improvement as it does not reflect the hydraulic characteristics of rivers. Therefore, in this study, discharge was predicted using the convenient calculation method with empirical mediating variables of the Manning and Chezy equations which were proposed by the author’s previous research as a new methodology for estimating discharge in an open channel. This was proven, based on the data measured in a meandering open channel system in a lab at the Mississippi River in the US and at the Columbia Del Dique Canal, and an accuracy level at a coefficient of 0.8 was demonstrated. Thus, this method, which reflects the hydraulic characteristics and predicts the discharge in a simple manner, is expected to be convenient in practice.     

Numerical investigation of hydraulic fracture network propagation in naturally fractured shale formations

Hydraulic fracture network (HFN) propagation in naturally fractured shale formations is investigated numerically using a 3D complex fracturing model based on the discrete element method. To account for the plastic deformation behavior of shales, the Drucker–Prager plasticity model is incorporated into the fracturing model. Parametric studies are then conducted for different Young's moduli, horizontal differential stresses, natural fracture (NF) properties, injection rates, and number and spacing of perforation clusters. Numerical results show that horizontal differential stress primarily determines the generation of a complex HFN. The plastic deformation of shale can reduce the stimulated reservoir volume; this is more obvious with Young's modulus of less than 20 GPa. In addition, a higher injection rate could largely increase the fracture complexity index (FCI). Moreover, increasing perforation cluster numbers per fracturing stage is beneficial for increasing the FCI, but it also increases the potential merging of neighboring fractures, which may lead to non-uniform development of HFN in far-wellbore regions. To achieve uniform development of HFN within a fracturing stage, the distribution of NFs should be fully considered. The results presented here may provide improved understanding of HFN generation and are favorable for optimizing fracturing treatment designs for shale formations.     

高坝明流泄洪洞掺气减蚀三维数值模拟分析

高坝明流泄洪洞的掺气减蚀问题倍受关注。以溪洛渡水电站3#泄洪洞为研究对象,结合大比尺单体模型试验成果,应用分段计算方法对泄洪洞掺气减蚀问题进行了三维紊流数值模拟分析。结果表明:泄洪洞各道掺气坎布置合理且有稳定掺气空腔,掺气坎体型基本合理,明流段洞顶余幅足够。泄洪洞反弧段及末端一段区域水力特性复杂,应重点关注此部位的设计及施工问题。对于高流速、大梯度泄洪洞水力学问题,采用分段计算能提高计算速度。     

利用探地雷达探测地下空洞

简要介绍了探地雷达工作原理及其探测地下空洞的方法技术，以广州某地基勘察中地下空洞的检测为例，分析了地下空洞的雷达波组特征，总结了探地雷达时深剖面图的异常特征。指出该方法具有速度快、探测精度高、可获得连续结果等特点。经开挖验证，探地雷达异常特征与地下空洞实际位置相吻合。     

根系和裂隙对土体水力和力学特性影响数值模拟

为探究根系和裂隙对土体水力及力学特性的影响,利用有限元软件计算降雨过程中裂隙和根土间隙对渗流场的影响,并以渗流计算结果为基础,分析降雨对根系固土作用的影响,采用分区强度折减法对降雨前后根?土复合体进行直剪试验模拟,同时考虑了侧根倾角的影响。结果表明,裂隙和根土间隙为雨水入渗提供优先通道,降雨影响深度随裂隙深度增加而增加;有根系时降雨影响深度由主根深度决定,侧根倾角对其影响较小,考虑根土间隙影响的降雨影响深度相较于无根系工况增加了93.3%;根系能显著提高土体的抗剪强度,相对于垂直主根方向的不同倾角,侧根增加土体抗剪强度由大至小依次为60o、45o、30o侧根和无根系;雨水入渗降低了土体强度,同时削弱了根系固土作用,使得降雨后根-土复合体抗剪强度大幅降低,是许多植被覆盖边坡仍发生浅层滑动原因之一。     

Numerical prediction of liquefied ground characteristics from back-analysis of lateral spreading centrifuge experiments

In this paper, liquefied and laterally spreading soils triggered by seismic shaking are modeled as viscoplastic Bingham media characterized by two rheological parameters: the undrained residual shear strength and the Bingham viscosity coefficient. Since the precise evaluation of these two characteristics directly by appropriate in situ or laboratory experimental tests remains a very difficult task, an identification procedure is developed to assign numerically realistic values to both rheological characteristics from back-calculation of liquefaction-induced lateral spreading using centrifuge experiments. The proposed numerical procedure is applied successfully to two series of reported centrifuge tests where lateral displacements data during shaking were available.     

Numerical investigation of slippage characteristics of normal and reverse faults under fluid injection and production

Nowadays, a great deal of petroleum geology and engineering projects associated with underground fluid injection and production (FIP) are widely conducted around the world. The FIP engineering may cause complex stress perturbation and trigger seismicity, which have been extensively reported and studied. In this paper, we investigated the fault slippage characteristics influenced by FIP. It reveals that for a fault (normal or reverse) that penetrating through the reservoir into the caprock and underburden, the footwall reservoir is the relatively stable one for fluid injection in a fluid FIP engineering. No matter it is a normal or a reverse fault, injecting fluid into footwall reservoir and producing fluid from the hanging wall reservoir can induce smaller fault slippage. After having determined the better fluid injection–production pattern, we studied the influence of three key factors of fault, i.e., dip, offset and depth, on fault stability. We found that, in our range of study, the influence of single action of fault dip, offset or depth on fault slippage in a FIP engineering was small. However, the influence of the combined effects of the factors may be large. Finally, we studied the effect of different pressure management scenarios on fault responses based on the specific fluid injection–production pattern. The results revealed that appropriate pressure management could effectively reduce fault slippage in a FIP engineering. However, inappropriate pressure management may cause much larger fault slippage. Given these concerns, it is therefore vital that the effect of pressure management scenario to be modeled prior to FIP.     

交叉裂隙水力学开度的计算及非线性水力特性研究

基于人工交叉裂隙模型,通过室内透水试验,得到了不同进口与出口组合时的流量和压力变化的关系曲线。同时结合串联裂隙水力学开度的计算公式,利用最小二乘法求解超静定方程组,计算交叉裂隙中每个裂隙单元的水力学开度。通过拓扑得到裂隙交叉点的几何尺寸,并求解Navier-Stokes方程,研究裂隙交叉对其水力特性的影响。结果表明,通过室内试验和串联裂隙水力学开度计算公式相结合的方法,可以准确计算每个裂隙单元的水力学开度。当雷诺数Re较小时,交叉点内部的流体成稳态层流流动;当雷诺数较大（比如Re≥100）时,可以观测到明显的漩涡,说明流体的惯性力远大于黏性力,经典的立方定律不再适用。对于1个进口2个出口的情况,出口的流量分配比率与雷诺数成二次函数关系,随着雷诺数的增大,流量分配比率的非线性越来越明显,其最大分配比率变化超过15%。出口的水力学开度e与初始水力学开度（即力学开度E0）的比值e/E0和雷诺数Re也具有二次函数关系。当Re10时,e/E0呈现出较强的非线性;利用该关系式可以得到修正的立方定律,从而进一步求解交叉裂隙的水力特性问题。     

Analytical investigation of hydraulic fracture-induced seismicity and fault activation

More recent public discourse has taken place regarding the potential correlation between seismic activity and hydraulic fracturing in shale gas reservoirs. Public fears about the risk of seismicity stem mainly from past earthquakes induced by conventional deep injections because the two types of projects share similar mechanisms of rock failure and fault activation. Although previous earthquake risks associated with fluid injection were not serious, the situation would be far more problematic if hydraulic fracturing in a shale gas reservoir triggered a similar-sized earthquake due to potential environmental issues. In fact, almost all documented injection-induced earthquakes have been associated with long-duration and high-volume injection rather than short-term (hours) pressurization (e.g., hydraulic fracturing). In general, hydraulic fracturing operations mostly induce microseismic events through rock failure and activation of small fractures. Although shale reservoirs in tectonically active zones pose a high risk of inducing large-magnitude seismic activities, the internal geological conditions and external stimulation conditions are impossible to be satisfied simultaneously to trigger activation of an entire fault and to result in a destructive earthquake during hydraulic fracturing operations.     

水平圆柱强迫振荡水动力特性的数值模拟

为了研究处于自由面以下完全淹没状态的水平圆柱在强迫振荡运动时的水动力特性,采用基于黏性流理论建立的二维两相流数值波浪水槽模型,对不同液相黏性条件下强迫振荡水平圆柱的受力进行计算,并对压力、黏性切力和圆柱运动之间的相位关系特征进行对比和分析,进而结合流场分析解释黏性影响机理。结果表明:黏性切力和涡旋压力对流体作用力的贡献差别是导致不同流体黏性下流体作用力结果差异的主要原因;涡旋运动相对圆柱振荡运动的滞后性受流体黏性影响显著,导致不同流体黏性下压力之间有相位差;流体水质点相对于圆柱的滞后运动在大黏性流体中更为显著,这导致了其黏性切力的相位超前现象。     

水平圆柱强迫振荡水动力特性的数值模拟

为了研究处于自由面以下完全淹没状态的水平圆柱在强迫振荡运动时的水动力特性,采用基于黏性流理论建立的二维两相流数值波浪水槽模型,对不同液相黏性条件下强迫振荡水平圆柱的受力进行计算,并对压力、黏性切力和圆柱运动之间的相位关系特征进行对比和分析,进而结合流场分析解释黏性影响机理。结果表明:黏性切力和涡旋压力对流体作用力的贡献差别是导致不同流体黏性下流体作用力结果差异的主要原因;涡旋运动相对圆柱振荡运动的滞后性受流体黏性影响显著,导致不同流体黏性下压力之间有相位差;流体水质点相对于圆柱的滞后运动在大黏性流体中更为显著,这导致了其黏性切力的相位超前现象。     

坡沟系统坡面径流流态及水力学参数特征研究

利用变坡度坡沟系统概化模型和人工模拟降雨试验,定量分析了在60、90和130mm/h降雨强度下坡沟系统坡面径流流态及水力学参数特征。结果表明,上方汇水和降雨强度的增大使坡沟系统水流雷诺数和弗劳德数呈明显增大,水流流态由缓流演变为急流,坡面水流阻力系数明显减小,从而使坡沟侵蚀产沙量显著增大。     

砂砾岩水力压裂裂缝扩展规律的数值模拟分析

砂砾岩储层一般具有岩性和渗透性变化大、孔隙度低、连通性差、孔隙结构复杂和非均匀性严重等特点,因此,在水力压裂过程中,裂缝扩展形态难以控制,大规模改造难度大。针对国内某典型砂砾岩油藏特征,采用数值计算方法对砂砾岩压裂裂缝的扩展规律进行了研究,包括地应力场、砾石含量和粒径等对裂缝扩展形态及压裂压力的影响。研究表明,砾石的存在增加了压裂裂纹扩展的复杂性,裂纹主要有止裂、偏转、穿透和吸附4种表现模式,但主应力差严格控制着裂纹的走向,随着主应力差的增大,裂纹由总体绕砾扩展转变为总体穿砾扩展,失稳压力随着主应力差的增大而明显减小;砾石含量的多少体现了砂砾岩试样宏观的非均匀性,含量越高均匀性越差,随着砾石含量的提高,裂纹与砾石的相互作用占据主导地位,失稳压力随砾石含量的增加而增大;当砾石体积含量一定时,砾石粒径对压裂压力的影响主要取决于砾石排列的随机性,失稳压力随砾石粒径的增大而略有增大。     

土壤饱和导水率的多尺度预测模型与转换关系

运用联合多重分形方法研究不同土层土壤饱和导水率与土壤基本物理特性的多尺度相关性,建立不同土层土壤饱和导水率的多尺度预测模型,构建不同土层土壤饱和导水率之间的转换关系。结果表明:不同土层土壤饱和导水率与土壤基本物理特性的相关程度排序不同;在单一尺度和多尺度上,0~20 cm土层土壤饱和导水率与土壤基本物理特性的相关程度排序相同,20~40 cm土层土壤饱和导水率与土壤基本物理特性的相关程度排序不同;土壤饱和导水率多尺度预测模型的预测精度较高,0~20 cm土层和20~40 cm土层拟合值的均方根误差分别为0.035 0和0.029 0;0~20 cm土层和20~40 cm土层土壤饱和导水率转换关系的计算精度较高,拟合值的均方根误差为0.037 5。     

Numerical modeling of the hydraulic signatures of horizontal and inclined faults

Pumping tests in discontinuous media are commonly interpreted using radial/non-radial flow analytical models, which do not account for a large range of multistage real signatures that combine both flow types. The term radial is employed in Barker??s sense, referring to the cross-flow area growth with a flow dimension (n) equal to 2, while multistage refers to a flow response of several time periods with successive n values. Synthetic domains are used to numerically investigate the flow behaviour of faulted media, specifically the effect of fault inclination on multistage signatures. Simulations are performed using single faults surrounded by a permeable matrix to characterize drawdown series and to identify hydrodynamic scenarios associated with some diagnostic multistage n sequences. The n sequence 2 ?C <?2 ?C 2 results from interactions between the matrix and a horizontal to weakly inclined fault, while 2 ?C 1.5 ?C 2 corresponds to a higher fault inclination and is a combination of the generalized horizontally and vertically faulted models. These results improve the interpretation of ambiguous drawdown series from long-term pumping tests. This provides a better hydraulic characterization of faulted environments, allowing assessment of fault and matrix properties such as fault inclination. The approach is successfully tested on real datasets.     

Numerical models of hydraulic fracturing and the interpretation of syntectonic veins

The hydraulic fracture is modelled as an ellipse in an infinite elastic medium with an internal fluid pressure and loaded under biaxial stresses at infinity. The available stress function for this model has been evaluated numerically, and the magnitudes of the stresses generated around the crack calculated for a variety of loading conditions and crack orientations. Fracture initiation is predicted from the Griffith maximum tensile stress criterion. The location of the maximum tensile stress around the crack is recorded and it is found that for many conditions of applied stresses and crack fluid pressure, the hydraulic shear fracture has a symmetrically developed maximum tensile stress and fracture initiation will occur by growth along the direction of the crack. It is also predicted that fracture initiation will occur when the ratio of fluid pressure to applied least principal stress is considerably less than one. The elastic strain energy fields around elliptical hydraulic flaws have been calculated, and in particular, the change in strain energy upon introduction of a small flaw, and the change in strain energy upon growth of this flaw, have been investigated. The results allow an evaluation of the second part of the Griffith criterion-that fracture growth is accompanied by a decrease in strain energy-for hydraulic fractures. Changes in strain energy with small increases in fluid pressure provide a physical basis for dilatancy hardening and fracture instability. Quasi-static growth from a flaw is modelled by calculating changes in strain energy for unit increases in half length. The distinction between fractures which show an increasing and a decreasing rate of change in strain energy with increasing length, and between fractures which may only extend spontaneously for short distances and those which may show extensive spontaneous growth on the basis of the rate of change of strain energy with length, is made. A gradual drop in crack fluid pressure once the threshold for fracture initiation has been passed may promote the extent to which spontaneous crack growth occurs.The formation of syntectonic veins, particularly in rocks being deformed under low grade metamorphic conditions, is often the most abundant evidence of natural hydraulic fracturing in rocks. Commonly observed geometric features of syntectonic veins-length, simple tapering, symmetric and asymmetric forking, branching, irregular zig-zag traces, en échelon patterns—are discussed primarily with reference to the strain energy model for growth established, and the geometric variation is interpreted in terms of variation in applied stress and fluid pressure conditions and the rate of change of stored strain energy with crack growth. In particular, terminal branching arises when the minimum stress changes from a symmetric to an asymmetric location at the tip of a growing shear fracture, and terminal forking results when there is an increase in the energy release rate during crack growth, and may be symmetric or asymmetric depending on the location of the minimum stress at the crack tip at the time of forking.     

土壤饱和导水率空间预测的不确定性分析

当土壤转换函数应用于土壤水力性质估计时,对于预测值的不确定性往往容易被忽视。为了有针对性地提出减少这种不确定性的方法和措施,提高土壤转换函数的实际应用能力,以两种现有的土壤转换函数(Vereecken和HYPRES模型)为例,将其应用于山东省平度市土壤饱和导水率的空间预测,并利用拉丁超立方抽样(LHS)方法对预测结果的不确定性进行了分析。结果表明,饱和导水率空间预测的不确定性主要来源于土壤基本性质的空间插值误差和土壤转换函数自身的预测误差。当Vereecken模型应用于饱和导水率空间预测时,预测结果的不确定性主要由土壤基本性质空间插值误差所决定,土壤转换函数预测误差的影响较小,而HYPRES模型则是受二者的双重影响。     

单孔岩样水压致裂的数值分析

水压致裂是改变岩体结构的一种天然行为和人为手段。采用F－RFPA2D软件,对水压致裂过程、裂缝扩展形态及注水孔形状和大小、应力条件和岩样强度等影响因素进行了研究。将开始出现声发射的水压称为微裂压力,将声发射急剧增多、裂缝非稳定扩展直至岩样破坏的水压称为破裂压力。岩样尺寸一定时,微裂压力和破裂压力随内孔面积增加而降低,方形孔岩样的微裂压力和破裂压力均小于同面积的圆形孔。微裂压力和破裂压力随围压或岩样强度增加而增加,且其差值随岩样强度增加而增加,理论破裂压力与模拟值趋势基本一致。方形孔的宏观裂纹起裂位置多在角点附近,而圆形孔比较随机。无围压时,宏观裂纹的延伸方向随机;有围压时,宏观裂纹扩展方向大致与主应力方向一致,且沿较大主应力方向的宏观裂纹扩展至岩样破坏,较小主应力方向宏观裂纹不完全发育。研究结果对水压致裂试验和工程实践有一定参考意义。     

人工掏挖坡面侵蚀微地貌演化及其水力学特性分析

人工掏挖是黄土高原地区夏闲地翻耕时广为应用的耕作方式,为探明其坡面侵蚀过程,利用三维激光扫描仪及ArcGIS软件,阐明其在人工模拟间歇降雨下坡面微地貌、侵蚀产沙及水力学参数演变规律。结果表明:①坡面侵蚀历经溅蚀-片蚀、断续细沟、连续细沟3个阶段;降雨截止时细沟平面密度、平均沟深、最大沟长和最大沟深分别增至初始的1.42倍、2.24倍、15.5倍和2.43倍。②地表糙度随降雨历时推移从1.706呈近似线性趋势减小至1.488;累积降雨量达80 mm之前,径流量、含沙量随地表糙度的减小增加缓慢,但之后随地表糙度减小增加剧烈。③降雨过程中流速呈现波动增加趋势,水流由层流快速过渡至紊流,但始终保持为缓流,水流阻力呈波动下降,且形态阻力一直居于主导地位。人工掏挖耕作坡面在雨强1.5 mm/min、降雨量80 mm以内可起到蓄水保土作用,但在连续强降雨下也更易引起细沟侵蚀。