History matching of facies distribution with the EnKF and level set parameterization

In this work, we develop a methodology to combine the Ensemble Kalman filter (EnKF) and the level set parameterization for history matching of facies distribution. With given prior knowledge about the facies of the reservoir geology, initial realizations are generated by commonly used software as the prior guesses of the unknown field. Furthermore, level set functions are used to reparameterize these initial realizations. In the reparameterization process, a representing node system is set up, on which the values of level set functions are assigned using Gaussian random numbers. The mean and the standard deviation of the Gaussian random numbers are designed according to the facies proportion, and the sign of the random numbers depends on the facies type at the representing nodes. The values of the level set functions at the other grid nodes are obtained by linear interpolation. The level set functions on the representing nodes are the model parameters of the EnKF state vector and are updated in the data assimilation process. On the basis of our numerical examples for two-dimensional reservoirs with two or three facies, the proposed method is demonstrated to be able to capture the main features of the reference facies distributions.     

一个新的海边光学湍流外尺度和C_n~2的廓线模式

2016年12月13日至2017年1月2日期间,在茂名博贺海洋气象科学实验基地,采用自行研制的湍流气象探空仪,获取了30份海边温、湿、压、风速、风向和C_n~2等探空数据.基于HMNSP99外尺度模式,利用海边的探空数据拟合得到一个茂名大气光学湍流外尺度经验公式.同时对实验测得的高空湍流廓线数据进行统计平均,然后基于Hufnagel-Valley模式拟合得到符合海边湍流廓线规律的统计平均模式(C_n~2sea model).根据Tatarski高空湍流参数化方案,将用茂名外尺度公式估算的C_n~2分别与探空测量的C_n~2以及用其他外尺度模式估算的C_n~2进行了比较.对其进行统计性分析发现,利用新拟合的茂名外尺度公式、HMNSP99,Dewan以及(Coulman等外尺度模式计算的log10C_n~2)与实测值的整体相关系数分别为0.924,0.848,0.763和0.651,在变化趋势和量级上都表现出较好的一致性;以上四种外尺度模式估算结果的误差都很小,其整体平均绝对误差和平均相对误差分别为0.514和2.963%,0.627和3.612%,0.943和5.439%,0.766和4.417%,新拟合的外尺度模式的误差最小.进一步验证了新的海边外尺度和C_n~2廓线模式的可靠性和有效性,此外还发现高空大气光学湍流的发生与风切变和温度梯度具有十分密切的关系,为光电工程在海边场景应用所需的大气光学湍流廓线模式提供支持.     

基于物理的变形曲线在结构形状优化中的应用

用各向异性模型定义NURBS曲线的变形能,基于所提出的变形能模型,用有限元法对NURBS曲线表达的设计边界进行模态计算,然后,将设计边界用模态向量的线性组合参数化表示。将这种基于边界特征向量的几何形状表示方法应用于优化参数定义．提出了一种适用于结构形状优化的自适应几何精化方法,它有效地将户型有限元分析、优化方法和设计边界的形状表达集成在一起。     

半干旱草原下垫面动量和感热总体输送系数参数化关系研究

本文利用锡林郭勒草原2008年春季近地层涡旋相关系统和铁塔的风、 温平均梯度观测资料, 分析了总体输送系数随梯度Richardson数的变化特征, 建立了动量总体输送系数随大气稳定度、近地层风速以及感热总体输送系数随大气稳定度和近地层气温的关系. 中性条件下, 半干旱草原植被下垫面动量总体输送系数与近地层大气动力状态之间存在明显的相互作用, 总体输送系数与近地面层风速之间满足二次曲线拟合关系; 风速较小时, 大气动力特征对地表粗糙度长度的改变不是很明显, 动量总体输送系数随气流增强而增大; 而当风速较大时, 强风速会使植被高度发生改变, 动量总体输送系数随气流增强而减小. 另外, 感热总体输送系数与近地层气温之间也存在二次曲线关系. 动量总体输送系数与近地层风速之间的关系、感热总体输送系数与近地层气温之间关系的建立为总体输送系数参数化提供了重要途径, 同时该方案避免了对动力学粗糙度长度和热力学粗糙度长度的求解. 关键词： 总体输送系数 参数化 湍流通量 相似性函数     

基于实验设计方法的高超声速飞机前缘型线优化分析

