Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.     

用双三次样条函数和GPS资料反演现今中国大陆构造形变场

将中国大陆现今构造变动视为一种连续的地壳变形,利用双三次样条函数模拟了近期GPS测定的大陆内部及周边地区412个测站速率,反演大陆地区自洽的构造变动速度场和应变率场.模拟结果显示:印度板块与欧亚板块的碰撞、挤压是构成中国大陆内部岩石层水平形变的主要驱动力.印度板块在东喜马拉雅构造结深深插入青藏高原,造成地壳大规模的缩短和抬升.青藏高原东南部的喜马拉雅带、拉萨和羌塘地块以及青藏高原东南边的川滇地区,内部构造活动强烈,其内部的构造变形包含地壳碎片的冲断、褶皱和侧向逃逸.大陆地壳(或岩石圈)的增厚,尤其是喜马拉雅山脉南北向的快速缩短和青藏高原东西向的缓慢拉张,大约吸收了印欧板块会聚量的85%,西藏中东地区东西向的拉张速率达到了(16±2.0)mm/a,且顺时针方向扭转明显.印度板块相对欧亚板块运动的欧拉极为(29.7°N, 19.3°E, 0.392°/Ma);华南地块相对于欧亚大陆向东(102°±7.4°)南的运动速率是(11±1.54)mm/a,华南块体相对欧亚板块运动的欧拉极为(62.25°N, 126.56°E, 0.141°/Ma);塔里木地块相对较稳定,其西部运动速度高于东部运动速度,作顺时针方向旋转.总体上讲,中国大陆运动方向为北偏东呈辐射状,从西部近南北方向的运动转向东部地区东南方向的运动,绕东喜马拉雅构造结有一顺时针方向的旋转.横穿喜马拉雅构造带及青藏内部的南北向压缩速率为(19±2.0)mm/a,横穿西天山构造带的南北向压缩平均速率为(13±1.5)mm/a,横穿东天山构造带的南北向压缩平均速率为(6.0±1.4)mm/a.阿尔金断裂带的左旋走滑速率为(6±1.2)mm/a.     

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.     

An operational multiscale system for hazards prediction, mapping, and response

By definition, a crisis is a situation that requires assistance to be managed. Hence, response to a crisis involves the merging of local and non-local emergency response personnel. In this situation, it is critical that each participant: (1) know the roles and responsibilities of each of the other participants; (2) know the capabilities of each of the participants; and (3) have a common basis for action. For many types of natural disasters, this entails having a common operational picture of the unfolding events, including detailed information on the weather, both current and forecasted, that may impact on either the emergency itself or on response activities. The Consequences Assessment Tool Set (CATS) is a comprehensive package of hazard prediction models and casualty and damage assessment tools that provides a linkage between a modeled or observed effect and the attendant consequences for populations, infrastructure, and resources, and, hence, provides the common operational picture for emergency response. The Operational Multiscale Environment model with Grid Adaptivity (OMEGA) is an atmospheric simulation system that links the latest methods in computational fluid dynamics and high-resolution gridding technologies with numerical weather prediction to provide specific weather analysis and forecast capability that can be merged into the geographic information system framework of CATS. This paper documents the problem of emergency response as an end-to-end system and presents the integrated CATS–OMEGA system as a prototype of such a system that has been used successfully in a number of different situations.     

Resources calculation of cobalt-rich crusts with the grid subdivision and integral method

On the basis of three geological models and several orebody boundaries,a method of grid subdivision and integral has been proposed to calculate and evaluate the resources of cobalt-rich crusts on the seamounts in the central Pacific Ocean.The formulas of this method are deduced and the interface of program module is designed.The method is carried out in the software "Auto mapping system of submarine topography and geomorphology MBChart".This method and program will possibly become a potential tool to calculate the resources of seamounts and determine the target diggings for China's next Five-year Plan.     

C网格嵌套技术及其在海洋波动传播模拟中的应用

采用粗细嵌套的ArakawaC网格模拟Klevin波和Rossby波沿赤道传播的过程，研究在粗细网格嵌套边界产生数值振荡和反射的原因及其消减方法，得到合理的并能用于复杂海洋模式的ArakawaC网格下的嵌套方案。数值实验结果表明：波形在粗细网格边界产生数值振荡的原因是波形在不同分辨率下的形态有差异，在粗网格下波形趋向于平坦化，而细网格下趋向于锐化。采用双向嵌套和粗细网格交界处加松弛的方法可以有效地消减数值振荡和反射。     

Parabolic Wave Propagation Modelling in Orthogonal Coordinate Systems

This paper presents a numerical model study of the propagation of water waves using the parabolic approximation of the mild-slope equation in the orthogonal coordinate system. Two types of coordinate systems are studied: (a) a general form of orthogonal coordinate system and (b) the conformal system, a special form of orthogonal coordinate system. Two typical examples, namely, expanded breakwaters and a circular channel, are studied to validate the model. First, the examples are studied by use of the general orthogonal coordinates. Then the same examples are computed by use of the confonnal system. The computational results show that the confonnal coordinate system generally gives better predictions than the general orthogonal system. A numerical technique for generating the conformal grid is combined with the numerical model to improve the practicability of the model. The comparison between the result from the numerical grid system and that from the analytical grid system shows that reliable computational results can be obtained by use of the numerical confonnal grid system.     

一种带自由面的三维船体贴体网络的数值生成

本文提出了一个生成自由面的三维船体贴体网格的数值方法。文中以静水面以上的船体型线作为描述非定常自由面的计算网格域，进一步完善了我们以往开发的无法考虑兴波与静水面以上的船体型线相互作用对船舶粘性流动影响的网格生成方法。静水面以上的初始自由面网格分别沿船面法线η及切线ζ方向的曲线长的指数进行分布，即可极方便地调节自由面在船面及静水面附近的网格间距，也保证了在数值离散三维Ｐｏｉｓｓｏｎ方程时具有相当好的收敛性。这种方法的特点是计算量相当小，尤其能极方便地控制网格在船面簿粘性层中及静水面附近的分布。以Ｓｅｒｉｅｓ６０船模作为计算算例，带自由面的网格分布是相当满意的     

ECOM模式在丁字湾的应用

应用河口、陆架和海洋模式 (ECOM模型 ) ,引入干湿网格法模拟潮滩涨落的改进 ,并建立丁字湾及近岸海域的三维变动边界潮流模型。该模型考虑了湾口拦门沙、湾内水道和人工围海等地形特点。计算结果与实测值比较符合良好 ,较好地刻画出丁字湾 M2 分潮潮流场的时空分布特点。     

Boussinesq-type modelling using an unstructured finite element technique

A model for solving the two-dimensional enhanced Boussinesq equations is presented. The model equations are discretised in space using an unstructured finite element technique. The standard Galerkin method with mixed interpolation is applied. The time discretisation is performed using an explicit three-step Taylor–Galerkin method. The model is extended to the surf and swash zone by inclusion of wave breaking and a moving boundary at the shoreline. Breaking is treated by an existing surface roller model, but a new procedure for the detection of the roller thickness is devised. The model is verified using four test cases and the results are compared with experimental data and results from an existing finite difference Boussinesq model.