Spatial averaging of hydraulic conductivity in three-dimensional heterogeneous porous media

Numerical models of groundwater flow require the assignment of hydraulic conductivities to large grid blocks discretizing the flow domain; however, conductivity data is usually available only at the much smaller scale of core samples. This paper describes a geostatistical model for hydraulic conductivity at both the core or point scale and that of grid blocks. Conductivity at the block scale is obtained empirically as a spatial power-average of point scale values. Assuming a multivariate Gaussian model for point log-conductivity, expressions are derived for the ensemble mean and variance of block conductivity. The expression for the ensemble mean of block scale conductivity is found to be similar to an expression for the ensemble effective conductivity of an infinite field derived analytically by earlier authors. Here, block conductivities obtained by power averaging are compared with effective conductivities obtained from a numerical flow model and are found to be in excellent agreement for a suitably chosen averaging exponent. This agreement deteriorates gradually as the log variance of conductivity increases beyond 2. For arbitrary flow field geometry and anisotropic conductivity covariances, the averaging exponent can be calibrated by recourse to numerical flow experiments. For cubic fields and an isotropic spatial covariance, the averaging exponent is found to be 1/3. In this particular case, it was found that flow field discretization at the block scale through local averaging of point conductivities gave similar results to those obtained directly using a point scale discretization of the flow field.     

变动流域滞时方程

流域滞时是流域汇流计算和预报中的一个重要因素。本文尝试将流域滞时表示成流域重心长度Lca(流域重心到流域出口断面的河道长度)与汇流速度v之比的函数,从而使流域滞时随汇流速度的不同而不同,而不再是一个常数。长度Lca反映了流域尺度和几何形态对滞时的影响,而汇流速度则反映流域动力条件对滞时的作用。5个小流域的实测资料检验表明结果较为满意。     

方差-协方差分量估计公式等价性的进一步证明

通过附有条件的间接平差模型进一步证明各类方差-协方差分量估计公式之间的等价性。     

由直线和曲线混合组成的地形符号在DCS中的处理

在ＤＣＳ输出的地籍图上，有许多由直线段和曲线段组成的地物符号，用传统方法处理难以取得理想的绘图效果。本文在分析ＤＣＳ数据结构和属性码的基础上，提出了曲线部分按三次参数样条曲线逼近，讨论了直线段和曲线段的各种组合情况，并且通过正确地选择边值条件使二者之间．的连接比较符合实际。最后扼要分析了该类线段按虚线绘出时的处理方法。     

中低纬过渡地区潜热影响下中尺度扰动的对称发展

本文引入第二科里奥利参数讨论了中低纬过渡地区含非绝热(潜热)作用的对称扰动波包发展问题。结果表明:这一地区扰动波包的发展不仅包含了中高纬度波包发展条件［1］，而且与潜热和基本风场梯度经向切变及波包传播的非垂直性有关。过渡地区基本风场的作用是通过第一和第二科里奥利参数完成的。     

VENTILATION FLOW IN A BAROCLINIC VORTEX RELATED TO TROPICAL CYCLONE MOTION

The tropical cyclone motion is numerically simulated with a quasi-geostrophic baroclinic model.The flow field of a tropical cyclone is decomposed into its axisymmetric and asymmetric components.The relation between the ventilation flow vector and the motion vector of the tropical cyclone is inves-tigated.The results of numerical experiments indicate:(1) There are both large-scale beta gyres andsmall-scale gyres in the asyrnmetric flow field.(2) The interaction between small-scale gyres andlarge-scale beta gyres leads to the oscillation of translation speed and translation direction for the tropi-cal cyclone.(3) There are the large deviations between the ventilation flow vector calculated bymeans of Fiorino and Elsberry's method and the motion vector of tropical cyclone.(4) The ventila-tion flow vector computed using the improved method closely correlates with the motion vector of thetropical cyclone.     

二苯硫脲－甲基异丁酮泡塑富集液珠萃取比色法测定化探样品中的金银

在１０％ＨＣＩ介质中，采用二苯硫脲（ＤＰＴＵ）－甲基异丁酮（ＭＩＢＫ）泡沫塑料富集，以金试剂液珠萃取比色法测定了化探样品中的Ａｕ和Ａｇ。方法的检出限：Ａｕ０．５×１０￣－９；Ａｓ５×１０￣－９（１０ｇ称样）。利用此方法制备了“野外快速分析箱”，可在野外简陋条件下进行现场测定。     