为探索前缘线变化对吸气式高超声速飞机气动性能的影响,基于一种旁侧进气布局翼身融合体构型,在飞行马赫数6,攻角4°和高度26 km的巡航飞行条件下,结合运用增量修正参数化设计方法、均匀实验设计方法和计算流体力学模拟,分析了飞行器前缘型线与其升阻力系数及纵向压心等性能参数间的关系.计算结果表明,前缘线形状对飞行器升阻力系数明显高于其对纵向压心影响,设计空间范围内升力系数变化约21.3%,阻力系数变化约31.8%,升阻比变化范围约10.63%,但相对压心变化范围仅为3.87%.在此基础上,通过对典型构型物面压力分布进行分析,发现前缘线形状适当弯曲可利用飞行器下表面侧壁压缩产生的高压气流,利用二者的耦合效应使飞行器获得额外的升力增量.      

五氧化二氮(N2O5)非均相摄取系数的定量和参数化

颗粒物污染是我国超大城市地区面临的主要污染问题,大气化学受到非均相反应的强烈影响。N₂O₅非均相反应过程的定量描述对厘清我国超大城市地区夜间大气氧化性、区域NO_x收支、臭氧生成等问题具有重要意义。本文从反应机理、摄取系数测量技术和不同气溶胶基体测量结果等方面系统总结了N₂O₅非均相反应动力学研究进展。N₂O₅非均相摄取是典型的反应性摄取过程,是气溶胶流动管测量系统的理想应用对象。实验室动力学研究从硫酸盐等模型颗粒物出发,结合外场观测获取的颗粒物信息,模拟使用的颗粒物越来越接近实际大气条件。N₂O₅在不同颗粒物表面摄取系数在0.001~0.2之间,其变化范围超过两个数量级,主要受到环境温度、湿度、颗粒物组分中NO₃^-、Cl^-、SO₄^2-、含水量、有机物、颗粒物混合状态和相态等因素影响;颗粒物中H₂O、Cl^-、SO₄^2-的浓度升高会提高非均相摄取速率,而NO₃^-、有机物则会对摄取过程产生抑制作用。欧美地区外场观测研究表明这些影响因素并不是线性叠加的,而是存在非常复杂的相互作用,基于实验室测量结果开发的参数化公式不能很好地拟合观测结果。由于N₂O₅和颗粒物的高值区共存于污染城市地区,下一步研究需要在包括我国超大城市地区在内的典型城市大气条件下开展更多的外场实测研究和建立更准确的参数化公式。     

Dark Parameterizations,Equivalent Partner Fields and Integrable Systems

After introducing dark parameters into the traditional physical models, some types of new phenomena may be found. An important difficult problem is how to directly observe this kind of physical phenomena. An alternative treatment is to introduce equivalent multiple partner fields. If use this ideal to integrable systems, one may obtain infinitely many new coupled integrable systemsconstituted by the original usual field and partner fields. The idea is illustrated via the celebrate KdV equation. From the procedure, some byproducts can be obtained: A new method to find exact solutions of some types of coupled nonlinear physical problems, say, the perturbation KdV systems, is provided; Some new localized modes such as the staggered modes can be found and some new interaction phenomena like the ghost interaction are discovered.     

MCMC Estimation of Restricted Covariance Matrices

This article is motivated by the difficulty of applying standard simulation techniques when identification constraints or theoretical considerations induce covariance restrictions in multivariate models. To deal with this difficulty, we build upon a decomposition of positive definite matrices and show that it leads to straightforward Markov chain Monte Carlo samplers for restricted covariance matrices. We introduce the approach by reviewing results for multivariate Gaussian models without restrictions, where standard conjugate priors on the elements of the decomposition induce the usual Wishart distribution on the precision matrix and vice versa. The unrestricted case provides guidance for constructing efficient Metropolis–Hastings and accept-reject Metropolis–Hastings samplers in more complex settings, and we describe in detail how simulation can be performed under several important constraints. The proposed approach is illustrated in a simulation study and two applications in economics. Supplemental materials for this article (appendixes, data, and computer code) are available online.     

Experimental mixing parameterization due to multiphase fluid–structure interactions

Experimental estimates of the rate at which energy is extracted from the barotropic tide at corrugated sloping topography are presented. To this end, a series of experimental simulations of the process of generation of internal tides (i.e., internal waves of the tidal frequency) over a corrugated slope in stratified fluid were performed. An oceanic interior is modeled as four-layer stratified fluid. The main focus in these studies was to obtain the relation for the potential energy available for mixing as a function of a slope of a corrugated slope. The available potential energy available for partial mixing to the topographic slope was compared with the available potential energy sufficient for complete mixing of the four layers. The experimental data were compared with the analytic results of a linear theory and found in agreement with a recent theoretically predicted scenario showing that the dominant contribution of the energy distribution in internal wave field resides in waves of the lowest allowed frequency.