A coupling of mixed and discontinuous Galerkin finite-element methods for poroelasticity

In this paper, we formulate a finite-element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. We approximate the flow variables by a mixed finite-element space and the displacement by a family of discontinuous Galerkin methods. Theoretical convergence error estimates are derived and, in particular, are shown to be independent of the constrained specific storage coefficient, c _o . This suggests that our proposed algorithm is a potentially effective way to combat locking, or the nonphysical pressure oscillations, which sometimes arise in numerical algorithms for poroelasticity.     

研究“源兆”的若干问题及途径分析

从位场理论的观点讨论了场和源的概念,建议将“场兆”改称为构造前兆,将“源兆”称为震源前兆。文中还分析了震源体与各种前兆观测方法探测的深度,并根据地磁与重力观测方法不同的物理基础和地质基础,指出震磁效应对于研究震源前兆的特殊意义,文中还介绍了１９７６唐山地震观测到震源前兆的部分实际地磁资料。     

The Epi－continental arc of Southeast China and relevant earthquakes

ＴｈｅＥｐｉ┐ｃｏｎｔｉｎｅｎｔａｌａｒｃｏｆＳｏｕｔｈｅａｓｔＣｈｉｎａａｎｄｒｅｌｅｖａｎｔｅａｒｔｈｑｕａｋｅｓＪＩＡ－ＷＥＩＸＵ（徐嘉炜）ＤｅｐａｒｔｍｅｎｔｏｆＲｅｓｏｕｒｃｅａｎｄＥｎｖｉｒｏｎｍｅｎｔａｌＳｃｉｅｎｃｅｓ,ＨｅｆｅｉＵ...     

A quantitative analysis of the prototype [WCL] star CPD-56° 8032

We present a detailed, quantitative study of the standard [WC10] Wolf-Rayet central star CPD-56^o 8032 based on new high resolution AAT UCLES observations and the Hillier (1990) WR standard model. Our analysis of CPD-56^o 8032 gives the wind properties (T _*=34500K, lg (L/L )=3.8, lg (M/M a^–1)=–5.4,v =225 km s^–1) and chemistry (C/He=0.5, O/He=0.1, by number), the latter suggesting an intimate relationship with the Ovi PN central stars and the PG 1159-035 objects. A comparison between the wind properties of CPD-56^o 8032 and Sk-66^o 40 (WN 10) indicates that low excitation, low wind velocity WR winds are common to both low mass PN central stars (WC sequence) and high mass post-LBV's (WN sequence).     

Blind prediction of the site effects at Ashigara Valley,Japan, and its comparison with reality

Weak and strong ground motions were numerically predicted for three stations of the Ashigara Valley test site. The prediction was based on the records from a rock-outcrop station, one weak-motion record from a surface-sediments station, and the standard geotechnical model. The data were provided by the Japanese Working Group on the Effects of Surface Geology as a part of an international experiment. The finite-difference method for SH waves in a 2-D linear viscoelastic medium (a causalQ model) was employed.Comparison with the real records shows that at two stations the predictions fit better than at the third one. Strangely, the two better predictions were for stations situated at larger distances from the reference rock station (one station was on the surface, the other in a borehole). The strong ground motion (the peak acceleration of about 200 cm s^–2) was not predicted qualitatively worse than the weak motion (8 cm s^–2). A less sophisticated second prediction (not submitted during the experiment), in which we did not attempt to fit the available weak-motion record at the sedimentary station, agrees with the reality significantly better.     

Automatic equipment to monitor groundwater radon

Prototype instrumentation, able to automatically measure groundwater radon content variations, is presented. The equipment is made of stainless steel and has spherical valves with automatic and pneumatic control. The deemanation of the gases from the water is obtained by evacuating a suitable expansion chamber. The instrumentation can make discrete sampling ranging from 1 per hour to 1 per 99 hours. The equipment was tested in the laboratory: the efficiency was measured by means of a²⁶⁶Ra solution. A mean value of (0.65±0.07) count/s/Bq was obtained. A calibration test was carried out by comparing countings from the automatic equipment with those obtained by the standard laboratory cell. Results of an operational check over a period of approximately one year indicate that variations in radon at the calibration site are attributable more to meteorological than to tectonic causes.     

