Numerical estimation of REV and permeability tensor for fractured rock masses by composite element method

The Monte Carlo method is used to generate parent stochastic discrete fracture network, from which a series of fractured rock samples of different sizes and orientations are extracted. The fracture network combined with a regular grid forms composite element mesh of the fractured rock sample, in which each composite element is composed of sub‐elements incised by fracture segments. The composite element method (CEM) for the seepage is implemented to obtain the nodal hydraulic potential as well as the seepage flow rates through the fractured rock samples. The application of CEM enables a large quantity of stochastic tests for the fractured rock samples because the pre‐process is facilitated greatly. By changing the sizes and orientations of the samples, the analysis of the seepage characteristics is realized to evaluate the variation of the permeability components, the existence of the permeability tensor and the representative element volume. The feasibility and effectiveness are illustrated in a numerical example. Copyright © 2008 John Wiley & Sons, Ltd.     

Vertical variations of wave-induced radiation stress tensor

INTRODUcrIONThe concept of radiation stress was deve1oPed by tonguet--Higgins and Stewart (1964 ),who intreduced the definition of radiation stress as the excess mornentum due to the presence ofwaves, on the basis of time-averaged laws of Newtonian fluid mechanics and the assmption ofa unifOrm velocity distribution over depth. Subequently, the theory has been applied success-fully in the investigation of phenomena such as wave set-up and set--down (Bowen et al.,l968), longshore currents …     

水成沉积与风成沉积及古土壤的磁组构特征

晚第四纪水成沉积、风成沉积、古土壤等的磁组构研究表明,不同成因的沉积物具有明显不同的磁组构特征。一般地,水成沉积的体积磁化率Ｋ受物源和水动力因素控制,Ｐ和Ｆ均＞１．０２,ｑ＜０．５;风成沉积的Ｋ受物源控制,内陆黄土的Ｋ明显高于沿海黄土及海岸风沙Ｋ,Ｐ和Ｆ均＜１．０２,ｑ的平均值＞０．５５;古土壤的Ｋ受母质沉积类型和成土环境的氧化还原条件控制,内陆古土壤的Ｋ偏高,沿海古土壤的Ｋ偏低,内陆古土壤的Ｋ远高于沿海古土壤的Ｋ,Ｐ和Ｆ均＜１．０２,ｑ的平均值在０．５左右。运用典型沉积的磁组构特征对比分析未知沉积物,可以得到未知沉积物的物源、沉积过程和环境信息。     

哀牢山-红河断裂带中段应变分析

对哀牢山—红河断裂带中段和平—水塘剖面、墨江—元江剖面和其它地段的岩石应变及磁组构进行了分析，表明应变强度的校正磁各向异性度PJ从断裂带向西至三叠系明显降低，变形强度向西迅速减弱。磁化率椭球体主轴展布反映出剪切带内、外变形方式的改变。剪切带内，以水平走滑运动为主；向西则以水平缩短为主。在应变分析中，对断裂带内的S—C组构两组面理夹角、杏仁体和石榴石应变标志体进行了测量，结果显示，剪切带内的磁组构和岩石组构间关系较差。此外，还对这两条剖面中的三叠系进行了应变和磁组构分析，找出了它们与磁化率椭球体对应轴率间的相关性。     

俯冲带榴辉岩的变形作用及其对俯冲-折返过程的意义

榴辉岩是大洋和大陆俯冲带的重要岩石类型,在研究俯冲带的形成过程、热结构、壳幔相互作用等方面有重要意义。通过对天然和实验样品中石榴子石、绿辉石等矿物的变形特征、变形机制、变形的影响因素等的综合分析,系统总结了高压变质带中榴辉岩矿物显微和超微变形研究的进展,探讨了榴辉岩的变形特征在恢复俯冲与折返过程研究中的意义及一些尚待解决的一些问题。     

Rupture process of large earthquakes in the northern Mexico subduction zone

