基于空间自相关的阿根廷滑柔鱼CPUE标准化研究

CPUE的观测往往不是独立的,而是存在空间相关性的。但是,大多数的CPUE标准化方法通常都假设名义CPUE在空间上是相互独立的。为此,本研究以西南大西洋阿根廷滑柔鱼为例,采用2000-2014年1-5月中国大陆鱿钓生产统计数据以及对应的海表温度和叶绿素浓度数据,选择广义线性模型（general linear model,GLM）为基础模型,将空间自相关加入到GLM中,比较标准GLM和4种加入空间自相关的空间GLM的CPUE标准化。根据最小信息准则（Akaike Information Criterion,AIC）及贝叶斯信息准则（Bayesian Information Criterion,BIC）,空间自相关的GLM的CPUE标准化结果优于标准GLM,其中指数模型的CPUE标准化结果最佳。同时,标准GLM与空间自相关的GLM相比,存在精确度过高估计的问题。因此,在CPUE标准化中,应充分考虑空间自相关这一因素。     

大别山及其周边现今垂直形变研究

基于1980年、2005年和2015年3期区域精密水准网观测资料,利用线性动态平差模型计算获取大别山地区水准网长期垂直运动速度场图像。研究发现,淮河平原地表下沉较为严重,大别山呈现弱隆升趋势运动,长江谷地边缘地区较大别山区呈现明显的隆升运动。跨郯庐断裂带水准剖面结果显示,垂直运动与地形呈负相关和弱相关。     

双线偏振雷达判别降水粒子类型技术及其检验

在综合分析国内外降水粒子相态识别方法的基础上,结合我国双线偏振雷达的应用需求,建立了一套利用双线偏振雷达探测资料识别降水粒子类型的识别模式,并根据美国KOUN雷达观测资料对该模式进行了分析验证,同时检验了模式中的几个偏振参数:差分反射率因子ZDR、差传播相移KDP、水平偏振和垂直偏振回波功率零滞后互相关系数ρHV(0)对不同降水粒子的识别效果。通过分析认为:利用该识别模式得到的结果基本合理,可以反映降水区域内降水粒子的相态结构,但还需要对资料作进一步的研究。同时在识别各类降水粒子方面,ZDR的效果要高于KDP和ρHV(0),在缺乏相关资料的情况下可以考虑用ZH和ZDR配对识别的结果作为真实的结果。     

Investigation of highly efficient algorithms for solving linear equations in the discontinuous deformation analysis method

Large‐scale engineering computing using the discontinuous deformation analysis (DDA) method is time‐consuming, which hinders the application of the DDA method. The simulation result of a typical numerical example indicates that the linear equation solver is a key factor that affects the efficiency of the DDA method. In this paper, highly efficient algorithms for solving linear equations are investigated, and two modifications of the DDA programme are presented. The first modification is a linear equation solver with high efficiency. The block Jacobi (BJ) iterative method and the block conjugate gradient with Jacobi pre‐processing (Jacobi‐PCG) iterative method are introduced, and the key operations are detailed, including the matrix‐vector product and the diagonal matrix inversion. Another modification consists of a parallel linear equation solver, which is separately constructed based on the multi‐thread and CPU‐GPU heterogeneous platforms with OpenMP and CUDA, respectively. The simulation results from several numerical examples using the modified DDA programme demonstrate that the Jacobi‐PCG is a better iterative method for large‐scale engineering computing and that adoptive parallel strategies can greatly enhance computational efficiency. Copyright © 2015 John Wiley & Sons, Ltd.     

Elastic fields in two joined transversely isotropic media of infinite extent as a result of rectangular loading

This paper presents the closed‐form solutions for the elastic fields in two bonded rocks induced by rectangular loadings. Each of the two bonded rocks behaves as a transversely isotropic linear elastic solid of semi‐infinite extent. They are completely bonded together at a horizontal surface. The rectangular loadings are body forces along either vertical or horizontal directions and are uniformly applied on a rectangular area. The rectangular area is embedded in the two bonded rocks and is parallel to the horizontal interface. The classical integral transforms are used in the solution formulation, and the elastic solutions are expressed in the forms of elementary harmonic functions for the rectangular loadings. The stresses and displacements in the rocks induced by both the horizontal and vertical body forces are also presented. The numerical results illustrate the important effect of the anisotropic bimaterial properties on the stress and displacement fields. The solutions can be easily implemented for numerical calculations and applied to problems encountered in rock mechanics and engineering. Copyright © 2011 John Wiley & Sons, Ltd.     

