流速逆变张量隐式求解方法及其在航道港池流场计算中的应用

运用高分辨率的边界适应网格进行流体动力学数值计算时,如何提高计算稳定性和减少计算量成为数值求解的关键性问题.在非正交的边界适应坐标系中,每个动量方程中同时出现了两个交叉方向的水位偏导数项,给隐式求解带来困难,而显式格式下的时间步长由于受与空间步长有关的Courant-Friedrichs-Lewy条件限制,计算量成倍增加.本文从广义曲线坐标系下浅海动力学方程组出发,导出了流速的逆变张量所满足的动量方程组,使方程中的水位偏导数项变成了沿某一协变基向量方向占优的形式,方便地采用了交替方向隐式差分格式,从而提高了计算稳定性并减小了计算量.本文通过对澳门海域航道和港池中流场的计算,证实了该模式是一种进行高分辩率数值计算的有效方法.     

面向目标自适应三维大地电磁正演模拟

本文将面向目标的自适应算法应用于三维大地电磁数值模拟.使用基于非结构网格的矢量有限单元法对起伏地表大地电磁正演模拟问题进行求解.使用利用垂向电流密度在物性界面上的连续性对后验误差进行估算的算法指导网格优化.由于全局自适应算法针对观测点优化网格的能力较差,本文通过求解正演问题的对偶问题计算后验误差的加权系数,并对相关加权系数进行改进,从而实现了面向目标的自适应算法.与传统基于结构化网格的电磁正演算法相比,采用非结构网格能够更好地拟合起伏地表和地下不规则异常体.由于使用了面向目标的自适应算法,本文能够使用更少的网格达到较高的计算精度.通过对比本文模拟结果与半空间响应和全局自适应算法计算结果,并通过对比使用改进前和改进后加权系数得到的网格剖分结果验证了本文算法的有效性.

     

Comparison of direct statistical and indirect statistical-deterministic frameworks in downscaling river low-flow indices

ABSTRACT

This work explores the ability of two methodologies in downscaling hydrological indices characterizing the low flow regime of three salmon rivers in Eastern Canada: Moisie, Romaine and Ouelle. The selected indices describe four aspects of the low flow regime of these rivers: amplitude, frequency, variability and timing. The first methodology (direct downscaling) ascertains a direct link between large-scale atmospheric variables (the predictors) and low flow indices (the predictands). The second (indirect downscaling) involves downscaling precipitation and air temperature (local climate variables) that are introduced into a hydrological model to simulate flows. Synthetic flow time series are subsequently used to calculate the low flow indices. The statistical models used for downscaling low flow hydrological indices and local climate variables are: Sparse Bayesian Learning and Multiple Linear Regression. The results showed that direct downscaling using Sparse Bayesian Learning surpassed the other approaches with respect to goodness of fit and generalization ability.

Editor D. Koutsoyiannis; Associate editor K. Hamed     

Approximating constant‐Q seismic propagation in the time domain

In this study, we investigate the accuracy of approximating constant‐Q wave propagation by series of Zener or standard linear solid (SLS) mechanisms. Modelling in viscoacoustic and viscoelastic media is implemented in the time domain using the finite‐difference (FD) method. The accuracy of numerical solutions is evaluated by comparison with the analytical solution in homogeneous media. We found that the FD solutions using three SLS relaxation mechanisms as well as a single SLS mechanism, with properly chosen relaxation times, are quite accurate for both weak and strong attenuation. Although the RMS errors of FD simulations using a single relaxation mechanism increase with increasing offset, especially for strong attenuation (Q = 20), the results are still acceptable for practical applications. The synthetic data of the Marmousi‐II model further illustrate that the single SLS mechanism, to model constant Q, is efficient and sufficiently accurate. Moreover, it benefits from less computational costs in computer time and memory.     

基于损伤部位向量法的钢框架损伤识别研究

对损伤部位向量(DLV)法作了简单介绍,并用该方法对钢框架进行了损伤识别和损伤定位。该方法假定结构损伤前后为线性,对结构损伤前后柔度矩阵差进行奇异值分解,将奇异值为零所对应的向量,作为静荷载施加在无损结构的测点位置,则应力为零的单元为可能损伤的单元。对3种不同工况的钢框架进行了振动模态试验,用前3阶模态参数构造框架的柔度矩阵,按照DLV法对其进行了损伤识别,识别结果与已知损伤情况相一致。从测试自由度不完备、噪声和振型质量归一化系数这3个方面对识别效果进行了分析,结果表明:当损伤使结构动力特性有微小改变时,使用该方法不易定位损伤,应结合局部损伤识别方法进行判定;当损伤使结构动力特性有较大改变时,该方法能有效识别损伤的单元。DLV方法概念简单,理论明确,不受结构类型的限制,不需要结构的数学模型和模型缩聚或扩展技术,只需获得结构损伤前后的前几个低阶模态参数,即可识别结构一处或多处损伤,实际应用时可操作性强。     

