长江日流量混沌变化特性研究——Ⅰ相空间嵌入滞时的确定

为了以新途径探索长江水量变化规律并在此基础上做出预测,全面而系统地研究了宜昌站日流量混沌变化特性并建立了新预测模型.重点论述相空间嵌入滞时的确定.分析了传统的自相关函数法和真实矢量场法确定嵌入滞时的优劣,并提出以广义相关函数合理确定嵌入滞时的方法.研究结果表明广义相关函数法是一个合理又可操作的途径.     

汪清7.2级深震地下流体动态异常探析

2002年6月29日汪清7．2级深震前，吉林省地下流体观测网部分井孔观测到一组从趋势到短期前兆异常。本文阐述这些异常出现的特征和时空演化过程。认为异常出现的时间与震中距密切相关，其空间展布具有迁移性，中期异常由外围向震中收缩，其速率很小，进入短期阶段，岩石处于失稳、临失稳状态，水位异常由震中向外围发散，速率迅速加大。同时发现了异常盲区即深震的震中区，它对预测深震中有非常实际的意义，同时对研究板块运动、地球动力学过程也有一定的参考价值。     

青海省3次强震前地下流体短临异常的对比分析

通过对发生于1999～2000年间库玛中东段地区3次5级以上地震前地下流体短临异常的分析对比，得出了该地区地下流体短临异常的时空强演化特征。     

Fluid mixing as the mechanism of formation of the Dajing Cu-Sn-Ag-Pb-Zn ore deposit,Inner Mongolia ——Fluid inclusion and stable isotope evidence

Since the 1990s, interest in the magmatic fluids and their relation to mineralization has been re-aroused[1—6]. Studies on stable isotopes of low-sulfidation deposits commonly show the predominance of meteoric water[7]. Paradoxically, the evidence for me…     

Sediment deposition from flow at low gradients into a buffer strip—a critical test of re-entrainment theory

The spatial and size distribution of sediment deposited from short periods of overland flow due to the effect of a simulated grass buffer strip was measured for low slopes of 1.6, 3.4 and 5.1%. These data were analysed so as to critically evaluate two alternative models of the process of re-entrainment of recently deposited sediment. A model of re-entrainment, previously thought to be appropriate only for a steady-state or equilibrium situation, was found to give better agreement with experiments than did a model previously used in the literature on this subject.     

Drainage in finite-sized unsaturated zones

After the initiation of gravity drainage, water is often assumed to be either (a) draining under unit gradient, or (b) at capillary/gravity equilibrium. Both of these simplifications can be useful, but the regimes of validity of each assumption must be delineated. Water pressures are measured versus time and distance as water drains out of a 1.6 m long sand column to determine the relative effects of capillary and gravitational forces during drainage. For medium sized sands (0.15–0.3 mm in diameter), the capillary pressure is constant in space in a large region of the column for over 12 days, and the water continues to flow under unit gradient for relatively long time scales. Similar results are seen for finer sands, but with a much faster approach to equilibrium. Numerical simulations and analytical estimates are presented and compare favorably to the measurements. Together, the experimental, theoretical and analytical results are used to calculate when capillary/gravity equilibrium is reached as a function of porous media properties and length of the unsaturated zone. The ratio of the length of the unsaturated zone to the bubbling pressure is a key parameter in determining the drainage regime, and that even for relatively short unsaturated zones the equilibrium time scale can be on the order of years.     

Parameters affecting maximum fluid transport in large aperture fractures

We present results of laboratory experiments to study the behavior of liquids moving in unsaturated wide-aperture fractures. A 5-mm-thick glass plate cut with a 1.7-mm aperture was used as a fractured rock analog to study behavior of film and capillary droplet flow modes. Flow rates ranged between 0.6 and 6.0 ml/min. Variability in the ambient barometric pressure, resulting from weather conditions, seemed to play a role in the natural selection of flow mode. For droplet mode, constant input conditions resulted in highly variable transport properties within the fracture. The advancing meniscus exhibited a dynamic contact angle that was a function of the droplet speed and much larger than the static contact angle. Flow rate was reduced due to the larger contact angle. Analytical expressions for droplet velocity and flow capacity are presented as a function of the dynamic rather than the static contact angle.     