The Cocos plate subducts beneath North America at the Mexico trench. The northernmost segment of this trench, between the Orozco and Rivera fracture zones, has ruptured in a sequence of five large earthquakes from 1973 to 1985; the Jan. 30, 1973 Colima event (M _s 7.5) at the northern end of the segment near Rivera fracture zone; the Mar. 14, 1979 Petatlan event (M _s 7.6) at the southern end of the segment on the Orozco fracture zone; the Oct. 25, 1981 Playa Azul event (M _s 7.3) in the middle of the Michoacan gap; the Sept. 19, 1985 Michoacan mainshock (M _s 8.1); and the Sept. 21, 1985 Michoacan aftershock (M _s 7.6) that reruptured part of the Petatlan zone. Body wave inversion for the rupture process of these earthquakes finds the best: earthquake depth; focal mechanism; overall source time function; and seismic moment, for each earthquake. In addition, we have determined spatial concentrations of seismic moment release for the Colima earthquake, and the Michoacan mainshock and aftershock. These spatial concentrations of slip are interpreted as asperities; and the resultant asperity distribution for Mexico is compared to other subduction zones. The body wave inversion technique also determines theMoment Tensor Rate Functions; but there is no evidence for statistically significant changes in the moment tensor during rupture for any of the five earthquakes. An appendix describes theMoment Tensor Rate Functions methodology in detail.The systematic bias between global and regional determinations of epicentral locations in Mexico must be resolved to enable plotting of asperities with aftershocks and geographic features. We have spatially shifted all of our results to regional determinations of epicenters. The best point source depths for the five earthquakes are all above 30 km, consistent with the idea that the down-dip edge of the seismogenic plate interface in Mexico is shallow compared to other subduction zones. Consideration of uncertainties in the focal mechanisms allows us to state that all five earthquakes occurred on fault planes with the same strike (N65°W to N70°W) and dip (15±3°), except for the smaller Playa Azul event at the down-dip edge which has a steeper dip angle of 20 to 25°. However, the Petatlan earthquake does prefer a fault plane that is rotated to a more east-west orientation—one explanation may be that this earthquake is located near the crest of the subducting Orozco fracture zone. The slip vectors of all five earthquakes are similar and generally consistent with the NUVEL-predicted Cocos-North America convergence direction of N33°E for this segment. The most important deviation is the more northerly slip direction for the Petatlan earthquake. Also, the slip vectors from the Harvard CMT solutions for large and small events in this segment prefer an overall convergence direction of about N20°E to N25°E.All five earthquakes share a common feature in the rupture process: each earthquake has a small initial precursory arrival followed by a large pulse of moment release with a distinct onset. The delay time varies from 4 s for the Playa Azul event to 8 s for the Colima event. While there is some evidence of spatial concentration of moment release for each event, our overall asperity distribution for the northern Mexico segment consists of one clear asperity, in the epicentral region of the 1973 Colima earthquake, and then a scattering of diffuse and overlapping regions of high moment release for the remainder of the segment. This character is directly displayed in the overlapping of rupture zones between the 1979 Petatlan event and the 1985 Michoacan aftershock. This character of the asperity distribution is in contrast to the widely spaced distinct asperities in the northern Japan-Kuriles Islands subduction zone, but is somewhat similar to the asperity distributions found in the central Peru and Santa Cruz Islands subduction zones. Subduction of the Orozco fracture zone may strongly affect the seismogenic character as the overlapping rupture zones are located on the crest of the subducted fracture zone. There is also a distinct change in the physiography of the upper plate that coincides with the subducting fracture zone, and the Guerrero seismic gap to the south of the Petatlan earthquake is in the wake of the Orozco fracture zone. At the northern end, the Rivera fracture zone in the subducting plate and the Colima graben in the upper plate coincide with the northernmost extent of the Colima rupture zone.     

离散裂隙渗流方法与裂隙化渗透介质建模

流体渗流模拟的连续介质方法通常适用于多孔地质体，并不一定适用于裂隙岩体，由于裂隙分布及其特征与孔隙差异较大。若流体渗流主要受裂隙的控制，对于一定尺寸的裂隙岩体，多孔介质假设则较难刻划裂隙岩体的渗流特征。离散裂隙渗流方法不但可直接用于模拟裂隙岩体非均质性和各向异性等渗流特征，而且可用其确定所研究的裂隙岩体典型单元体及其水力传导（渗透）张量大小。主要讨论了以下问题：（1）饱和裂隙介质中一般的离散流体渗流模拟；（2）裂隙岩体中的REV（典型单元体）及其水力传导（渗透）张量的确定；（3）利用离散裂隙网络流体渗流模型研究裂隙方向几何参数对水力传导系数和REV的影响；（4）在二维和三维离散裂隙流体渗流模型中对区域大裂隙和局部小裂隙的处理方法。调查结果显示离散裂隙流体渗流数学模型可用来评价不同尺度上的裂隙岩体的水力特征，以及裂隙方向对裂隙化岩体的水力特征有着不可忽视的影响。同时，局部小裂隙、区域大裂隙应当区别对待，以便据其所起的作用及水力特征，建立裂隙化岩体相应的流体渗流模型。