山东省大中型水库水化学性状研究──几种主要水化因子的变化及其在评价水库鱼产力中的作用

本文通过对山东省３８座代表性大、中型水库总氮（ＴＮ）、总磷（ＴＰ）、化学耗氧量、硬度及电导率（λ）五种主要水化因子连续两年（１９８９～１９９０）的调查与研究，探讨了水化因子与水库鱼产力的关系。用反映鲢、鳙生长状况的综合生长指数（ＧＩ）作为水库鱼产力的指标。它与水化因子相关分析表明，总磷是山东省３８座水库鱼产力的主要影响因子。提出ＧＩ与总磷、总氮及电导率等水化因子的逐步回归方程关系非常显著。证明利用水化因子作为水库鱼产力的评价指标是可行的。并指出：同一水库的相同水化因子存在年间差异，其主要原因在于降雨量和水库进水量的不同，鉴于此．同一水库鱼产力也存在年间差异。     

The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C

The recent improvements in the Gravity Recovery And Climate Experiment (GRACE) tracking data processing at GeoForschungsZentrum Potsdam (GFZ) and Groupe de Recherche de Géodésie Spatiale (GRGS) Toulouse, the availability of newer surface gravity data sets in the Arctic, Antarctica and North-America, and the availability of a new mean sea surface height model from altimetry processing at GFZ gave rise to the generation of two new global gravity field models. The first, EIGEN-GL04S1, a satellite-only model complete to degree and order 150 in terms of spherical harmonics, was derived by combination of the latest GFZ Potsdam GRACE-only (EIGEN-GRACE04S) and GRGS Toulouse GRACE/LAGEOS (EIGEN-GL04S) mean field solutions. The second, EIGEN-GL04S1 was combined with surface gravity data from altimetry over the oceans and gravimetry over the continents to derive a new high-resolution global gravity field model called EIGEN-GL04C. This model is complete to degree and order 360 and thus resolves geoid and gravity anomalies at half- wavelengths of 55 km at the equator. A degree-dependent combination method has been applied in order to preserve the high accuracy from the GRACE satellite data in the lower frequency band of the geopotential and to form a smooth transition to the high-frequency information coming from the surface data. Compared to pre-CHAMP global high-resolution models, the accuracy was improved at a spatial resolution of 200 km (half-wavelength) by one order of magnitude to 3 cm in terms of geoid heights. The accuracy of this model (i.e. the commission error) at its full spatial resolution is estimated to be 15 cm. The model shows a reduced artificial meridional striping and an increased correlation of EIGEN-GL04C-derived geostrophic meridional currents with World Ocean Atlas 2001 (WOA01) data. These improvements have led to select EIGEN-GL04C for JASON-1 satellite altimeter data reprocessing. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Methodology and use of tensor invariants for satellite gravity gradiometry

Although its use is widespread in several other scientific disciplines, the theory of tensor invariants is only marginally adopted in gravity field modeling. We aim to close this gap by developing and applying the invariants approach for geopotential recovery. Gravitational tensor invariants are deduced from products of second-order derivatives of the gravitational potential. The benefit of the method presented arises from its independence of the gradiometer instrument’s orientation in space. Thus, we refrain from the classical methods for satellite gravity gradiometry analysis, i.e., in terms of individual gravity gradients, in favor of the alternative invariants approach. The invariants approach requires a tailored processing strategy. Firstly, the non-linear functionals with regard to the potential series expansion in spherical harmonics necessitates the linearization and iterative solution of the resulting least-squares problem. From the computational point of view, efficient linearization by means of perturbation theory has been adopted. It only requires the computation of reference gravity gradients. Secondly, the deduced pseudo-observations are composed of all the gravitational tensor elements, all of which require a comparable level of accuracy. Additionally, implementation of the invariants method for large data sets is a challenging task. We show the fundamentals of tensor invariants theory adapted to satellite gradiometry. With regard to the GOCE (Gravity field and steady-state Ocean Circulation Explorer) satellite gradiometry mission, we demonstrate that the iterative parameter estimation process converges within only two iterations. Additionally, for the GOCE configuration, we show the invariants approach to be insensitive to the synthesis of unobserved gravity gradients.     

Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition

We develop and apply an efficient strategy for Earth gravity field recovery from satellite gravity gradiometry data. Our approach is based upon the Paige-Saunders iterative least-squares method using QR decomposition (LSQR). We modify the original algorithm for space-geodetic applications: firstly, we investigate how convergence can be accelerated by means of both subspace and block-diagonal preconditioning. The efficiency of the latter dominates if the design matrix exhibits block-dominant structure. Secondly, we address Tikhonov-Phillips regularization in general. Thirdly, we demonstrate an effective implementation of the algorithm in a high-performance computing environment. In this context, an important issue is to avoid the twofold computation of the design matrix in each iteration. The computational platform is a 64-processor shared-memory supercomputer. The runtime results prove the successful parallelization of the LSQR solver. The numerical examples are chosen in view of the forthcoming satellite mission GOCE (Gravity field and steady-state Ocean Circulation Explorer). The closed-loop scenario covers 1 month of simulated data with 5 s sampling. We focus exclusively on the analysis of radial components of satellite accelerations and gravity gradients. Our extensions to the basic algorithm enable the method to be competitive with well-established inversion strategies in satellite geodesy, such as conjugate gradient methods or the brute-force approach. In its current development stage, the LSQR method appears ready to deal with real-data applications.     