1983年菏泽5.9级地震前后沂水泉氡的变化

沂水泉位于山东省沂水一汤头断裂带上。多年来的观测分析表明,影响沂水泉氡值变化的最大干扰是降雨。由于这种干扰存在着“记忆”滞后影响,所以用一般的相关分析难以将这种影响排除干净。本文对最小二乘法为建立模型的判据,用动态灰箱分析法拟合并推估观测系统的变化。在分析中用了四个状态变量：长趋势成分Ｌ（ｔ）、年周期变化成分Ｐ（ｔ）、降雨影响部分Ｒ（ｔ）和随机变化部分Ｓ（ｔ）。经过分析表明,沂水泉氡的长趋势部分１     

Possible creep-related tilt precursors obtained in the Central Apennines (Italy) and in the Southern Caucasus (Georgia)

Daily averaged tilt component data from two sites of the Central Apennines (Italy) and of the Southern Caucasus (Georgia), respectively, revealed intermediate-term tilts as possible precursors to earthquakes (M=3.0÷4.7) which occurred in the above-mentioned seismic areas within a distance of 50 km from the sites. A good temporal correlation as well as a fair spatial correspondence between these residual tilts (with amplitude and duration of some microradians and months, respectively) and main shocks were pointed out, by removing both secular trends and seasonal thermoelastic effects from the raw tilts. An attempt was made to justify the above-mentioned results, based on the assumption that the observed intermediate-term preseismic tilts are the manifestation of aseismic creep episodes of comparable duration in the fault materials of thrust faults close to the tilt sites. The mechanism refers to a strain field slowly propagating from the preparation focal area to the tilt site, through crustal blocks separated by weak transition zones. This propagation is thought to be the cause of the local aseismic fault slip recorded by the tiltmeters. Previously, both discrete structures and strain propagation effects were revealed in the Central Apennines and are thought also to exist in the Southern Caucasus. As in the past, the rheological properties of fault materials are revealed as viscoelastic ones. In fact, creep equations obtained by applying several viscoelastic models on our data, proved to fit quite well some of the observed tilt precursors, producing viscosity and rigidity values very similar to those reported in literature.Professor Petr Viktorovich Manjgaladze died during the writing of this paper     

The circumstellar shells of carbon miras

The relation between mass loss rate and pulsation period in carbon Miras is discussed. The dust mass loss rate is very low (about 2 × 10^–10 M/yr) up to aboutP = 380 days, where there is a sudden increase. ForP > 400 days there is a linear relation between logM andP. The change in the mass loss rate near 380 days may be related to radiation pressure on dust becoming effective in driving the outflow.     

A combined study of macroseismic data and focal mechanisms applied to the West-Bohemian earthquake,Czechoslovakia

The main shock of the West-Bohemian earthquake swarm, Czechoslovakia, (magnitudem=4.5, depthh=10 km) exhibits an irregular areal distribution of macroseismic intensities 6° to 7° MSK-64. Four lobes of the 6° isoseismal are found and the maximum observed intensity is located at a distance of 8 km from the instrumentally determined epicentre. This distribution can be explained by the energy flux of the directS wave generated by a circular source, the hypocentral location and focal mechanism of which are taken from independent instrumental studies. The theoretical intensity, which is assumed to be logarithmically proportional to the integrated squared ground-motion velocity (i.e.,I=const+log v ² (t)dt), fits the observed intensity with an overall root-mean-square error less than 0.5°. It is important that the present intensity data can also be equally well explained by the isotropic source. The fit was attained by means of a horizontally layered model though large fault zones and an extended sedimentary basin suggest a significant lateral heterogeneity of the epicentral region. The results encourage a broader application of the simple modelling technique used.