GPS单点测速的误差分析及精度评价 (英文)

Error sources which decrease the accuracy of GPS in absolute velocity determination have been changed since SA was turned off. Firstly, quantities of all kinds of error sources that influence velocity determination are analyzed. The potential accuracy of GPS absolute velocity determination is derived from both theory and field GPS data simulation. After that, two tests were carried out to evaluate the performance of GPS absolute velocity determination in the case of a static and an airborne GPS receiver and INS （Inertial Navigation System） instrument in kinematic mode. In static mode, the receiver velocity has been estimated to be several mm/s with the carrier-phase derived Doppler measurements, and several cm/s with the receiver generated Doppler measurements. In kinematic mode, GPS absolute velocity estimates are compared with the synchronized measurements from the high accuracy INS. The root mean square statistics of the velocity discrepancies between GPS and INS come up to dm/s. Moreover, it has a strong correlation with the acceleration or jerk of the aircraft.     

浅谈煤田地震勘探中影响纵向分辨率的主要因素

针对野外数据采集参数对分辨率的影响,对资料处理过程中静校正、反褶积、速度分析、Q补偿及谱白化等几个重要模块原理进行分析,指出在高速岩层中数据采集,应采用小药量激发,小道距,小排列,全三维,不组合、单点数字检波器接收。高密度多点采集,能够充分保护有效地质信息的高频成份,在更大范围内满足高分辨地震勘探的要求。资料处理过程中,在保证同相叠加和高信噪比的前提下,采用提高分辨率的模块以弥补地震波在传播过程中损失的高频成份。实践证明通过高保真、宽频带数据采集和高信噪比前提下的高分辨率数据处理,完全能够为矿井采掘提供高精度地震勘探资料。     

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.     

Earthquake rupture complexity due to dynamic nucleation and interaction of subsidiary faults

We numerically study the dynamic interaction of propagating cracks. It is assumed that propagating cracks can nucleate and drive subsidiary cracks because of shear strain enhancement near the propagating crack tips. The critical strain fracture criterion is assumed in the analysis. Intense interaction is expected to occur among the cracks. All the cracks are assumed to be parallel and antiplane strain deformation is assumed in the computation.In the interaction of two non-coplanar cracks, a strain shadow is formed in the neighborhood of each crack because of the strain release by the introduction of the crack. The growth of each crack is accelerated when the propagating tips of each crack are outside of the strain shadow of the other crack. In general, the crack tips enter the strain shadow, and the crack tips decelerate. The calculation shows that only one of the two cracks can continue to grow, and the other's growth is decelerated and arrested. If we can assume that the suite of cracks interact in a pairwise manner only, then this may suggest that only a limited number of cracks can continue to grow during the final stage of the rupture process. Hence the crack interaction causes complexity in dynamic earthquake faulting. The concepts of barrier and asperity have been employed by many researchers for the interpretation of complex seismic wave data. However, the physical realities of such concepts are obscure. Our calculations show that dynamic crack interactions can produce barriers and asperities in some cases; the crack tip deceleration or arrest due to the interactions among non-coplanar cracks can be interpreted as being due to a barrier. The dynamic coalescence among the coplanar cracks can be regarded as an asperity.Umeda found a localized area that strongly radiates high-frequency seismic waves in the epicentral areas of some large shallow earthquakes. He defined this as an earthquake bright spot. Our analysis implies that only a limited number of cracks continue to grow when many interactive cracks nucleate, and that all other cracks stop extending soon after nucleation. Hence, if the nucleation and termination of several cracks occur in a localized area, it will be observed seismologically as an earthquake bright spot. This is because it is theoretically known that the sudden termination of crack growth and dynamic crack coalescence efficiently emits high-frequency elastic waves.     

南海海盆区莫霍面分布规律及其对深部钻探的意义

钻遇莫霍面是人类一直以来的梦想。深海海底是地球上离莫霍面最近的地方，目前有研究推测南海是世界上莫霍面深度最浅的海域之一，但缺乏足够的直接证据。深反射地震探测可以直接揭示岩石圈的构造形态，是莫霍面探测的重要手段。本文基于长达15000 km的深反射多道地震剖面的解释、处理、制图和分析，结合前人的研究，形成了南海海盆区莫霍面反射特征和空间分布的初步认识。① 南海东部次海盆南部早期经历了较快速扩张，岩浆供应充足，受扩张停止后岩浆活动影响较小，基底平坦，地质构造相对简单，同时洋壳地震速度结构不存在异常，且有较强的广角莫霍面反射波和可识别的地幔顶部折射波，具备莫霍面钻探的基本条件。② 南海海盆不同区域的莫霍面反射强度存在较大差异。其中东部次海盆莫霍面反射最为强烈且清晰，西北次海盆次之，西南次海盆仅有零星出现的清晰莫霍面反射且可信度不高。③ 识别南海海盆区莫霍面地震反射长度超过3500 km，首次形成了海盆区深度域莫霍面地震反射空间分布图。与重力反演的莫霍面深度相比，利用深反射多道地震计算的莫霍面深度细节更为丰富，并且可以在垂向上清晰刻画莫霍面的结构。整体上，南海海盆区莫霍面地震反射强烈和可信度高的区域中，深度较浅的区域之一是东部次海盆南部，最浅处仅约9. 5 km，其中水深4. 01 km，洋壳厚度仅5. 54 km。综合判断，东部次海盆南部是南海重要的莫霍面钻探备选区，这对南海莫霍面钻探选址具有重要意义。     