Study on Recovering the Earth＇s Potential Field Based on GOCE Gradiometry

Given the second radial derivative Vrr（P） ｜δs of the Earth＇s gravitational potential V（P） on the surface δS corresponding to the satellite altitude, by using the fictitious compress recovery method, a fictitious regular harmonic field rrVrr（P）^＊ and a fictitious second radial gradient field V：（P） in the domain outside an inner sphere Ki can be determined, which coincides with the real field V（P） in the domain outside the Earth. Vrr^＊（P）could be further expressed as a uniformly convergent expansion series in the domain outside the inner sphere, because rrV（P）^＊ could be expressed as a uniformly convergent spherical harmonic expansion series due to its regularity and harmony in that domain. In another aspect, the fictitious field V^＊（P） defined in the domain outside the inner sphere, which coincides with the real field V（P） in the domain outside the Earth, could be also expressed as a spherical harmonic expansion series. Then, the harmonic coefficients contained in the series expressing V^＊（P） can be determined, and consequently the real field V（P） is recovered. Preliminary simulation calculations show that the second radial gradient field Vrr（P） could be recovered based only on the second radial derivative V（P）｜δs given on the satellite boundary. Concerning the final recovery of the potential field V（P） based only on the boundary value Vrr （P）｜δs, the simulation tests are still in process.     

三峡库区高切坡致灾因素及防护对策——以重庆市万州区为例

三峡库区百万移民的搬迁建设导致大量高切坡的形成。未及时治理或治理不当的高切坡,不同程度地发生破坏,影响了人们的正常生活,甚至引发了地质灾害。因此,科学认识及治理高切坡十分必要。重庆市万州区高切坡主要有风化及差异风化、崩塌和土体滑塌等破坏形式。影响破坏的因素主要为工程地质条件、气候及时间效应。针对高切坡不同的破坏形式,以安全、经济及环境和谐的人文理念作为设计指导思想,充分论证了高切坡的防护方案,即：削坡、排水以及格构防护较适合万州区高切坡的防护。这些措施,既能很好地解决高切坡的防护问题,又能兼顾环保、降低造价。研究结果可以作为三峡地区高切坡类地质灾害防护设计的参考。也可为其它流域类似问题借鉴。     

三峡库区消落区土壤磷释放的环境影响因子

以三峡库区消落区万州段为试验基地,选取释磷能力较强的紫色冲积土,根据三峡水库消落带的"干湿交替"空间和时间特征,进行万州断面土壤磷释放的影响因子的实验室模拟试验和万州江面淹没的对照试验。研究发现,随着淹水时间的变化,各种形态的磷在前5~10d各形态磷变化有相当的差异,随后变化趋势趋于稳定。TP有降低的趋势,土壤磷有一定的释放。野外研究表明Olsen-P在淹水10周左右到达最大值,以后缓慢降低。淹水时土壤Olsen-P增加,干燥后降低。多次淹水时每次淹水后土壤的有效磷水平都略有增加,落干后相较上次落干后有效磷水平降低,最后一次淹水后相较与初始时Olsen-P水平低16.7%。随着淹水深度增加,土壤的Olsen-P水平在淹水时由淹水前20.53mg/kg提高到43.23mg/kg,增加110.6%。当上覆水磷浓度较低时(<2mg/L),磷吸附到达平衡的时间较短,约需要6周。当上覆水磷浓度较高(>2mg/L)时,磷吸附到达平衡的时间较长。微生物活动对淹水土壤的磷释放有一定影响,有微生物时磷的释放高于无微生物者0.048mg/L。种植植物的土壤在淹水后Olsen-P含量大于土壤直接淹没时的释放量,种植狗牙根(Cynodondactylon)和野地瓜藤(Ficustikoua)的土壤中Olsen-P分别较未种植物土壤释放量高出21.5%和12.7%。     

Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.     

渤海湾盆地埕岛油田缝洞型潜山油气藏构造及储层特征

埕岛油田古潜山油气藏为典型的缝洞型碳酸盐岩油气藏。埕岛地区经历了多期次的构造运动,并经受了多期次的断裂与剥蚀作用,为油气成藏创造了良好条件。埕岛油田古潜山可以划分为三个构造单元,即埕北20断层下降盘的西排山,上升盘的中排山和东排山,三个潜山带成帚状敛于凸起最南端,为几个二级构造单元交汇处,构造应力较发育。埕岛油田有三种储集空间类型:晶间孔隙、溶蚀孔洞和裂缝。     

鄂尔多斯盆地靖边气田马五_(1+2)气藏储集单元研究

综合地震、地质、测井等信息,考虑沉积、成岩、地貌、生产动态等多因素,在沟槽识别和致密带研究基础上,通过井间干扰、油层压力变化和试井边界响应分析,按照局部、就近原则对靖边气田马五1+2气藏进行了储集单元划分。共划分出111个储集单元,其中独立型69个,复合型42个。这些储集单元的面积多介于10～50km2之间,单个储集单元的井数一般少于6口。在对气藏储层特征综合认识的基础上,确定流动单元划分思路与流程,并认为关键环节包括古沟槽识别、致密带研究、动静结合确定储集单元等3个方面。储集单元一方面受控于被致密岩相包围的有利成岩相带,宏观上与弱—中等充填泥粉晶白云岩相对应性好;另一方面又受各级古沟槽夹持,在两个一级古沟槽之间往往共生多个储集单元,而这些储集单元边界又严格地受二级或三级沟槽切割控制。