Deep crustal growth of quartz, kyanite and garnet into large-aperture, fluid-filled fractures, north-eastern Connecticut, USA

Abstract A petrographic and petrological analysis of exceptionally well-preserved hydrothermal veins from the Merrimack synclinorium, north-eastern Connecticut, has been carried out in order to place new field-based constraints on fracture aperture dimensions and porosity in the lower continental crust. The veins preserve substantial open space today in outcrop, and contain mineral assemblages including subhedral to euhedral crystals of quartz, kyanite and almandine-rich garnet. Textural evidence indicates unequivocally that the vein minerals grew into macroscopic (mm- to cm-scale) open space between the vein walls. The veins are interpreted to have been large-aperture fractures along which significant advective fluid infiltration and chemical reaction occurred. The porosity of the rock mass due to open space between fracture walls today is c . 0.3%, but it could have been as large as several percent when the flow system was active. Quantitative thermobarometry results from vein mineral assemblages indicate that the fractures formed at pressures and corresponding crustal depths of c. 0.8 GPa and c. 30km, and temperatures of 550–600° C. The depth of fracture formation corresponds to published estimates of the maximum burial depth of the Merrimack synclinorium during the Acadian orogeny. The formation of large-aperture fractures could increase significantly the transient permeability of the deep crust, and therefore influence metamorphic heat and mass transfer.     

Strain-related differences in the crystal growth of white mica and chlorite: a TEM and XRD study of the development of metapelitic microfabrics in the Southern Uplands thrust terrane, Scotland

Abstract TEM and XRD techniques were used to study crystal growth characteristics of the fabric-forming phyllosilicates which developed in response to low-grade metamorphism and tectonic imbrication in part of the Southern Uplands thrust terrane. Prograde regional metamorphism, ranging from late diagenesis through the anchizone to the epizone, was accompanied by the development of a slaty cleavage which is commonly bedding-parallel. TEM-measured mean thicknesses of white mica and chlorite crystallite populations increase with advancing grade and correlate with XRD-measured crystallinity indices. Analytical TEM data show that prograde changes in composition lead to a net loss of Si, Ca and minor Fe from the fabric-forming phyllosilicates. White micas are paragonite-poor phengites with a mean b lattice parameter of 9.037 Å, and indicate an intermediate pressure series of metamorphism with a field gradient of <25° C km^-1. Chlorite compositions evolved from diabantite (with intergrown corrensite) to ripidolite over an estimated temperature range of 150–320° C. Field gradient and temperature estimates suggest that crystal growth and fabric development occurred at burial depths ranging from 6 km to at least 13 km in the thrust terrane. During late diagenesis, crystal growth of white mica and chlorite was predominantly a consequence of polytypic and phase transitions, and resulted in similar size distributions which resemble typical Ostwald ripening curves. Under anchizonal and epizonal conditions, white mica grew more rapidly than chlorite because of its greater ability to store strain energy and recover from subgrain development; as a result crystal thickness distributions are not typical of Ostwald ripening. In contrast, chlorite crystals which grew under these conditions developed subgrain boundaries at high strain rates which were only partially recovered at low strain rates; these retained dislocations reduce the crystallite thicknesses detected by TEM and XRD, compared with those of white mica. These differences in strain-induced crystal growth indicate that white mica (illite) and chlorite crystallinity indices are likely to show significant differences where low-grade metamorphism is closely associated with tectonic fabric development.     

Generalized cross-validation for covariance model selection

A weighted cross-validation technique known in the spline literature as generalized cross-validation (GCV), is proposed for covariance model selection and parameter estimation. Weights for prediction errors are selected to give more importance to a cluster of points than isolated points. Clustered points are estimated better by their neighbors and are more sensitive to model parameters. This rational weighting scheme also provides a simplifying significantly the computation of the cross-validation mean square error of prediction. With small- to medium-size datasets, GCV is performed in a global neighborhood. Optimization of usual isotropic models requires only a small number of matrix inversions. A small dataset and a simulation are used to compare performances of GCV to ordinary cross-validation (OCV) and least-squares filling (LS).