大地电磁自适应正则化反演算法

针对大地电磁正则化反演中正则化因子的选取困难问题提出了自适应正则化反演算法（Adaptive Regularized Inversion Algorithm, ARIA）. 在该算法中, ①提出了一种新的数据方差处理方法：数据方差规范化,使得数据方差的大小只对数据的拟合发生影响,不对数据目标函数和模型约束目标函数的权重产生影响,从而减少了正则化因子取值的影响因素;②提出了粗糙度核矩阵的概念,并给出了由基本结构插值基函数计算粗糙度核矩阵的公式,使得模型目标函数的构建更为简便、直接;③根据数据目标函数、模型约束目标函数和正则化因子之间的关系,提出了两种正则化因子自适应调节方法. 本文详细阐述了最平缓模型约束下的大地电磁一维连续介质反演的ARIA实现,以几个算例的分析比较来说明ARIA的有效性.     

多尺度阵列嵌套组合反演宾川气枪源区横波速度结构

本文基于多尺度阵列嵌套组合的方式,利用频率-贝塞尔变换法（Frequency-Bessel,F-J方法）提取背景噪声面波频散信息,通过多个阵列融合的频散曲线反演得到宾川气枪发射台周边不同深度的横波速度结构.结果显示：浅层一阶面波频散信息的加入,使得基阶反演结果更加收敛,反演深度加深到8 km;深度在8 km以下的结构的研究利用多尺度阵列（密集台阵-宾川气枪台网-云南区域地震台网）嵌套组合的方式,面波基阶低频信息从0.55 Hz拓宽到0.008 Hz,使横波速度结构的反演深度显著增加,同时对反演过程提供约束,使得70 km深度以上的横波速度更收敛.由此本文所得的横波速度结构为该区地下结构的探测提供基础,多尺度阵列嵌套组合频散谱的研究方式也为以后区域结构的研究提供一种新的方法和思路.

     

Reducing the dimensionality of compositional data sets

The high-dimensionality of many compositional data sets has caused geologists to look for insights into the observed patterns of variability through two dimension-reducing procedures: (i)the selection of a few subcompositions for particular study, and (ii)principal component analysis. After a brief critical review of the unsatisfactory state of current statistical methodology for these two procedures, this paper takes as a starting point for the resolution of persisting difficulties a recent approach to principal component analysis through a new definition of the covariance structure of a composition. This approach is first applied for expository purposes to a small illustrative compositional data set and then to a number of larger published geochemical data sets. The new approach then leads naturally to a method of measuring the extent to which a subcomposition retains the pattern of variability of the whole composition and so provides a criterion for the selection of suitable subcompositions. Such a selection process is illustrated by application to geochemical data sets.     

Evaluation of grain shear stresses required to initiate movement of particles in natural rivers

