用传输函数构建的大气重力波传播理论模式

本文根据考虑大气热传导和黏滞的重力波复色散关系,采用传输函数的概念,基于重力波的线性理论,构建了用于研究对流层内重力波激发源与电离层响应之间的传输函数数值模式.在相空间中讨论了传输函数振幅的分布特性,并以地面单位脉冲源为例,分析了从地面到300 km高空的响应,得到了物理量的时空分布特征.结果表明：（1）对内重力波的传播而言,大气相当于一个滤波器,只有波动周期在15～30 min,水平波长在200～450 km之间的重力波扰动才最容易到达300 km电离层高度;（2）电离层的响应主要在与地面的激发源之间相隔较远的水平距离上发生;（3）黏滞和热传导系数在低层对上传重力波的影响较小,随着高度的增加它们对重力波的影响越来越大;（4）在低层计算的波动频率与Row理论的计算结果比较一致,然而到了高层却相差较大.     

A consistent derivation of the wave-energy equation from basic hydrodynamic principles

Based on a decomposition of the velocity into mean flow, turbulent and wave components, momentum and hereafter a wave-energy equation is derived. It contains a turbulent energy dissipation term which is closed by applying a wave-related mixing length model and linear wave theory solutions. This closure produces a non-linear turbulent wave-energy dissipation including the wave energy in a 5/2 power law. The theory is able to predict correctly the shape of deep-water wave spectra according to Phillips' similarity law.Responsible Editor: Hans Burchard     

基于GSLS模型TI介质衰减拟声波方程

随着计算机硬件技术的发展以及高分辨率勘探需求的增加,我们希望能够更准确地模拟地下介质,得到更丰富的地层信息.然而,传统的声学假设并不能描述实际地层所存在各向异性和黏滞性,使得成像分辨率较低.为了实现深部储层的高精度成像,本文同时考虑了介质的各向异性和黏滞性,从TI介质弹性波的基本理论出发,结合各向异性GSLS理论,并通过声学近似方法导出基于GSLS模型的各向异性衰减拟声波方程.数值模拟表明该方程既能准确地描述各向异性介质下的准P波运动学规律,又能体现地层的吸收衰减效应;模型逆时偏移结果表明,在实现成像过程中考虑各向异性和黏滞性的影响,能对高陡构造清晰成像,且剖面振幅相对均衡,分辨率较高.     

基于连续方程的河道水流双变量耦合演算模型

针对现有的河道水流洪水演算模型只能模拟单一变量(流量或水位)的问题,以水流连续方程和河段蓄水量的两种不同表达形式(蓄水量等于平均过水断面面积与河段长乘积,蓄水量等于河段平均流量与传播时间的乘积)为基础,对马斯京根模型进行了通用性改进,提出了双变量耦合通用演算模型.选取了四大水系(包括内陆河流和入海河流)的16个河段汛期洪水资料进行模型检验,模型验证考虑了地理范围、不同的河段特征和水力特征、洪水量级等因素,全面地检验了模型结构的合理性和模拟实际洪水的有效性.将双变量耦合通用演算模型与传统的马斯京根法进行了效果比较,结果表明双变量耦合通用演算模型的模拟精度高于马斯京根法,模拟效果比马斯京根法稳定一些,而且具有较好的通用性.     

一个中纬电离层E层理论模式

本文建立了一个中纬电离层E层理论模式. 该模式从NO+,O２＋,O+和N２＋这四种主要离子的连续性方程出发,通过数值模拟得到中纬电离层E层电子和各种离子密度随时间和高度的变化情况. 计算结果能较好地反映出E层电子密度峰值（NmE）或E层临界频率（foE）的日变化、季节变化以及随太阳活动的变化趋势. 将模式的计算结果与武汉地区测高仪的观测数据进行比较,结果证明模式能够较为客观地反映中纬电离层E层的实际形态. 针对以往电离层E层理论模式存在的主要问题,本文还进一步讨论了几种重要因素,包括二次离化源,λ<150?谱段的辐射通量,吸收截面以及NO分布对于模式计算结果的影响.     

Changes of electrical conductivity in the vicinity of the curie temperature of basaltic rocks — theoretical model

Summary Samples of basalts, haematites and magnetites display either temporary or permanent (magnetite) weakening of the increase of the electrical conductivity with increasing temperature in the vicinity of their Curie temperature. Using the second quantization, this paper explains the observed pattern of the electrical conductivity adequately for magnetite and approximately for the other rocks by means of a quantum theory model of ferromagnetic minerals. This theory describes only the electron component of the electrical conductivity, which is responsible for the Curie temperature effect.Dedicated to RNDr. Jan Pícha, CSc., on his 60th Birthday     

Conductivity model for pyrite-bearing laminated and dispersed shaly sands based on a differential equation and the generalized Archie equation

