抽水引起的含水层水平应变──地裂缝活动新机理

基于Darcy -Gersevanov广义渗流关系式 ,本文首先导出了单井定流量抽水引起的泰斯承压含水层水平运动和应变解析表达式 ,并据此提出了抽水诱发地裂缝活动的含水层水平应变新机理 ,其主要观点为 :受井栅的约束阻挡作用 ,抽水活动将在含水层内部形成近井径向挤压区和远井径向拉张区两个应变区域 ,随着抽水时间的增长 ,近井径向挤压区将逐渐向外扩展 ;当地裂缝及其下伏断层或裂隙位于抽水活动引起的径向拉张应变区域时 ,地裂缝将出现加速活动 ;反之 ,当地裂缝处于抽水活动引起的径向挤压应变区域时 ,地裂缝活动反而受到抑制。上述地裂缝活动机理通过大同机车工厂现场抽水实验得到了初步证实。     

Pseudodynamic tests on rubber base isolators with numerical substructuring of the superstructure and strain‐rate effect compensation

A pseudodynamic testing procedure has been applied by which the seismic response of a base‐isolated building is obtained by using as specimen the isolators, while the superstructure is numerically simulated. The procedure also takes advantage of the continuous pseudodynamic testing capabilities of the ELSA laboratory, which increase the accuracy of the results and reduce the strain‐rate effect of the rubber bearings. A simple proportional correction of the measured forces compensates the remaining strain‐rate effect due to the unrealistic speed of the test. The correction factor is obtained by means of a characterizing test on the specific rubber isolators. The developed method has been successfully applied to the prediction of the seismic response of a base‐isolated four‐storey building submitted to several specified accelerograms. The results for those earthquakes as well as the effects of some changes of the parameters of the system are discussed. Copyright © 2002 John Wiley & Sons, Ltd.     

遗传算法在边坡数值计算中的应用

改进了进化方向的遗传算法与有限元数值法，结合并研制了相应的软件。应用该软件对多类型岩土边坡进行弹性模量、内聚力、内摩擦角等参数反演分析，显示误差很小，收敛速度也很快，这说明改进进化方向遗传算法这种新型的优化算法在多类型岩土参数优化估计中具有独特的优势。     

An experimental study of grain scale melt segregation mechanisms in two common crustal rock types

Creation of pathways for melt to migrate from its source is the necessary first step for transport of magma to the upper crust. To test the role of different dehydration‐melting reactions in the development of permeability during partial melting and deformation in the crust, we experimentally deformed two common crustal rock types. A muscovite‐biotite metapelite and a biotite gneiss were deformed at conditions below, at and above their fluid‐absent solidus. For the metapelite, temperatures ranged between 650 and 800 °C at P_c=700 MPa to investigate the muscovite‐dehydration melting reaction. For the biotite gneiss, temperatures ranged between 850 and 950 °C at P_c=1000 MPa to explore biotite dehydration‐melting under lower crustal conditions. Deformation for both sets of experiments was performed at the same strain rate (ε.) 1.37×10^?5 s^?1. In the presence of deformation, the positive ΔV and associated high dilational strain of the muscovite dehydration‐melting reaction produces an increase in melt pore pressure with partial melting of the metapelite. In contrast, the biotite dehydration‐melting reaction is not associated with a large dilational strain and during deformation and partial melting of the biotite gneiss melt pore pressure builds more gradually. Due to the different rates in pore pressure increase, melt‐enhanced deformation microstructures reflect the different dehydration melting reactions themselves. Permeability development in the two rocks differs because grain boundaries control melt distribution to a greater extent in the gneiss. Muscovite‐dehydration melting may develop melt pathways at low melt fractions due to a larger volume of melt, in comparison with biotite‐dehydration melting, generated at the solidus. This may be a viable physical mechanism in which rapid melt segregation from a metapelitic source rock can occur. Alternatively, the results from the gneiss experiments suggest continual draining of biotite‐derived magma from the lower crust with melt migration paths controlled by structural anisotropies in the protolith.     

Numerical Analysis of Space–Time Finite Element Formulations for Miscible Displacements

A finite element formulation is proposed to approximate a nonlinear system of partial differential equations, composed by an elliptic subsystem for the pressure–velocity and a transport equation (convection–diffusion) for the concentration, which models the incompressible miscible displacement of one fluid by another in a rigid porous media. The pressure is approximated by the classical Galerkin method and the velocity is calculated by a post-processing technique. Then, the concentration is obtained by a Galerkin/least-squares space–time (GLS/ST) finite element method. A numerical analysis is developed for the concentration approximation. Then, stability, convergence and numerical results are presented confirming the a priori error estimates.     

