任意曲线曲面上重磁位场转换的B样条函数法

提出用B样条函数求解曲线、曲面上重磁位场的向上延拓,水平、垂向导数计算,磁异常分量互换的方法。该方法的特点是:原理简明,程序通用性强,计算精度高。     

Estimation of Water Particle Velocities by Empirical Transfer Function

-In previous and this studies it appears that the linear and nonlinear wave theory can notaccurately and easily predict the water particle velocities.Therefore,different from the theoretical consider-ations,in this study we have attempted to determine the transfer function empirically.Laboratory experi-ments were performed under various wave conditions.The empirical formulas of the transfer function ofthe wave height,angular frequency and water particle velocity were obtained on the basis of these test databy dimensional analysis and regression analysis.In intermediate and deep water depth conditions,thetransfer function was only a function of a nondimensional parameter which is composed of the angular fre-quency,the depth of the velocity gauge under the still water level,water depth and the acceleration of grav-ity.Finally,the empirical formulas were compared with experimental data and observational data formpresent and Cavaleri's(1978)studies.The empirical formulas were found to be in sufficient correl     

长江口三维潮流数值计算及动力分析

采用三维动力学数学模型对长江口潮流场进行了计算，对长江口潮流流态在一个潮周期内逐时进行了分析，并对长江口南支盐水倒灌从动力学角度进行了初步探讨。计算结果表明：文中模型较好地重演长江口潮流场，可以用于大型复杂河口工程实际中的潮流场的计算和分析。     

中国西部地区基于经验模型的土壤水分反演

地表土壤水分含量的时空分布信息是十分重要的,常常作为水文模型、气候模型、生态模型的输入参数,同时,也是干旱预报、农作物估产等工作的重要指标。被动微波遥感是监测土壤含水量最有效的手段之一。相比红外与可见光,它具有波长长,穿透能力强的优势。相比主动微波雷达,被动微波辐射计具有监测面积大、周期短,受粗糙度影响小,对土壤水分更为敏感,算法更为成熟的优势。目前,已研究出许多反演土壤水分的方法.本课题的主要内容是借助AMSR-E土壤水分影像数据、MODIS归一化植被指数(NDVI)影像数据和MODIS分类影像数据,利用ENVI软件进行遥感图像数据处理,运用统计分析方法建立NDVI与土壤水分的经验模型,研究中国西部地区稀疏植被覆盖区土壤水分的反演。     

正交距离圆曲线拟合方法

在分析圆曲线拟合准则基础上,提出采用正交最小二乘法拟合圆曲线。该方法以圆曲线正交距离残差平方和极小为准则。基于间接平差原理详细推导了相关模型和算法,并以实际算例说明了该方法的效果。     

水平井含水率上升影响因素

利用塔里木油田塔中4油田(TZ4)底水油藏相关的地质、流体数据建立数值模型。在所建模型的基础上,应用数值模拟计算的累积产油、产水和产液量回归俞启泰水驱特征曲线,以反映水平井见水特征的参数b。以参数b为研究对象,采用正交试验的方法研究不同因素对b值的影响,筛选影响水平井见水特征的主要因素,认为原油粘度、油层厚度、非均质性及水平井水平段在油藏中的位置是影响水平井含水上升趋势的主要指标。最后建立主要因素与b值的关系式,结合俞启泰曲线b值图版,提出预测水平井见水规律的公式——图版法(F-b法)。应用F-b法对塔里木水平井含水率进行预测,并与其他相关方法及实际生产数据对比,认为F-b法可作为预测塔里木油田水平井含水率、估算可采储量的一种有效方法。     

经验模式分解的尺度阈值滤波在变形分析中的应用

本文将信号处理领域的经验模式分解算法应用于变形信息的提取中,通过引入阈值函数,建立了基于经验模式分解的尺度阈值滤波模型,采用优化模型确定了经验模式分解的次数。分别通过模拟试验和实测数据与小波阈值法和多项式拟合法进行了比对,分析表明:在低噪声情况下,三种方法都有一定的滤波效果;在高噪声情况下,经验模式分解的尺度阈值滤波法具有与小波阈值去噪法相等的精度,而且瞬时强噪声识别能力更好,优于多项式拟合法。     

Rupture process of large earthquakes in the northern Mexico subduction zone

