回归估计精度的局部影响

从模型微小扰动的角度研究数据对回归估计精度的局部影响。首先提出刻划这种影响的４个统计量，并求出它们在３种不同扰动方式下的最大局部变化方向，其最大绝对分量对应于强影响点。最后给出一个数值例子来说明这种强影响点的诊断。     

联合多点地质统计学与序贯高斯模拟的随机反演方法

岩相和储层物性参数是油藏表征的重要参数,地震反演是储层表征和油气藏勘探开发的重要手段.随机地震反演通常基于地质统计学理论,能够对不同类型的信息源进行综合,建立具有较高分辨率的储层模型,因而得到广泛关注.其中,概率扰动方法是一种高效的迭代随机反演策略,它能综合考虑多种约束信息,且只需要较少的迭代次数即可获得反演结果.在概率扰动的优化反演策略中,本文有效的联合多点地质统计学与序贯高斯模拟,并结合统计岩石物理理论实现随机反演.首先,通过多点地质统计学随机模拟,获得一系列等可能的岩相模型,扰动更新初始岩相模型后利用相控序贯高斯模拟建立多个储层物性参数模型;然后通过统计岩石物理理论,计算相应的弹性参数;最后,正演得到合成地震记录并与实际地震数据对比,通过概率扰动方法进行迭代,直到获得满足给定误差要求的反演结果.利用多点地质统计学,能够更好地表征储层空间特征.相控序贯高斯模拟的应用,能够有效反映不同岩相中储层物性参数的分布.提出的方法可在较少的迭代次数内同时获得具有较高分辨率的岩相和物性参数反演结果,模型测试和实际数据应用验证了方法的可行性和有效性.     

弱二次切变流轴对称涡旋演变的解析研究

运用格林函数方法,给出了弱二次切变流轴对称涡旋演变的解析解.结果表明:当环境流中弱二次切变强度增加(减少)时,扰动动能的最大值增加(减小),扰动动能的变化加速(减速);而扰动涡旋的范围对扰动动能的大小与变化也存在一定的影响.     

台湾岛对气流扰动的影响

本文根据一个三维气流过山的模型,用数值快速傅里叶变换方法,讨论理想东北季风气流经过台湾岛产生的扰动气压场、扰动垂直速度场和水平流场.经与实际东北季风下台湾岛的气候特征进行比较,结果令人满意.     

Conditional nonlinear optimal perturbation: Applications to stability,sensitivity, and predictability

Conditional nonlinear optimal perturbation (CNOP) is a nonlinear generalization of linear singular vector (LSV) and features the largest nonlinear evolution at prediction time for the initial perturbations in a given constraint. It was proposed initially for predicting the limitation of predictability of weather or climate. Then CNOP has been applied to the studies of the problems related to predictability for weather and climate. In this paper, we focus on reviewing the recent advances of CNOP’s applications, which involves the ones of CNOP in problems of ENSO amplitude asymmetry, block onset, and the sensitivity analysis of ecosystem and ocean’s circulations, etc. Especially, CNOP has been primarily used to construct the initial perturbation fields of ensemble forecasting, and to determine the sensitive area of target observation for precipitations. These works extend CNOP’s applications to investigating the nonlinear dynamical behaviors of atmospheric or oceanic systems, even a coupled system, and studying the problem of the transition between the equilibrium states. These contributions not only attack the particular physical problems, but also show the superiority of CNOP to LSV in revealing the effect of nonlinear physical processes. Consequently, CNOP represents the optimal precursors for a weather or climate event; in predictability studies, CNOP stands for the initial error that has the largest negative effect on prediction; and in sensitivity analysis, CNOP is the most unstable (sensitive) mode. In multi-equilibrium state regime, CNOP is the initial perturbation that induces the transition between equilibriums most probably. Furthermore, CNOP has been used to construct ensemble perturbation fields in ensemble forecast studies and to identify sensitive area of target observation. CNOP theory has become more and more substantial. It is expected that CNOP also serves to improve the predictability of the realistic predictions for weather and climate events plays an increasingly important role in exploring the nonlinear dynamics of atmospheric, oceanic and coupled atmosphere-ocean system. Supported by National Basic Research Program of China (Grant Nos. 2006CB403606, 2007CB411800), National Natural Science Foundation of China (Grant Nos. 40830955, 40675030, 40505013), Institute of Atmospheric Physics, Chinese Academy of Sciences (Grant No. IAP07202), and LASG State Key Laboratory Special Fund     

ENSEMBLE PREDICTION EXPERIMENTS OF TYPHOON TRACK BASED ON THE STOCHASTIC TOTAL TENDENCY PERTURBATION