Abrupt behaviours of streamflow and sediment load variations of the Yangtze River basin,China

Monthly sediment load and streamflow series spanning 1963–2004 from four hydrological stations situation in the main stem of the Yangtze River, China, are analysed using scanning t‐test and the simple two‐phase linear regression scheme. Results indicate significant changes in the sediment load and streamflow from the upper reach to the lower reach of the Yangtze River. Relatively consistent positive coherency relations can be detected between streamflow and sediment load in the upper reach and negative coherency in the middle and lower reaches. Interestingly, negative coherency is found mainly for larger time scales. Changes in sediment load are the result mainly of human influence; specifically, the construction of water reservoirs may be the major cause of negative coherency. Accentuating the human influence from the upper to the lower reach results in inconsistent correlations between sediment load and streamflow. Decreasing sediment load being observed in recent years has the potential to alter the topographical properties of the river channel and the consequent development and recession of the Yangtze Delta. Results of this study are of practical significance for river channel management and evaluation of the influence of human activities on the hydrological regimes of large rivers. Copyright © 2012 John Wiley & Sons, Ltd.     

Preconditioned IDR(s) iterative solver for non‐symmetric linear system associated with FEM analysis of shallow foundation

Non‐associated flow rule is essential when the popular Mohr–Coulomb model is used to model nonlinear behavior of soil. The global tangent stiffness matrix in nonlinear finite element analysis becomes non‐symmetric when this non‐associated flow rule is applied. Efficient solution of this large‐scale non‐symmetric linear system is of practical importance. The standard Krylov solver for a non‐symmetric solver is Bi‐CGSTAB. The Induced Dimension Reduction [IDR(s)] solver was proposed in the scientific computing literature relatively recently. Numerical studies of a drained strip footing problem on homogenous soil layer show that IDR(s = 6) is more efficient than Bi‐CGSTAB when the preconditioner is the incomplete factorization with zero fill‐in of global stiffness matrix K_ep (ILU(0)‐K_ep). Iteration time is reduced by 40% by using IDR(s = 6) with ILU(0)‐K_ep. To further reduce computational cost, the global stiffness matrix K_ep is divided into two parts. The first part is the linear elastic stiffness matrix K_e, which is formed only once at the beginning of solution step. The second part is a low‐rank matrix Δ, which is re‐formed at each Newton–Raphson iteration. Numerical studies show that IDR(s = 6) with this ILU(0)‐K_e preconditioner is more time effective than IDR(s = 6) with ILU(0)‐K_ep when the percentage of yielded Gauss points in the mesh is less than 15%. The total computation time is reduced by 60% when all the recommended optimizing methods are used. Copyright © 2013 John Wiley & Sons, Ltd.     

An application of Karmarkar's interior-point linear programming algorithm for multi-reservoir operations optimization

Optimization of multi-reservoir systems operations is typically a very large scale optimization problem. The following are the three types of optimization problems solved using linear programming (LP): (i) deterministic optimization for multiple periods involving fine stage intervals, for example, from an hour to a week (ii) implicit stochastic optimization using multiple years of inflow data, and (iii) explicit stochastic optimization using probability distributions of inflow data. Until recently, the revised simplex method has been the most efficient solution method available for solving large scale LP problems. In this paper, we show that an implementation of the Karmarkar's interior-point LP algorithm with a newly developed stopping criterion solves optimization problems of large multi-reservoir operations more efficiently than the simplex method. For example, using a Micro VAX II minicomputer, a 40 year, monthly stage, two-reservoir system optimization problem is solved 7.8 times faster than the advanced simplex code in MINOS 5.0. The advantage of this method is expected to be greater as the size of the problem grows from two reservoirs to multiples of reservoirs. This paper presents the details of the implementation and testing and in addition, some other features of the Karmarkar's algorithm which makes it a valuable optimization tool are illuminated.     

A formulation for diagnostic anomaly models

It is shown that the anomaly model governing the difference between two integrations of the equations of motion can be written in a form that is formally linear in the anomalies, by choosing the basic state to be the mean of the results of the two integrations.     

Asymptotic Rayleigh instantaneous unit hydrograph

The instantaneous unit hydrograph for a channel network under general linear routing and conditioned on the network magnitude,N, tends asymptotically, asN grows large, to a Rayleigh probability density function. This behavior is identical to that of the width function of the network, and is proven under the assumption that the network link configuration is topologically random and the link hydraulic and geometric properties are independent and identically distributed random variables. The asymptotic distribution depends only on a scale factor, , where is a mean link wave travel time.