Shear stresses were evaluated at different sites on two rivers. The first (the Rulles) is characterized by a pebbly bedload and a meandering bed with riffles and pools. The second (the Rouge Eau) has mainly a sandy rippled bed where meandering is well developed but also flat gravelly sectors without meandering system. Shear stresses calculated from friction velocities (τ*) using a redefined y₁ roughness height parameter were compared with total shear stresses calculated from the energy grade line and the hydraulic radius (τ), Divergence between these shear stresses seems to increase in the presence of bedforms and large-scale irregularities of the channel. The τ*/τ ratio is close to 0·5 in the gravelly sector of the Rouge Eau and reaches 0·65 in the riffles of the Rulles (generally located at the inflexion point of the meanders), while it is less than 0·3 in the pools of the same river (located in the loops) and only 0·2 in the sandy rippled sector of the Rouge Eau. Grain and bedform shear stresses were evaluated at these same sites by different methods. The grain shear stress (τ') represents on average 30 per cent of the total shear stress in the riffles of the Rulles and the gravelly sector of the Rouge Eau, but less than 15 per cent in the pools in the Rulles and the sandy sectors of the Rouge Eau. However, it emerges from experiments conducted with marked pebbles and in situ observations of erosion and transport of sandy and gravelly particles, that the grain shear stresses are underestimated and cannot explain the movements and modifications actually observed. Conversely, shear stresses calculated from friction velocities at the sites where erosion actually occurred (or failed to occur despite very high velocities) provide a better explanation of the observed movements.     

大黄鱼(Larimichthys crocea)新品种“东海1号”体长相关的DArT标记筛选

大黄鱼(Larimichthys crocea)是我国重要的海产经济鱼类。因多年未加选育的累代养殖,养殖大黄鱼生长缓慢,性早熟、免疫力低下,亟需对养殖大黄鱼进行遗传改良。多样性芯片技术(Diversity arrays technology,DAr T)具有高通量和低成本的显著特点,不需要明确物种的基因组DNA序列信息,因而广泛应用于动植物遗传图谱制作和遗传多样性分析。本研究旨在采用DAr T技术鉴定与"东海1号"大黄鱼体长相关分子标记。首先随机测得"东海1号"大黄鱼199个个体体长,均值为13.45cm,符合正态分布(P0.05),"极端大群体"和"极端小群体"之间差异显著(P0.01)。然后采用7组限制性内切酶组合(Pst I/Alu I、Pst I/Ban II、Pst I/Bsp1286I、Pst I/Bst NI、Pst I/Hae III、Pst I/Rsa I、Pst I/Taq I)分别进行酶切,获得基因组代表性DNA片段并用于文库构建。根据多态性克隆数和多态性率确定Pst I/Rsa I酶切组合为最优降低基因组复杂度方法。之后从Pst I/Rsa I基因组代表性DNA片段文库中扩增获得3360个片段,以此为多样性芯片探针进行芯片点制,杂交筛选获得18个大黄鱼体长相关DAr T候选标记。经过新群体样本再次验证,仍有8个DAr T标记与体长相关。本研究有助于选育生长性状优良的大黄鱼群体。     

风浪状态对海面拖曳系数的影响

通过对四种参数化方法的比较,认为风海雷诺数RB较风速、波龄和波陡能更好地描述拖曳系数CDN。利用分组平均法对CDN与RB关系进行处理,得到最佳的CDN参数化方案。利用COARE3算法测试了四种依赖海况的粗糙长度,与实测结果进行了比较,结果表明CDN与RB的关系式更真实地反映了海-气界面动量交换过程。     

一个标量场的客观分析方法和资料点空间结构的设计

这是一个根据Ｇａｕｓｓ－Ｍａｒｋｏｖ原理建立的标量场客观分析方法，并在资料的模拟中得到了应用。在这个方法的推导过程中还给出了计算ｒ．ｍ．ｓ误差的公式。通过试错调整格点的空间结构，可以设计出有效的空间取样阵列。本文使用这个方法对温度场、气压场和高度场分别进行了分析．试验结果表明这个方法简单、方便，分析结果令人满意。     

Scalar Concentration Profiles in the Canopy and Roughness Sublayer

The mean flow and scalar concentration profiles within and above a tall canopy are well known to violate the standard boundary-layer flux-gradient relationships. We present a theory for the scalar concentration profile that is comprised of a canopy exchange model coupled to a modified surface-layer model. The coupling between the two components and the modifications to the surface-layer profiles are formulated through the mixing-layer analogy for the flow at canopy top. This analogy provides an additional length scale—the vorticity thickness—upon which the profiles depend and a set of criteria that allows a reduction in the empiricism associated with earlier forms in the literature. Predictions of the mean scalar concentration profiles are shown to match observations over a wide range of diabatic stabilities for both potential temperature and water vapour.     

Aerodynamic Parameters of Regular Arrays of Rectangular Blocks with Various Geometries

The aerodynamic effects of various configurations of an urban array were investigated in a wind-tunnel experiment. Three aerodynamic parameters characterising arrays—the drag coefficient (C _d ), roughness length (z _o) and displacement height (d)—are used for analysis. C _d is based on the direct measurement of the total surface shear using a floating element, and the other two parameters are estimated by logarithmic fitting of the measured wind profile and predetermined total drag force. The configurations of 63 arrays used for measurement were designed to estimate the effects of layout, wind direction and the height variability of the blocks on these parameters for various roughness packing densities. The results are summarised as follows: (1) The estimated C _d and z _o of the staggered arrays peak against the plan area index (λ _p ) and frontal area index (λ _f ), in contrast with values for the square arrays, which are less sensitive to λ _p and λ _f . In addition, the square arrays with a wind direction of 45° have a considerably larger C _d , and the wind direction increases z _o/H by up to a factor of 2. (2) The effect of the non-uniformity of roughness height on z _o is more remarkable when λ _f exceeds 20%, and the discrepancy in z _o is particularly remarkable and exceeds 200%. (3) The effect of the layout of tall blocks on C _d is stronger than that of short blocks. These results indicate that the effects of both wind direction and the non-uniformity of the heights of buildings on urban aerodynamic parameters vary greatly with λ _p and λ _f ; hence, these effects should be taken into account by considering the roughness packing density.     

Note on a Parametric Relation for Separating Flow over a Rough Hill

A parametrization method used to account for the effects of flow separation and wall roughness on the lower boundary condition for turbulent boundary layers is investigated against direct numerical simulation and laser Doppler anemometry data. The numerical simulation represents flow over a smooth, flat surface with a prescribed external adverse pressure gradient. The water-channel experiments cover flow over smooth and rough hills for two specified Reynolds numbers. Global optimization algorithms based on four different direct search methods are used to assess the parametrization function, C, in terms of local mean velocity profiles and the parametrization parameters u _* (friction velocity), ∂ _x p (local pressure gradient), z ₀ (effective roughness) and d (zero-plane displacement). The study investigates regions of attached and reversed flows, and forty-two velocity profiles are compared with the proposed expression for the function C, including two profiles that satisfy the solution of Stratford.