Two solution methods for dynamic game in reservoir operation

During the last decade, a number of models have been developed to consider the conflict in dynamic reservoir operation. Most of these models are discrete dynamic models which are developed based on game theory. In this study, a continuous model of dynamic game and its corresponding solutions are developed for reservoir operation. Two solution methods are used to solve the model of continuous dynamic game, namely the Ricatti equations and collocation methods. The Ricatti equations method is a closed form solution, requiring less computational efforts compared with discrete models. The collocation solution method applies Newton's method or a quasi-Newton method to find the problem solution. These approaches are able to generate operating policies for dynamic reservoir operation. The Zayandeh-Rud river basin in central Iran is used as a case study and the results are compared with alternative water allocation models. The results show that the proposed solution methods are quite capable of providing appropriate reservoir operating policies, while requiring rather short computational times due to continuous formulation of state and decision variables. Reliability indices are used to compare the overall performance of the proposed models. Based on the results from this study, the collocation method leads to improved values of the reliability indices for total reservoir system and utility satisfaction of water users, compared to the Ricatti equations method. This is attributed to the flexible structure of the collocation model. When compared to alternative water allocation models, lower values of reliability indices are achieved by the collocation method.     

Lattice Boltzmann simulations of the transient shallow water flows

A two-dimensional lattice Boltzmann model (LBM) is presented for transient shallow water flows. The model is based on the shallow water equations coupled with the large eddy simulation model. In order to obtain accurate results efficiently, a multi-block lattice scheme is applied at the area where a local finer grid is needed for strong change in physical variables. The model is verified by applying to five cases with transient processes: (a) a tidal wave over steps; (b) a perturbation over a submerged hump; (c) partial dam break flow; (d) circular dam break flow; (e) interaction between a dam break surge and four square cylinders. The objectives of this study are to validate the two-dimensional LBM in transient flow simulation and provide the detailed transient processes in shallow water flows.     

二维均匀介质中积分方程的快速精确算法

采用稳定型双共轭梯度快速Fourier变换(BCGS-FFT)算法精确计算二维均匀介质中的积分方程.采用一种新的插值函数作为基函数和试探函数对积分方程进行弱化离散,离散后的积分方程采用稳定型双共轭梯度迭代方法进行求解,从而得到异常体内电场的分布.计算时采用快速Fourier变换技术将积分方程内Green函数与电场的乘积表示成褶积形式以加快计算速度.数值计算举例说明了算法的精确性和有效性.     

线性方程组在ＧＰＳ坐标转换中的应用

手持GPS的出现给各个行业导航定位带来了极大的便利。而不同的坐标系统需要在使用中随时变化,利用统计法线性方程组求解手持GPS平面坐标参数,是一条解决手持GPS坐标系统的便捷方法。     

Ground-motion prediction equations for the intermediate depth Vrancea (Romania) earthquakes

We present the regional ground-motion prediction equations for peak ground acceleration (PGA), peak ground velocity (PGV), pseudo-spectral acceleration (PSA), and seismic intensity (MSK scale) for the Vrancea intermediate depth earthquakes (SE-Carpathians) and territory of Romania. The prediction equations were constructed using the stochastic technique on the basis of the regional Fourier amplitude spectrum (FAS) source scaling and attenuation models and the generalised site amplification functions. Values of considered ground motion parameters are given as the functions of earthquake magnitude, depth and epicentral distance. The developed ground-motion models were tested and calibrated using the available data from the large Vrancea earthquakes. We suggest to use the presented equations for the rapid estimation of seismic effect after strong earthquakes (Shakemap generation) and seismic hazard assessment, both deterministic and probabilistic approaches.     

双相各向异性介质中偶数阶精度有限差分数值模拟

To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.     

An Automatic Method of Direct Interpretation of Residual Gravity Anomaly Profiles due to Spheres and Cylinders

We have developed a least-squares minimization approach to determine the depth and the amplitude coefficient of a buried structure from residual gravity anomaly profile. This approach is basically based on application of Werner deconvolution method to gravity formulas due to spheres and cylinders, and solving a set of algebraic linear equations to estimate the two-model parameters. The validity of this new method is demonstrated through studying and analyzing two synthetic gravity anomalies, using simulated data generated from a known model with different random error components and a known statistical distribution. After being theoretically proven, this approach was applied on two real field gravity anomalies from Cuba and Sweden. The agreement between the results obtained by the proposed method and those obtained by other interpretation methods is good and comparable. Moreover, the depth obtained by the proposed approach is found to be in very good agreement with that obtained from drilling information.     

A weighted surface-depth gradient method for the numerical integration of the 2D shallow water equations with topography

A finite volume MUSCL scheme for the numerical integration of 2D shallow water equations is presented. In the framework of the SLIC scheme, the proposed weighted surface-depth gradient method (WSDGM) computes intercell water depths through a weighted average of DGM and SGM reconstructions, in which the weight function depends on the local Froude number. This combination makes the scheme capable of performing a robust tracking of wet/dry fronts and, together with an unsplit centered discretization of the bed slope source term, of maintaining the static condition on non-flat topographies (C-property). A correction of the numerical fluxes in the computational cells with water depth smaller than a fixed tolerance enables a drastic reduction of the mass error in the presence of wetting and drying fronts. The effectiveness and robustness of the proposed scheme are assessed by comparing numerical results with analytical and reference solutions of a set of test cases. Moreover, to show the capability of the numerical model on field-scale applications, the results of a dam-break scenario are presented.     

High order numerical code for hyperbolic mild-slope equations with nonlinear dispersion relation

Based on the hyperbolic mild-slope equations derived by Copeland (1985), a numerical model is established in unstag- gered grids. A composite 4 th-order Adam-Bashforth-Moulton (ABM) scheme is used to solve the model in the time domain. Terms involving the first order spatial derivatives are differenced to O ( Δx )4accuracy utilizing a five-point formula. The nonlinear dispersion relationship proposed by Kirby and Dalrymple (1986) is used to consider the nonlinear effect. A numerical test is performed upon wave propagating over a typical shoal. The agreement between the numerical and the experimental results validates the present model. Biodistribution and applications are also summarized.     

绿洲生态系统的发展和气候的关系

在植被和裸地共存的情况下,考虑了植物的某些生态条件,利用热量平衡方程,对于不同气候和生态条件计算了绿洲总的蒸散率以及植被和土壤的温度。研究表明,在某些条件下,蒸散率的解可出现准分岔和多平衡态的形式,它们主要发生在植被覆盖小的地区。