The Cocos plate subducts beneath North America at the Mexico trench. The northernmost segment of this trench, between the Orozco and Rivera fracture zones, has ruptured in a sequence of five large earthquakes from 1973 to 1985; the Jan. 30, 1973 Colima event (M _s 7.5) at the northern end of the segment near Rivera fracture zone; the Mar. 14, 1979 Petatlan event (M _s 7.6) at the southern end of the segment on the Orozco fracture zone; the Oct. 25, 1981 Playa Azul event (M _s 7.3) in the middle of the Michoacan gap; the Sept. 19, 1985 Michoacan mainshock (M _s 8.1); and the Sept. 21, 1985 Michoacan aftershock (M _s 7.6) that reruptured part of the Petatlan zone. Body wave inversion for the rupture process of these earthquakes finds the best: earthquake depth; focal mechanism; overall source time function; and seismic moment, for each earthquake. In addition, we have determined spatial concentrations of seismic moment release for the Colima earthquake, and the Michoacan mainshock and aftershock. These spatial concentrations of slip are interpreted as asperities; and the resultant asperity distribution for Mexico is compared to other subduction zones. The body wave inversion technique also determines theMoment Tensor Rate Functions; but there is no evidence for statistically significant changes in the moment tensor during rupture for any of the five earthquakes. An appendix describes theMoment Tensor Rate Functions methodology in detail.The systematic bias between global and regional determinations of epicentral locations in Mexico must be resolved to enable plotting of asperities with aftershocks and geographic features. We have spatially shifted all of our results to regional determinations of epicenters. The best point source depths for the five earthquakes are all above 30 km, consistent with the idea that the down-dip edge of the seismogenic plate interface in Mexico is shallow compared to other subduction zones. Consideration of uncertainties in the focal mechanisms allows us to state that all five earthquakes occurred on fault planes with the same strike (N65°W to N70°W) and dip (15±3°), except for the smaller Playa Azul event at the down-dip edge which has a steeper dip angle of 20 to 25°. However, the Petatlan earthquake does prefer a fault plane that is rotated to a more east-west orientation—one explanation may be that this earthquake is located near the crest of the subducting Orozco fracture zone. The slip vectors of all five earthquakes are similar and generally consistent with the NUVEL-predicted Cocos-North America convergence direction of N33°E for this segment. The most important deviation is the more northerly slip direction for the Petatlan earthquake. Also, the slip vectors from the Harvard CMT solutions for large and small events in this segment prefer an overall convergence direction of about N20°E to N25°E.All five earthquakes share a common feature in the rupture process: each earthquake has a small initial precursory arrival followed by a large pulse of moment release with a distinct onset. The delay time varies from 4 s for the Playa Azul event to 8 s for the Colima event. While there is some evidence of spatial concentration of moment release for each event, our overall asperity distribution for the northern Mexico segment consists of one clear asperity, in the epicentral region of the 1973 Colima earthquake, and then a scattering of diffuse and overlapping regions of high moment release for the remainder of the segment. This character is directly displayed in the overlapping of rupture zones between the 1979 Petatlan event and the 1985 Michoacan aftershock. This character of the asperity distribution is in contrast to the widely spaced distinct asperities in the northern Japan-Kuriles Islands subduction zone, but is somewhat similar to the asperity distributions found in the central Peru and Santa Cruz Islands subduction zones. Subduction of the Orozco fracture zone may strongly affect the seismogenic character as the overlapping rupture zones are located on the crest of the subducted fracture zone. There is also a distinct change in the physiography of the upper plate that coincides with the subducting fracture zone, and the Guerrero seismic gap to the south of the Petatlan earthquake is in the wake of the Orozco fracture zone. At the northern end, the Rivera fracture zone in the subducting plate and the Colima graben in the upper plate coincide with the northernmost extent of the Colima rupture zone.     

Efficient generation of conditional simulations by chebyshev matrix polynomial approximations to the symmetric square root of the covariance matrix

Consider the problem of generating a realization y₁ of a Gaussian random field on a dense grid of points ₁ conditioned on field observations y₂ collected on a sparse grid of points ₂. An approach to this is to generate first an unconditional realization y over the grid =_{1 2}, and then to produce y₁ by conditioning y on the data y₂. As standard methods for generating y, such as the turning bands, spectral or Cholesky approaches can have various limitations, it has been proposed by M. W. Davis to generate realizations from a matrix polynomial approximations to the square root of the covariance matrix. In this paper we describe how to generate a direct approximation to the conditional realization y₁, on ₁ using a variant of Davis' approach based on approximation by Chebyshev polynomials. The resulting algorithm is simple to implement, numerically stable, and bounds on the approximation error are readily available. Furthermore we show that the conditional realization y₁ can be generated directly with a lower order polynomial than the unconditional realization y, and that further reductions can be achieved by exploiting a nugget effect if one is present. A pseudocode version of the algorithm is provided that can be implemented using the fast Fourier transform if the field is stationary and the grid ₁ is rectangular. Finally, numerical illustrations are given of the algorithm's performance in generating various 2-D realizations of conditional processes on large sampling grids.     

Reconstructing hydrological variability from testate amoebae analysis in Carpathian peatlands

Peatlands offer the potential for high resolution records of water balance over Holocene timescales, yet this potential is under-exploited in many areas of the world. Within Europe, peatlands are mostly confined to areas north of 55° N, but several areas of southern and eastern Europe contain small peatlands which may be suitable for palaeoclimatic reconstruction. In this paper we test the potential of peatlands in the Carpathian region for deriving quantified estimates of water table changes using testate amoebae analysis. A training set for palaeohydrological reconstruction from testate amoebae assemblages was obtained by collecting surface samples from 13 peatlands, including 9 from Hungary and 4 from Transylvania (Romania). Using a simple measure of mean annual water tables estimated from staining of PVC tape, we found that some peatlands were heavily influenced by runoff and groundwater, and were therefore not suitable as modern analogues of ombrotrophic climatically sensitive sites. The relationship between the testate amoebae assemblages in the modern samples and the environmental variables was explored using CCA. The CCA biplot showed that the most important variables are depth to water table and moisture content, confirming that hydrology is a key control on taxon distribution. pH was a secondary gradient. A transfer function for % moisture and depth to water table was established and applied to fossil assemblages from a sequence from Fenyves-tető, Transylvania, Romania. The reconstructed water table shows a number of variations which have parallels with other palaeoclimatic records from Europe and the North Atlantic prominent phases of higher water tables are associated with the periods 8000–8300 cal BP, 3000–2500 cal BP and after 600 cal BP. We suggest that these were periods of particular intensification of westerly airflow which affected eastern Europe as well as western and central Europe.