上海市长江口及邻近海域地质调查现状及展望

海洋区域地质调查近年来逐渐得到重视,海洋地质调查技术的发展也为开展海洋区域地质调查奠定了基础。上海地区社会经济发展对地质工作的需求,为摸清上海市长江口及邻近海域的地质现状提出了必然要求。本文在系统整理长江口及邻近海域地质资料的基础上,重点分析了区内存在的重点地质问题和研究存在的问题,根据海洋地质调查的技术发展,结合国内相关地区的调查经验,对长江口及邻近海域内区域地质调查进行了展望,提出了采用综合物探技术、地质取样和测试分析进行综合地质环境调查的建议,并对具体采用的调查技术和调查内容进行了论述。     

Dynamical evolution of comets and the problem of cometary fading

Possibilities to explain the observed 1/a-distribution are discussed in the light of improved understanding of the dynamical evolution of long-period comets. It appears that the fading problem applies both to single-injection and continuous-injection models. Although uncertainties due to nongravitational effects do not allow detailed results to be drawn from the observed 1/a-distribution at small perihelion distance q, that for q 1.5 AU shows that a constant fading probability cannot explain all the features of the observed distribution. Assuming that comets can reappear following a period of fading, values for the assumed constant fading and renewal probabilities, and the total cometary flux have been estimated for q > 1.5 AU.     

Infinitesimal orbits around lagrange points in the elliptic,restricted three-body problem

This study presents a method of obtaining asymptotic approximations for motions near a Lagrange point in the planar, elliptic, restricted three-body problem by using a von Zeipel-type method. The calculations are carried out for a second-order escape solution in the proximity of the equilateral Lagrange point, L ₄, where the primaries' orbital eccentricity is taken as the small parameter .     

Periodic orbits of an electric charge in a magnetic dipole field

Periodic orbits in the Stormer problem are studied using the symmetry lines of the Poincaré map introduced by De Vogelaere. Many known facts are explained by mean of these lines. The dynamics of four special symmetry lines when the Stormer parameter ₁ changes is presented, and we obtain a clear global view of the structure of the simple periodic orbits and their bifurcations, including the asymmetrical ones. New asymmetrical multiple periodic orbits are obtained.     

On Szebehely's problem for holonomic systems involving generalized potential functions

We consider the problem of finding the generalized potential function V = U _i(q ¹, q ²,..., q ⁿ)q ⁱ + U(q ¹, q ²,...;q ⁿ) compatible with prescribed dynamical trajectories of a holonomic system. We obtain conditions necessary for the existence of solutions to the problem: these can be cast into a system of n – 1 first order nonlinear partial differential equations in the unknown functions U ₁, U ₂,...;, U _n, U. In particular we study dynamical systems with two degrees of freedom. Using adapted coordinates on the configuration manifold M ₂ we obtain, for potential function U(q ¹, q ²), a classic first kind of Abel ordinary differential equation. Moreover, we show that, in special cases of dynamical interest, such an equation can be solved by quadrature. In particular we establish, for ordinary potential functions, a classical formula obtained in different way by Joukowsky for a particle moving on a surface.Work performed with the support of the Gruppo Nazionale di Fisica Matematica (G.N.F.M.) of the Italian National Research Council.     

Reservoir description using dynamic parameterisation selection with a combined stochastic and gradient search

There is a correspondence between flow in a reservoir and large scale permeability trends. This correspondence can be derived by constraining reservoir models using observed production data. One of the challenges in deriving the permeability distribution of a field using production data involves determination of the scale of resolution of the permeability. The Adaptive Multiscale Estimation (AME) seeks to overcome the problems related to choosing the resolution of the permeability field by a dynamic parameterisation selection. The standard AME uses a gradient algorithm in solving several optimisation problems with increasing permeability resolution. This paper presents a hybrid algorithm which combines a gradient search and a stochastic algorithm to improve the robustness of the dynamic parameterisation selection. At low dimension, we use the stochastic algorithm to generate several optimised models. We use information from all these produced models to find new optimal refinements, and start out new optimisations with several unequally suggested parameterisations. At higher dimensions we change to a gradient-type optimiser, where the initial solution is chosen from the ensemble of models suggested by the stochastic algorithm. The selection is based on a predefined criterion. We demonstrate the robustness of the hybrid algorithm on sample synthetic cases, which most of them were considered insolvable using the standard AME algorithm.     

Optimization of nonconventional wells under uncertainty using statistical proxies

The determination of the optimal type and placement of a nonconventional well in a heterogeneous reservoir represents a challenging optimization problem. This determination is significantly more complicated if uncertainty in the reservoir geology is included in the optimization. In this study, a genetic algorithm is applied to optimize the deployment of nonconventional wells. Geological uncertainty is accounted for by optimizing over multiple reservoir models (realizations) subject to a prescribed risk attitude. To reduce the excessive computational requirements of the base method, a new statistical proxy (which provides fast estimates of the objective function) based on cluster analysis is introduced into the optimization process. This proxy provides an estimate of the cumulative distribution function (CDF) of the scenario performance, which enables the quantification of proxy uncertainty. Knowledge of the proxy-based performance estimate in conjunction with the proxy CDF enables the systematic selection of the most appropriate scenarios for full simulation. Application of the overall method for the optimization of monobore and dual-lateral well placement demonstrates the performance of the hybrid optimization procedure. Specifically, it is shown that by simulating only 10% or 20% of the scenarios (as determined by application of the proxy), optimization results very close to those achieved by simulating all cases are obtained.     

Optimization of SAGD Well Elevation

The application of steam-assisted gravity drainage (SAGD) to recover heavy oil sands is becoming increasingly important in the northern Alberta McMurray Formation because of the vast resources/reserves accessible with this mechanism. Selecting the stratigraphic elevations of SAGD well pairs is a vital decision for reservoir evaluation and planning. The inherent uncertainty in the distribution of geological variables significantly influences this decision. Geostatistical simulation is used to capture geological uncertainty, which is used can be transformed into a distribution of the best possible well pair elevations. A simple exhaustive calculation scheme is used to determine the optimum stratigraphic location of a SAGD well pair where the recovery R is maximized. There are three basic steps to the methodology: (1) model the uncertainty in the top continuous bitumen (TCB) and bottom continuous bitumen (BCB) surfaces, (2) calculate the recovery at all possible elevation increments within the TCB and BCB interval, and (3) identify the elevation that maximizes R. This is repeated for multiple TCB/BCB pairs of surfaces to assess uncertainty. The methodology is described and implemented on a subset of data from the Athabasca Oilsands in Fort McMurray, Alberta.