The conductance of pyrite-bearing laminated and dispersed shaly sands is not well understood and resistivity models for pyrite-bearing shaly sands are nonexistent. Thus, we first synthesize clean pyrite-matrix samples, and quartz-matrix samples with variable laminated shale, dispersed shale, and pyrite content and then perform petrophysics experiments to assess the effect of pyrite content on the conductivity of pyrite-bearing shaly sands. Second, based on the differences in conductivity and conduction pathways and geometries because of the variable composition of the pyrite-bearing laminated and dispersed shaly sands, we divide the shaly sands into their components, i.e., laminated shale, quartz grains, pyrite grains, hydrocarbon, dispersed shale, microscopic capillary water, and mobile water. A generalized resistivity model is proposed to describe the conductivity of pyrite-bearing laminated and dispersed shaly sands, based on the combined conductivity differential equation and generalized Archie equation. In the generalized resistivity model, the conductivity differential equation is used to describe the conductivity of dispersed inclusions in a host, whereas the generalized Archie equation is used to describe the conductivity of two conducting phases. Moreover, parallel conductance theory is used to describe the conductivity of dispersed shaly sands and laminated shale. Theoretical analysis suggests that the proposed model satisfies the physical constraints and the model and experimental results agree. The resistivity and resistivity index of shaly sands decrease with increasing conductivity and pyrite. Finally, the accuracy of the resistivity model is assessed based on experimental data from 46 synthetic core samples with different oil saturation. The model can describe the conductivity of clean pyrite-matrix samples, and quartz-matrix samples with different volumes of laminated shale, dispersed shale, and pyrite. An accurate saturation model of pyrite-bearing laminated and dispersed shaly sands is thus obtained and the log data interpretation in complex shaly sands can improve with the proposed model.     

A two-stage model of fracture of rocks

In this paper we propose a two-stage model of rock fracture. In the first stage, cracks or local regions of failure are uncorrelated and occur randomly throughout the rock in response to loading of pre-existing flaws. As damage accumulates in the rock, there is a gradual increase in the probability that large clusters of closely spaced cracks or local failure sites will develop. Based on statistical arguments, a critical density of damage will occur where clusters of flaws become large enough to lead to larger-scale failure of the rock (stage two). While crack interaction and cooperative failure is expected to occur within clusters of closely spaced cracks, the initial development of clusters is predicted based on the random variation in pre-existing flaw populations. Thus the onset of the unstable second stage in the model can be computed from the generation of random, uncorrelated damage. The proposed model incorporates notions of the kinetic (and therefore time-dependent) nature of the strength of solids as well as the discrete hierarchic structure of rocks and the flaw populations that lead to damage accumulation. The advantage offered by this model is that its salient features are valid for fracture processes occurring over a wide range of scales including earthquake processes. A notion of the rank of fracture (fracture size) is introduced, and criteria are presented for both fracture nucleation and the transition of the failure process from one scale to another.     

A constitutive model for water flow in unsaturated fractured rocks

A conceptual model for describing effective saturation in fractured hard rock is presented. The fracture network and the rock matrix are considered as an equivalent continuum medium where each fracture is conceptualized as a porous medium of granular structure and the rock matrix is assumed to be impermeable. The proposed model is based on the representation of a rough‐walled fracture by an equivalent porous medium, which is described using classical constitutive models. A simple closed‐form equation for the effective saturation is obtained when the van Genuchten model is used to describe saturation inside fractures and fractal laws are assumed for both aperture and number of fractures. The relative hydraulic conductivity for the fractured rock is predicted from a simple relation derived by Liu and Bodvarsson. The proposed constitutive model contains three independent parameters, which may be obtained by fitting the proposed effective saturation curve to experimental data. Two of the model parameters have physical meaning and can be identified with the reciprocal of the air entry pressure values in the fractures of minimum and maximum apertures. Effective saturation and relative hydraulic conductivity curves match fairly well the simulated constitutive relations obtained by Liu and Bodvarsson. Copyright © 2008 John Wiley & Sons, Ltd.     

The influence of dry and water saturated cracks on seismic velocities of crustal rocks - A comparison of experimental data with theoretical model

