中国大陆科学钻探主孔声发射法现今地应力状态的确定

对中国大陆科学钻探主孔的岩心进行了声发射测量,确定了301～1531m深度的最大主应力。并与钻孔崩落法(深度1269～1655m)测量结果进行了对比,结果表明,声发射测量所得测值基本上落在钻孔崩落法测值的趋势线上,两种方法所得结果一致,说明测量结果可信。测量结果表明科学钻探主孔地应力大小随深度增加,在浅部301m最大主应力为13.4MPa,在深部1655m为55.2MPa。随深度的增加率为0.0279MPa/m。最大主应力方向为N54°±3.3°E,且方向不随深度变化。     

山东省临邑县临盘镇地裂冒喷水油灾害初步分析

1998年7月临邑县临盘镇西十二里村发生地裂冒喷水油灾害：对此次灾害研究认为，灾害发生前地应力已有明显改变，这种改变应是自然因素造成的，油田开发难以直接引发断层的活动。因此，地裂冒喷水、油灾害是因活动性断昙——临邑断层在局部薄弱部位活动而引起的，断层活动使地下流体顺断层涌出地表，是自然地质灾害。灾害的诱发因素是多方面的，其中地面不均衡沉降耦合使得地下应力布不均习，导致在断层的薄弱部位产生应力集中，并引发灾害发生。灾害发生已使得地下不均衡应力得以释放，近期内再次发生灾害的可能性不大。     

邯峰矿区古构造应力场与构造演化

采用共轭剪节理应力反演方法,恢复了邯郸-峰峰矿区晚古生代以来的3期古构造应力场,进而探讨了煤田构造的演化历史,将其分为4大阶段：①中生代早期近NS向挤压,煤系后期改造初动期;②中生代晚期SE-NW向挤压,奠定煤田构造格架的基础;③中生代末至古近纪NW-SE向拉张,煤田构造格架定型;④新近纪以来近东西向拉张,煤田构造的现代活动.     

三峡深水高土石围堰工程的渗流研究

为了对三峡工程二期高土石围堰防渗设施的布置方案及其阻渗效果进行比较,采用有限元法对二期高土石围堰在不利的运行工况下的渗流场进行数值分析.比较了双排混凝土防渗墙方案、单排塑性混凝土低防渗墙加土工膜斜墙和单排厚塑性混凝土防渗墙等3种方案,分别采用恒定与非恒定模型计算,对立面二维和三维绕渗及防渗墙局部开裂等不利工况分别进行了数值模拟.结果表明:三方案均可有效抑制渗流场,双排混凝土防渗墙的防渗效果最好;墙体的局部开裂仅对局部区域的流场有影响.非恒定数值分析表明,堰体、基础不均匀沙石料及基坑抽水速度对渗流场影响极大,为保证堰体稳定,应限制基坑水位降落速度小于2 m/d.     

Outlier detection algorithms and their performance in GOCE gravity field processing

The satellite missions CHAMP, GRACE, and GOCE mark the beginning of a new era in gravity field determination and modeling. They provide unique models of the global stationary gravity field and its variation in time. Due to inevitable measurement errors, sophisticated pre-processing steps have to be applied before further use of the satellite measurements. In the framework of the GOCE mission, this includes outlier detection, absolute calibration and validation of the SGG (satellite gravity gradiometry) measurements, and removal of temporal effects. In general, outliers are defined as observations that appear to be inconsistent with the remainder of the data set. One goal is to evaluate the effect of additive, innovative and bulk outliers on the estimates of the spherical harmonic coefficients. It can be shown that even a small number of undetected outliers (<0.2 of all data points) can have an adverse effect on the coefficient estimates. Consequently, concepts for the identification and removal of outliers have to be developed. Novel outlier detection algorithms are derived and statistical methods are presented that may be used for this purpose. The methods aim at high outlier identification rates as well as small failure rates. A combined algorithm, based on wavelets and a statistical method, shows best performance with an identification rate of about 99%. To further reduce the influence of undetected outliers, an outlier detection algorithm is implemented inside the gravity field solver (the Quick-Look Gravity Field Analysis tool was used). This results in spherical harmonic coefficient estimates that are of similar quality to those obtained without outliers in the input data.     

A parametric study on the impact of satellite attitude errors on GOCE gravity field recovery