Fluid permeability of sedimentary rocks in a complete stress–strain process

The permeability of sedimentary rocks during triaxial compression tests was investigated to relate it to the complete strain–stress process. It was found that the permeability was not constant, but varied with the stress and strain states in the rocks. Prior to the peak strength, the permeability decreases with increasing load. A dramatic increase in permeability occurs during the strain softening period. In the present study, in situ measurements of fluid flow and pressure in floor strata was carried out in a double longwall mining face in the Yangzhuang colliery. These measurements show that both the strata pressure and the position with respect to the mining face influence the hydrogeologic properties. The permeability increased in the floor strata behind the mining face because those mining induced fractures opened as the strata pressure decreased. To better understand this change in hydraulic behavior around the mining faces, 3-D numerical modeling was carried out. The model provides the general picture of the stress distribution and failure zone both in the floor and roof strata. The field and model results demonstrate the importance of changes in the stress and strain states on the hydrogeology of a site.     

地震前兆时空非均匀性指标Cv值的实验检验

在双向加载条件下，对两种构造和介质不同的岩石标本进行变形实验，测量了标本表面不同部位的应变．在此基础上，利用王晓青和陈学忠等提出的描述地震前兆非均匀分布的参量————Cv值，分析了标本变形失稳过程中应变异常的分布，以期从实验角度检验Cv值方法并探讨其物理意义．研究表明，Cv值的变化与岩石变形特征的变化有关，是描述前兆分布非均匀性的一种有效指标；Cv值在失稳发生前先异常上升，然后下降，并在Cv值恢复过程中或恢复后伴随着失稳事件的发生．因此，Cv值是一种有效的预报指标．      

由GPS观测结果推导中国大陆现今水平应变场

以中国大陆及周边近400个GPS测站的水平运动速率为基础，给出了现今地壳水平应变场结果表明:①中国大陆水平应变为西强东弱，剪应变数值大于正应变数值(绝对值)，应变量级一般为10-8/a，局部区域达到10-7/a，但应变分布不均匀；②南北向应变最突出的部位为中国西南部西段的喜马拉雅条带、西北部的36N~42N段及柴达木断块的北缘；③东西向应变西边缘变化最大．此外，由西向东还具有正负交替的变化特征；④REN(东-北向剪切应变)与Rmax(最大剪切应变)数值较大的区域分别是喜马拉雅条带、西北部的36N～42N段、柴达木断块的西部、川滇菱形块体，以及阿拉善、祁连及塔里木断块的交界区；⑤青藏块体周边以面收缩为主，内部则以面膨胀为主．其以北的地区以面收缩为主．西界数值最大，东部数值最小(除燕山构造带外)；⑥西部西区主压应变为南北向，主张应变为东西向．西部东缘区主压应变为近东西向，主张应变为近南北向．川滇菱形块体主应变的方向发生了很大的变化，北部地区为东西压南北张，南部地区则恰好相反；⑦中国大陆的应变模式可能是断块模式与连续模式的组合．此外，小尺度优势应变可能是剪切应变．造成上述结果除与印度板块的碰撞及边界耦合有关外，还与深部物质的活动及地壳介质的物性有密切的关系．必须指出，由于GPS测站在空间上分布的不均匀性，那么，由此而来的应变场，其应变尺度也不一样．      

用应变积累释放模型研究中国大陆地块分区地震活动

依据活动地块假说及活动边界的划分，对中国大陆地震活动进行分区．在此基础上，应用应变积累释放模型，通过地震应变积累释放的计算，研究了中国大陆各活动地块分区的地震活动性，初步探讨了各地块的地震活动阶段划分．此外，还对应变积累释放模型的应用前提条件进行了必要的讨论，并对模型结果中可能出现的问题给予一定分析．      

Improved treatments for evaluating horizontal magnetic components through the 3-D FDM in E -polarization induction problems

To improve the accuracy of the numerical evaluation through the 3-D finite difference method, the surface boundary conditions are added to modify the old program. The author has tested the new program by making calculations for the model constructed by Wanamaker, et al (1984). The comparison between the numerical results obtained from this paper and those by Wannamaker, et al (1984) indicates that a pronounced improvement is realized in the evaluation of the horizontal magnetic components. Moreover, better calculations for the vertical magnetic components are also obtainable by using the new program.