含煤层地质环境下地震波场的数值模拟

本文对含低速煤层地质环境下弹性波场多波多分量地震资料进行了二维数值模拟研究,对人工边界反射进行了有效处理,频散效应得到了有效的压制,对几种不同激发与观测排列方式下的弹性波资料进行了模型计算与分析。     

Scattering of elastic waves by a fracture zone containing randomly distributed cracks

We theoretically study the scattering ofP, SV andSH waves by a zonal distribution of cracks, which simulates a fault fracture zone. An investigation is conducted how the geometrical properties of the crack distribution and the frictional characteristics of the crack surface are reflected in the attenuation and dispersion of incident waves, as well as in the amplitudes of the transmitted and reflected waves from the zone. If the crack distribution within the fault zone changes temporally during the preparation process of the expected earthquake, it will be important for earthquake prediction to monitor it, utilizing the scattering-induced wave phenomena.We consider the two-dimensional problem. Aligned cracks with the same length are assumed to be randomly distributed in a zone with a finite width, on which elastic waves are assumed to be incident. The distribution of cracks is assumed to be homogeneous and sparse. The crack surface is assumed to be stress-free, or to undergo viscous friction; the latter case simulates fluid-filled cracks. The opening displacement of the crack is assumed to be negligibly small. The idea of the mean wave formalism is employed in the analysis, and Foldy's approximation is assumed.When the crack surface is stress-free, it is commonly observed for every wave mode (P, SV andSH) that the attenuation coefficientQ ^–1 peaks aroundka1, the phase velocity is almost independent ofk in the rangeka<1 and it increases monotonically withk in the rangeka>1, wherek is the intrinsicS wavenumber anda is the half length of the crack. The effect of the friction is to shift the peak ofQ ^–1 and the corner of the phase velocity curve to the low wavenumber range. The high wavenumber asymptote ofQ ^–1 is proportional tok ^–1 independently of model parameters and the wave modes. If the seismological observation thatQ ^–1 ofS waves has a peak at around 0.5 Hz in the earth's crust is combined with our results, the upper limit of crack size within the crust is estimated about 4 km. The information regarding the transmitted and reflected waves, such as the high wavenumber limit of the amplitude of the transmitted wave etc., allows estimation of the strength of the friction.     

Transmission and reflection of Rayleigh waves at a thin low velocity vertical layer: A laboratory and theoretical study

A three-dimensional ultrasonic modelling was used to study the Rayleigh waves reflected from a thin low velocity vertical layer and those transmitted through the same zone. The reflection and transmission coefficients of the Rayleigh waves have been calculated by an approximate method based on the Green function technique. Nonrigid contact boundary conditions were used for the low velocity layer modelling. A rather good correspondence between the theoretical and experimental absolute values of the reflection and the transmission coefficients was obtained. Some discrepancies between experimental and theoretical results, in particular with regard to phase shifts, may be explained by a certain inadequacy of the theoretical model to the low velocity layer investigated by the ultrasonic modelling.     

Continuous air rate measurement in flotation cells: Some fundamental considerations

A robust and reliable sensor to measure gas (air) superficial velocity (J_g) continuously in flotation systems is introduced. It is based on the sampling of bubbles by buoyancy with the accumulating air allowed to exit through an orifice. At steady state, pressure drop is measured and related to the J_g by prior calibration. The continuous device is readily automated and extended to a multi-unit set-up. The sensor and data collections are described. The governing gas flow equation and models of the dynamic response to air flow rate set point change and fluctuations in froth depth are derived. Model predictions are confirmed against plant data. By taking a moving average, the J_g remains valid in the face of typical plant disturbances.     

Mais comment s'écoule donc un glacier ? Aperçu historique

Ice and snow have often helped physicists understand the world. On the contrary it has taken them a very long time to understand the flow of the glaciers. Naturalists only began to take an interest in glaciers at the beginning of the 19th century during the last phase of glacier advances. When the glacier flow from the upslope direction became obvious, it was then necessary to understand how it flowed. It was only in 1840, the year of the Antarctica ice sheet discovery by Dumont d'Urville, that two books laid the basis for the future field of glaciology: one by Agassiz on the ice age and glaciers, the other one by canon Rendu on glacier theory. During the 19th century, ice flow theories, adopted by most of the leading scientists, were based on melting/refreezing processes. Even though the word ‘fluid’ was first used in 1773 to describe ice, more the 130 years would have to go by before the laws of fluid mechanics were applied to ice. Even now, the parameter of Glen's law, which is used by glaciologists to model ice deformation, can take a very wide range of values, so that no unique ice flow law has yet been defined. To cite this article: F. Rémy, L. Testut, C. R. Geoscience 338 (2006).