In the recent design of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) satellite mission, the gravity gradients are defined in the gradiometer reference frame (GRF), which deviates from the actual flight direction (local orbit reference frame, LORF) by up to 3–4°. The main objective of this paper is to investigate the effect of uncertainties in the knowledge of the gradiometer orientation due to attitude reconstitution errors on the gravity field solution. In the framework of several numerical simulations, which are based on a realistic mission configuration, different scenarios are investigated, to provide the accuracy requirements of the orientation information. It turns out that orientation errors have to be seriously considered, because they may represent a significant error component of the gravity field solution. While in a realistic mission scenario (colored gradiometer noise) the gravity field solutions are quite insensitive to small orientation biases, random noise applied to the attitude information can have a considerable impact on the accuracy of the resolved gravity field models.     

GGM02 – An improved Earth gravity field model from GRACE

A new generation of Earth gravity field models called GGM02 are derived using approximately 14 months of data spanning from April 2002 to December 2003 from the Gravity Recovery And Climate Experiment (GRACE). Relative to the preceding generation, GGM01, there have been improvements to the data products, the gravity estimation methods and the background models. Based on the calibrated covariances, GGM02 (both the GRACE-only model GGM02S and the combination model GGM02C) represents an improvement greater than a factor of two over the previous GGM01 models. Error estimates indicate a cumulative error less than 1 cm geoid height to spherical harmonic degree 70, which can be said to have met the GRACE minimum mission goals. Electronic Supplementary Material Supplementary material is available in the online version of this article at     

Global height datum unification: a new approach in gravity potential space

The problem of “global height datum unification” is solved in the gravity potential space based on: (1) high-resolution local gravity field modeling, (2) geocentric coordinates of the reference benchmark, and (3) a known value of the geoid’s potential. The high-resolution local gravity field model is derived based on a solution of the fixed-free two-boundary-value problem of the Earth’s gravity field using (a) potential difference values (from precise leveling), (b) modulus of the gravity vector (from gravimetry), (c) astronomical longitude and latitude (from geodetic astronomy and/or combination of (GNSS) Global Navigation Satellite System observations with total station measurements), (d) and satellite altimetry. Knowing the height of the reference benchmark in the national height system and its geocentric GNSS coordinates, and using the derived high-resolution local gravity field model, the gravity potential value of the zero point of the height system is computed. The difference between the derived gravity potential value of the zero point of the height system and the geoid’s potential value is computed. This potential difference gives the offset of the zero point of the height system from geoid in the “potential space”, which is transferred into “geometry space” using the transformation formula derived in this paper. The method was applied to the computation of the offset of the zero point of the Iranian height datum from the geoid’s potential value W ₀=62636855.8 m²/s². According to the geometry space computations, the height datum of Iran is 0.09 m below the geoid.     

转光膜在棉田中的应用试验与示范

1998-1999年对转光膜在和田市定点试验和大田示范，转光膜具有比普通膜更明显的增温及提高棉花苗期抗低温能力、促进棉花生长发育的作用。棉田使用转光膜后，可提高棉花品质，增产增收效果显著。     

Laboratory Studies Of Wind Stress Over Surface Waves

Simultaneous laboratory observations of wind speed, wind stress, and surfacewind-wave spectra are made under a variety of wind forcing patterns using cleanwater as well as water containing an artificial surfactant. Under typical experimentalconditions, more than half of the total stress is supported by the wave-induced stressrather than by the surface viscous stress. When the surfactant reduces the shortwind-wave spectra, the wind stress also decreases by as much as 20–30% at agiven wind speed. When the wind forcing is modulated in time, the wind stresstends to be higher under decreasing wind than under increasing wind at a givenwind speed, mainly because the response of short wind-wave spectra to varyingwind forcing is delayed in time. These examples clearly demonstrate that therelationship between the wind speed and the wind stress can be significantlymodified if the surface wave field is not in equilibrium with the wind forcing.Next, we examine whether the wind stress is estimated accurately if the wave-inducedstress by all surface wave components is explicitly evaluated by linear superpositionand is added to the surface viscous stress. It is assumed that the surface viscous stressis uniquely related to the wind speed, and that the wind input rate is determined by thelocal, reduced turbulent stress rather than the total stress. Our wind stress estimatesincluding the wave contributions agree well with observed wind stress values, evenif the surface wave field is away from its equilibrium with the wind in the presenceof surface films and/or under time-transient wind forcing. These observations stronglysuggest that the wind stress is accurately evaluated as a sum of the wave-induced stressand the surface viscous stress. At very high winds, our stress estimates tend to be lowerthan the observations. We suspect that this is because of the enhancement of wind stressover very steep (or breaking) short wind-waves.