The GRAPES-TCM is used to make ensemble prediction experiments for typhoon track. Three kinds of ensemble schemes are designed for the experiments. A total of 109 experiments are made for the nine typhoons in 2011 and the integral time is 72 h. The experiment results are shown as follows. In the three ensemble schemes, on the whole, scheme 1 has the best track prediction. Its average absolute track error and overall deviations of typhoon moving speed and moving direction are all the smallest in the three schemes. For both scheme 1 and scheme 2, they are all smaller than those of their control predictions. Both of their ensemble predictions show superiority to their deterministic predictions. Overall, compared with the observations, the typhoon moving directions of the three schemes mainly skew to the right, and in the late integration they mainly tend to be relatively slow. In the three schemes, the track dispersion of scheme 1 is the largest and that of scheme 3 the smallest. In scheme 1 it is much larger than in schemes 2 and 3. The difference of dispersion between scheme 2 and scheme 3 is small. The track dispersions of the three schemes are all much smaller than their rational dispersions. Compared with the eight domestic and overseas operational numerical weather prediction (NWP) models, scheme 1 has better predictions than the other seven operational models except ECMWF NWP model. Scheme 1 has the value of operational application.     

集合预报研究及应用进展

从初值扰动的生成方法、数值预报模式的建立、输出产品的提炼应用等方面简要概括了集合预报的最新进展。指出，运用BGM法生成初值扰动效果较好，多模式集合预报正成为天气预报的发展趋势；集合概率预报是集合预报的最佳结果，但发布还需要做大量的普及推广和培训工作；集合预报正在被越来越多的国家投入准业务运行和业务运行．有可能成为新的主流预报方式。     

中国北方春季沙尘暴频数与北半球500hPa高度场的SVD分析

选取1957-2000年中国北方地区春季沙尘暴发生次数资料和北半球500hPa秋、冬、春季的平均高度场资料,对沙尘暴和高度场作SVD分析。结果表明,我国北方春季沙尘暴次数各地具有比较一致的变化趋势,表明可能与大尺度气候背景的变化有联系。变化的最敏感区域为内蒙古中西部地区、新疆西部、青海西部和东北地区。我国北方春季沙尘暴发生次数与同期及前期500hPa高度场有较好的相关关系,前期环流形势对春季沙尘暴频数有一定的指示和预测意义,冬季环流场尤其具有预报指示意义,因为前一年冬季北大西洋涛动对我国春季北方沙尘暴发生次数有影响。     

A semianalytical approach for solving first-order perturbation equations of dissolution-timescale reactive infiltration instability problems in fluid-saturated rocks

This paper presents a semianalytical approach for solving first-order perturbation (FOP) equations, which are used to describe dissolution-timescale reactive infiltration instability (RII) problems in fluid-saturated rocks. The proposed approach contains two parts because the chemical dissolution reaction divides the whole problem domain into two subdomains. In the first part, the interface-condition substitution strategy is used to derive the analytical expressions of purely mathematical solutions for the FOP equations in the upstream subdomain, where the dissolution chemical reaction is ceased and the FOP equations are weakly coupled. In the second part, the finite element method (FEM) is used to derive the analytical expressions of numerical solutions for the FOP equations in the downstream subdomain, where the dissolution chemical reaction needs to be considered and the FOP equations are strongly coupled so that it is impossible to derive purely mathematical solutions for them. Particular attention is paid to the development of the element-by-element forward marching strategy, which is associated with the use of the FEM for solving this new kind of scientific problem. The related analytical results demonstrated that (1) both the dynamic characteristic of a reactive infiltration system and the dimensionless wavenumber can have pronounced influences on the distribution of the FOP dimensionless acid concentration within the entire domain of the dissolution-timescale RII problems in fluid-saturated rocks and (2) the FOP dimensionless acid concentration distribution exhibits two significantly different patterns in the upstream and downstream subdomains of the dissolution-timescale RII system.     

Secular resonances in the asteroid belt: Theoretical perturbation approach and the problem of their location

In this paper a theoretical perturbation approach to the problem of the dynamics in secular resonance is exposed. This approach avoids any expansion of the main term of the Hamiltonian (linear term in the masses) with respect to the eccentricity or the inclination of the asteroid, in order to achieve results valid for any value of these variables. Moreover suitable action-angle variables are introduced to take properly into account the dynamics related to the motion of the argument of perihelion of the asteroid, which is relevant at high inclination. A class of secular resonances wider than that usually considered is found. An explicit computation of the location of the main secular resonances, estimating also the contribution of the quadratic term in the masses by means of classical series expansion, is reported in the last sections. The accuracy of computations obtained by series expansion is discussed in the paper.