The pressure dependence of P- and S-wave velocities, velocity anisotropy, shear wave splitting and crack-porosity has been investigated in a number of samples from different crustal rock types for dry and wet (water saturated) conditions. At atmospheric pressure, P-wave velocities of the saturated, low-porosity rocks (< 1%) are significantly higher than in dry rocks, whereas the differences for S-wave velocities are less pronounced. The effect of intercrystalline fluids on seismic properties at increased pressure conditions is particularly reflected by the variation of the Poisson's ratio because P-wave velocities are more sensitive to fluids than S-wave velocities in the low-porosity rocks. Based on the experimental data, the respective crack-density parameter (), which is a measure of the number of flat cracks per volume unit contained within the background medium (crack-free matrix), has been calculated for dry and saturated conditions. There is a good correlation between the calculated crack-densities and crack-porosities derived from the experimentally determined volumetric strain curves. The shear wave velocity data, along with the shear wave polarisation referred to a orthogonal reference system, have been used to derive the spatial orientation of effective oriented cracks within a foliated biotite gneiss. The experimental data are in reasonable agreement with the self consistent model of O'Connell and Budiansky (1974). Taking the various lithologies into account, it is clear from the present study, that combined seismic measurements ofV _p andV _s , using theV _p V _s -ratio, may give evidence for fluids on grain boundaries and, in addition, may provide an estimate on the in-situ crack-densities.     

A theoretical model of water and trade

Water is an essential input for agricultural production. Agriculture, in turn, is globalized through the trade of agricultural commodities. In this paper, we develop a theoretical model that emphasizes four tradeoffs involving water-use decision-making that are important yet not always considered in a consistent framework. One tradeoff focuses on competition for water among different economic sectors. A second tradeoff examines the possibility that certain types of agricultural investments can offset water use. A third tradeoff explores the possibility that the rest of the world can be a source of supply or demand for a country’s water-using commodities. The fourth tradeoff concerns how variability in water supplies influences farmer decision-making. We show conditions under which trade liberalization affect water use. Two policy scenarios to reduce water use are evaluated. First, we derive a target tax that reduces water use without offsetting the gains from trade liberalization, although important tradeoffs exist between economic performance and resource use. Second, we show how subsidization of water-saving technologies can allow producers to use less water without reducing agricultural production, making such subsidization an indirect means of influencing water use decision-making. Finally, we outline conditions under which riskiness of water availability affects water use. These theoretical model results generate hypotheses that can be tested empirically in future work.     

基于粒子群优化的理论变异函数拟合方法研究

变异函数是地统计学中区域化变量空间结构分析和空间局部插值的主要分析工具.理论变异函数模型的获取是地质统计学中的基础性工作,它是了解区域化变量的变异特征、进一步对地质统计学计算的必要环节.针对现有的理论变异函数的拟合方法,如人工拟合法、线性规划拟合法、加权多项式拟合法、目标规划拟合法等的不足之处,充分利用粒子群优化算法在求解非线性优化问题时具有的全局寻优的特点,提出基于粒子群优化的理论变异函数拟合方法.在实例应用中,分别利用粒子群优化算法和加权多项式拟合方法进行理论变异函数拟合,交叉验证结果表明粒子群优化算法预测精度较高,具有较强的稳健性.     

基于变形P-L模型的矩阵方程迭代精细横波预测

横波速度预测问题的关键有两个,一是如何建立合理的岩石物理模型,二是针对建立的横波预测目标函数,如何准确高效地求解.针对第一个问题,对Pride模型和Lee模型(P-L模型)进行变形,提出拟固结指数的概念,将干岩石模量和岩石基质模量相联系,变形后的P-L模型在没有降低P-L模型准确度的情况下简化了问题的复杂度,建立起了饱和流体岩石弹性模量与干岩石模量、岩石基质模量、混合流体模量之间的关系,进而计算理论上的纵波速度,并通过比较实测纵波速度与计算的理论纵波速度大小,最终建立了横波预测的目标函数.针对第二个问题,借鉴了地震反演的思路,将该目标函数的最优化问题转化为线性矩阵方程组迭代求解问题,通过几步迭代就可以求解出合适的拟固结指数,进而得到预测横波速度.实际验证和应用表明,该横波预测方法具有很好的稳定性和准确性,并且岩石物理模型的构建和目标函数的求解思路可用于其他储集类型地层的横波预测.     

Hydrodynamic derivation of storage parameters of the Muskingum model

The St. Vénant equations for unsteady flow in open channels and the Muskingum method are written both in their conventional forms and in the state-space formulation. The hydrodynamic equation of motion is solved by the method of state trajectory variation and the result for the first-order variation in the state-space variables is used as a basis of linking the parameters of the Muskingum model with the hydraulic parameters of the open channel reach. The results are applicable to any shape of cross-section and to any type of friction law.     

A multifractal model of crack coalescence in rocks

Geometric characteristics of fractal sets of cracks are investigated from the standpoint of their coalescence. A 1D computer model, in which a realistic character of forms of crack sets is achieved through the inducing of sets with the use of the multiplicative cascade procedure, is considered. The investigation is aimed at the determination of conditions of coalescence of cracks organized as a superposition of fractal subsets of a unified set of cracks. The process of crack coalescence is investigated for three different values of the coalescence criterion. The geometric characteristics of the sets of cracks that are necessary for the transition of the crack coalescence process into an avalanche-like stage are estimated.