长江口三维潮流数值计算及动力分析   总被引：10，自引：1，他引：10  

于凤香  宋志尧  李瑞杰 《海洋湖沼通报》2003,(3):14-23

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

正交距离圆曲线拟合方法   总被引：2，自引：0，他引：2  

丁克良  刘全利  陈翔 《测绘科学》2008,(Z1)

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

水平井含水率上升影响因素   总被引：1，自引：1，他引：0  

凌宗发  王丽娟  李保柱  夏静  朱玉新 《现代地质》2008,22(1):86-90

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

On the multivariate total least-squares approach to empirical coordinate transformations. Three algorithms   总被引：2，自引：1，他引：1  

Burkhard Schaffrin  Yaron A. Felus 《Journal of Geodesy》2008,82(6):373-383

The multivariate total least-squares (MTLS) approach aims at estimating a matrix of parameters, Ξ, from a linear model (Y−E _Y = (X−E _X ) · Ξ) that includes an observation matrix, Y, another observation matrix, X, and matrices of randomly distributed errors, E _Y and E _X . Two special cases of the MTLS approach include the standard multivariate least-squares approach where only the observation matrix, Y, is perturbed by random errors and, on the other hand, the data least-squares approach where only the coefficient matrix X is affected by random errors. In a previous contribution, the authors derived an iterative algorithm to solve the MTLS problem by using the nonlinear Euler–Lagrange conditions. In this contribution, new lemmas are developed to analyze the iterative algorithm, modify it, and compare it with a new ‘closed form’ solution that is based on the singular-value decomposition. For an application, the total least-squares approach is used to estimate the affine transformation parameters that convert cadastral data from the old to the new Israeli datum. Technical aspects of this approach, such as scaling the data and fixing the columns in the coefficient matrix are investigated. This case study illuminates the issue of “symmetry” in the treatment of two sets of coordinates for identical point fields, a topic that had already been emphasized by Teunissen (1989, Festschrift to Torben Krarup, Geodetic Institute Bull no. 58, Copenhagen, Denmark, pp 335–342). The differences between the standard least-squares and the TLS approach are analyzed in terms of the estimated variance component and a first-order approximation of the dispersion matrix of the estimated parameters.  相似文献   

Rupture process of large earthquakes in the northern Mexico subduction zone     

Larry J. Ruff  Angus D. Miller 《Pure and Applied Geophysics》1994,142(1):101-172

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     

C. R. Dietrich  G. N. Newsam 《Mathematical Geology》1995,27(2):207-228

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.  相似文献   

Modelling net erosion responses to enviroclimatic changes recorded upon multisecular timescales     

Nazzareno Diodato   《Geomorphology》2006,80(3-4):164-177

Land use change has been recognized throughout the Earth as one of the most important factors influencing the occurrence of rainfall-driven geomorphological processes. However, relating the occurrence of historical soil erosion rates is difficult because of the lack of long-term research projects in river basins. Also, complex models are not adequate to reconstruct erosion rate changes because they require significant input data not always available on long timescales. Given the problems with assessing sediment yield using complex erosion models, the objective of this study is to explore a parsimonious scale-adapted erosion model (ADT) from the original Thornes and Douglas algorithms, which aims at reconstruction of annual net erosion (ANE) upon multisecular timescales. As a test site, the Calore River basin (3015 km² in southern Italy) provides a peculiar and unique opportunity for modelling erosion responses to climate and land cover changes, where input-data generation and interpretation results were also supported by documented hydrogeomorphological events that occurred before and after land deforestation. In this way, ANE_ADT-values were reconstructed for the period 1675–2004 by using precipitation indexes, complemented by recent instrumental records, and by using land cover statistics from documented agrarian sources. Pulses of natural sedimentation in the predeforestation period have been related to Vesuvius volcanic activity and changes in rainstorm frequency. After deforestation, the basin system became unstable with sudden fluctuations in the hydrogeomorphological regime contributing significantly to increased erosion and, in turn, sediment transport sequences via river drainage towards the Tyrrhenian coast.  相似文献   

Reconstructing hydrological variability from testate amoebae analysis in Carpathian peatlands     

C. Schnitchen  D. J. Charman  E. Magyari  M. Braun  I. Grigorszky  B. Tóthmérész  M. Molnár  Zs. Szántó 《Journal of Paleolimnology》2006,36(1):1-17

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.  相似文献   

Probabilistic local waveform inversion for moment tensor and hypocentral location   总被引：3，自引：0，他引：3  

Zoltán Wéber 《Geophysical Journal International》2006,165(2):607-621

  相似文献   

近45 a东北地区春季降水异常的气候特征   总被引：8，自引：0，他引：8  

杨素英  孙凤华  陈岩 《南京气象学院学报》2005,28(2):197-204

利用1959-2003年我国东北地区93站春季降水资料，将降水场分成5个区域，并在此基础之上分析了春季降水的时空变化特征，发现：降水量年际变化及长期趋势有明显的区域差异，呈东多西少的分布特征；西部是旱涝易发生区，近45a来降水量略有增多；降水量的周期振荡存在明显的区域差异。  相似文献