Laplace-transform analytic element solution of transient flow in porous media

A Laplace-transform analytic element method (LT-AEM) is described for the solution of transient flow problems in porous media. Following Laplace transformation of the original flow problem, the analytic element method (AEM) is used to solve the resultant time-independent modified Helmholtz equation, and the solution is inverted numerically back into the time domain. The solution is entirely general, retaining the mathematical elegance and computational efficiency of the AEM while being amenable to parallel computation. It is especially well suited for problems in which a solution is required at a limited number of points in space–time, and for problems involving materials with sharply contrasting hydraulic properties. We illustrate the LT-AEM on transient flow through a uniform confined aquifer with a circular inclusion of contrasting hydraulic conductivity and specific storage. Our results compare well with published analytical solutions in the special case of radial flow.     

Flow and fracture maps for basaltic rock deformation at high temperatures

Understanding how the strength of basaltic rock varies with the extrinsic conditions of stress state, pressure and temperature, and the intrinsic rock physical properties is fundamental to understanding the dynamics of volcanic systems. In particular it is essential to understand how rock strength at high temperatures is limited by fracture. We have collated and analysed laboratory data for basaltic rocks from over 500 rock deformation experiments and plotted these on principal stress failure maps. We have fitted an empirical flow law (Norton’s law) and a theoretical fracture criterion to these data. The principal stress failure map is a graphical representation of ductile and brittle experimental data together with flow and fracture envelopes under varying strain rate, temperature and pressure. We have used these maps to re-interpret the ductile–brittle transition in basaltic rocks at high temperatures and show, conceptually, how these failure maps can be applied to volcanic systems, using lava flows as an example.     

Petrography and uplift history of the Quaternary Takidani Granodiorite: could it have hosted a supercritical (HDR) geothermal reservoir?

The Quaternary Takidani Granodiorite (Japan Alps) is analogous to the type of deep-seated (3–5 km deep) intrusive-hosted fracture network system that might support (supercritical) hot dry/wet rock (HDR/HWR) energy extraction. The I-type Takidani Granodiorite comprises: porphyritic granodiorite, porphyritic granite, biotite-hornblende granodiorite, hornblende-biotite granodiorite, biotite-hornblende granite and biotite granite facies; the intrusion has a reverse chemical zonation, characterized by >70 wt% SiO₂ at its inferred margin and <67 wt% SiO₂ at the core. Fluid inclusion evidence indicates that fractured Takidani Granodiorite at one time hosted a liquid-dominated, convective hydrothermal system, with <380°C, low-salinity reservoir fluids at hydrostatic (mesothermal) pressure conditions. ‘Healed’ microfractures also trapped >600°C, hypersaline (35 wt% NaCl_eq) fluids of magmatic origin, with inferred minimum pressures of formation being 600–750 bar, which corresponds to fluid entrapment at 2.4–3.0 km depth. Al-in-hornblende geobarometry indicates that hornblende crystallization occurred at about 1.45 Ma (7.7–9.4 km depth) in the (marginal) eastern Takidani Granodiorite, but later (at 1.25 Ma) and shallower (6.5–7.0 km) near the core of the intrusion. The average rate of uplift across the Takidani Granodiorite from the time of hornblende crystallization has been 5.1–5.9 mm/yr (although uplift was about 7.5 mm/yr prior to 1.2 Ma), which is faster than average uplift rates in the Japan Alps (3 mm/yr during the last 2 million years). A temperature–depth–time window, when the Takidani Granodiorite had potential to host an HDR system, would have been when the internal temperature of the intrusive was cooling from 500°C to 400°C. Taking into account the initial (7.5 mm/yr) rate of uplift and effects of erosion, an optimal temperature–time–depth window is proposed: for 500°C at 1.54–1.57 Ma and 5.2±0.9 km (drilling) depth; and 400°C at 1.36–1.38 Ma and 3.3±0.8 km (drilling) depth, which is within the capabilities of modern drilling technologies, and similar to measured temperature–depth profiles in other active hydrothermal systems (e.g. at Kakkonda, Japan).