基于Hilbert-Huang变换的地震瞬时属性提取方法及应用(英文)

希尔伯特—黄变换(Hilbert-Huang transform,HHT)是一种新的适合非平稳和非线性信号的分析方法,由于地震信号一般呈现出非平稳与非线性特性,因此HHT非常适合地震信号的分析。本文首先介绍了HHT中关于经验模态分解(Empirical Mode Decomposition,EMD)的实现过程,在此基础上分析了几种基于EMD获得本征模态函数(Intrinsic Mode Functions,IMF)来计算瞬时频率的算法,其中利用了两个采样间隔瞬时频率的平均来计算瞬时频率,较好地反映了地震信号频率成分随时间变化的特征。将该方法应用于四川东北部某地区海相碳酸盐岩地层三维地震叠后偏移数据处理,提取"三瞬"地震属性,与传统的希尔伯特变换提取的"三瞬"地震属性进行对比,结果表明基于HHT的"三瞬"地震属性结果具有更高的分辨率,IMF2的瞬时相位能够较好地刻画台地边缘生物礁相,IMF2的瞬时频率亦具有较好的分带性。将IMF2的"三瞬"地震属性与钻井等资料结合分析,能够更好地识别沉积相的分布。     

Decomposing productivity and efficiency changes in the Alaska head and gut factory trawl fleet

Fishing fleets are subject to numerous factors that affect economic performance, making identification and attribution of such impacts difficult. This paper separately identifies the effects of changing input and output prices, fishery management, and quota allocations on total factor productivity using a Lowe Index. Indices account for technical change and decompose productivity estimates into its technical, environmental, and scale-mix components. This results in measures that reflect shifts in the production frontier, and movements by vessels toward and around the frontier, to capture economies of scale and mix after a policy shift to a catch share program that includes fishing cooperatives and a limited access fishery. The difference between cooperative and limited access vessels is exploited to compare the changes in economic performance between the groups after the introduction of the shift to catch shares and cooperative management, which allowed the vessels to improve the timing and coordination across multi-species fisheries and to decrease incidental catch of quota-limited bycatch species that had closed the target fisheries prematurely in the past. Results indicate that total factor productivity increased significantly after the move to a catch share program, largely due to increases in technical change that shifted out the production frontier of the fishery.     

I Zwicky 1: Decomposition and dynamics of the nearby QSO host

We address current theories of research on morphology, environment and evolution of active galaxies in the form of a detailed case study of the nearby QSO host I Zw 1. This study is based on sub-kpc resolution ¹²CO(1-0) observations with the BIMA mm-interferometer and on near-infrared imaging with ISAAC at the VLT. The ¹²CO(1-0) maps reveal a circumnuclear molecular gas ring with a radius of 0.9 kpc. The imaging data in the J-band are analyzed with respect to a disk-bulge decomposition. Together with a model of the gas rotation curve, the radial profile of the dynamical J-band mass-to-light ratio (M/L) is discussed. The J-band images give new evidence for an interaction between I Zw 1 and its nearby companion. First results from a sequence of N-body simulations for an extensive test of the parameter space of the interaction are presented. This revised version was published online in August 2006 with corrections to the Cover Date.     

The Q value variations in the preparing pro－cess of rock rupture

ＴｈｅＱｖａｌｕｅｖａｒｉａｔｉｏｎｓｉｎｔｈｅｐｒｅｐａｒｉｎｇｐｒｏ┐ｃｅｓｏｆｒｏｃｋｒｕｐｔｕｒｅＢＩＮＷＡＮＧ１）（王彬），ＺＨＡＯ－ＹＯＮＧＸＵ１）（许昭永），ＪＩＮ－ＭＩＮＧＺＨＡＯ１）（赵晋明）ＹＩ－ＬＩＨＵ２）（胡毅力）１）Ｓｅｉ...     

注水诱发地震的震源应力场和尾波衰减特征

利用区域性台网记录测定了任邱油田注水诱发地震的震源应力场和尾波衰减特征,结果表明地壳上层有很大水平应力,断层面解受注水井附近局部原生断层影响与华北现代应力场的优势破裂方向和错动方式有差别。地震应力降分布范围为0.2—3.0巴。平均尾波Q_0值为75,属于地壳上部沉积岩层的正常情况,发现在时间上表现出与注水量增减有关的波动。在空间上表现出震群密集区Q_0值较外围区偏低的特点