Upscaling of elastic properties for large scale geomechanical simulations

Large scale geomechanical simulations are being increasingly used to model the compaction of stress dependent reservoirs, predict the long term integrity of under‐ground radioactive waste disposals, and analyse the viability of hot‐dry rock geothermal sites. These large scale simulations require the definition of homogenous mechanical properties for each geomechanical cell whereas the rock properties are expected to vary at a smaller scale. Therefore, this paper proposes a new methodology that makes possible to define the equivalent mechanical properties of the geomechanical cells using the fine scale information given in the geological model. This methodology is implemented on a synthetic reservoir case and two upscaling procedures providing the effective elastic properties of the Hooke's law are tested. The first upscaling procedure is an analytical method for perfectly stratified rock mass, whereas the second procedure computes lower and upper bounds of the equivalent properties with no assumption on the small scale heterogeneity distribution. Both procedures are applied to one geomechanical cell extracted from the reservoir structure. The results show that the analytical and numerical upscaling procedures provide accurate estimations of the effective parameters. Furthermore, a large scale simulation using the homogenized properties of each geomechanical cell calculated with the analytical method demonstrates that the overall behaviour of the reservoir structure is well reproduced for two different loading cases. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical prediction of blast‐induced stress wave from large‐scale underground explosion

This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.     

Optimization of grid-drilling using computer simulation

A computer simulation method has been developed to find efficient drilling grids for mineral deposits. A well-known ore deposit is used as a model to develop an efficient pattern for undiscovered ore bodies in the same area or in other prospects where similar geometry is suspected. The model for this study is the Austinville, Virginia deposit, a Mississippi Valley-type deposit composed of 17 ore bodies totaling 34 million short tons (30 million metric tons). The method employs a computer program that simulates drilling the model deposit with different patterns, including various levels of follow-up drilling. Follow-up holes are drilled in fences at one half the original spacing around holes in the grid that show ore-grade mineralization. Each pattern is drilled 100 times from random starting locations to provide a range of outcomes of drilling, including the best, worst, and most likely. For this study, patterns of 100 drill holes were composed of 10 fences spaced 1000–5000 feet (305–1524 m) apart, each with 10 holes spaced 200–1000 feet (61–305 m) apart. In all, 25 grids were used with zero to three levels of follow-up drilling. The 600/2000 grid, with drill holes spaced 600 feet (183 m) apart in fences spaced 2000 feet (610 m) apart, was compared with the 200/5000 grid because they represented contrasting outcomes. The 600/2000 grid penetrated many ore bodies consistently but with few multiple hits to individual ore bodies; whereas the 200/5000 grid inconsistently penetrated few ore bodies with many multiple hits. The 600/2000 grid was more efficient than the 200/5000 grid at hitting large ore bodies of 1,000,000 short tons or greater (900,000 metric tons or greater) and was made more effective by adding one cycle of follow-up drilling. The 600/2000 grid had a 97% chance of hitting one or more large ore bodies with at least one drill hole per ore body, and the 200/5000 grid had a 64% chance. Once hit, there was an 82% chance that the largest ore body would be penetrated by three or more holes when using the 600/2000 grid and an 88% chance using the 200/5000 grid.     

Is Precipitation the Dominant Controlling Factor of High Inorganic Nitrogen Content in the Changjiang River and Its Mouth？

The main reasons for the high content of inorganic N and its increase by several times in the Changjiang River and its mouth during the last 40 years were analysed in this work. The inorganic N in precipitation in the Changjiang River catchment mainly comes from gaseous loss of fertilizer N, N resulting from the increases of population and livestock, and from high temperature combustions of fossil fuels. N from precipitation is the first N source in the Changjiang River water and the only direct cause of high content of inorganic N in the Changjiang River and its mouth. The lost N in gaseous form and from agriculture non-point sources fertilizer comprised about 60% of annual consumption of fertilizer N in the Changjiang River catchment and were key factors controlling the high content of inorganic N in the Changjiang River mouth. The fate of the N in precipitation and other N sources in the Changjiang River catchment are also discussed in this paper.     

介质横向不均匀对震后地表位移场的影响

利用三维有限元模拟两层地球模型(包括弹性层和粘弹性层)介质的横向不均匀和断层错动的复杂性对震后地表位移场的影响，其数值结果包括弹性解和粘弹性解。结果表明，震后位移场的分布特征不仅与粘弹性层的流变特性有关，还与弹性层介质的剪切模量以及断层的同震错动有关。     

海拉尔盆地乌尔逊凹陷古流体动力场数值模拟

首先对地下水动力场进行数学描述,以水动力为主线,建立适合该地区的数学模型,选择正确的地质模型参数,主要采用数值计算和模拟等手段,在恢复出乌尔逊凹陷主要目的层在不同地质历史时期古水动力的基础上又恢复出古流体势的演化历史。整体上,在乌尔逊凹陷,各主要地质时期的油势演化特征为西部地区大于东部地区,乌南大于乌北,且低势区的发育有较好的继承性。最后,依据油气聚集理论,进行油气运移聚集系统的划分、油气运移方向的分析,从而得出有关研究区油气运移聚集条件的认识,主要从流体动力的角度考虑,指出有利聚集区,认为在乌尔逊凹陷苏4～苏18井区、苏132～苏17～铜2井区及巴1～巴2～巴4井区为有利油气运移聚集区。     

首都圈地区现今断层活动的数学模拟及地震危险性讨论

以首都圈地区现今活动断层上近２０年的位移测量资料为依据，用二维线弹性有限元对该区断层的活动特征进行了拟合，并结合有关资料讨论了该区的地震危险性。结果表明：１９７７－１９８６年间该区主压应力优势方位与华北较一致，约为Ｎ４５°－８０°Ｅ；１９８６－１９９０年其主压应力方向向北偏转，约为Ｎ５°－６０°Ｅ；１９９０年至现在其主压应力又向东偏转，角度大于第一阶段，约为Ｎ８０°－９５°Ｅ。该区的张家口－延庆一带近期有发生中强震的可能；丰镇－阳高－大同地区和凉城－古营盘地区应力较高，也应引起注意。     

Application of stress release models to historical earthquakes from North China

Univariate and multivariate stress release models are fitted to historical earthquake data from North China. It is shown that a better fit is obtained by treating separately the Eastern part of the region, including the North China Plain and Bohai Sea, and the Western part of the region, including the Ordos Plateau and its Eastern boundary. Further improvement is obtained by fitting the large events (M7.6) and smaller events in the Western region by different stress release models. The comparisons are made by computing the likelihoods of the fitted models and discounting the number of parameters used by Akaike's AIC criterion. The models are used to develop long-term risk scenarios for the East and West regions.     

RETRACTION: Numerical simulation of a wind-sand flow under different conditions of atmospheric stability

A theoretical model for wind‐sand flow is developed by considering the coupling between wind flow and sand particle motion, the latter subject to the Magnus effect, under different atmospheric stability conditions. Using this model, the characteristics of the wind‐sand flow are discussed in detail. The results show that the atmospheric stability and the Magnus effect both have a strong influence on wind profiles and on the trajectories of sand particles. This approach produces results with characteristics that differ from those previously reported; the latter only applying to atmospheric conditions of neutral stability. The saltating sand reaches a greater height under non‐neutral stability than under neutral stability, while the maximum horizontal distance is greater under unstable conditions and is smaller under stable conditions than under conditions of neutral stability.