基于姿态变化的航空频率域电磁法仪器偏置实时校正方法研究

磁偶极子的航空频率域电磁法仪器在飞行测量的过程中由于仪器偏置的存在,且仪器偏置会随着外部气压、温度等环境因素以及收发线圈晃动的影响而呈现非线性变化,使得观测数据出现误差,因此需要对仪器偏置进行校正.而传统的在测线飞行前后将仪器抬至高空的"零场值"标定方法具有成本高、受测区环境影响大以及采用线性插值获取测线飞行过程中仪器偏置的精度低等缺点.本文根据仪器偏置与仪器姿态角变化无关的特性,通过测得仪器的姿态角信息,在满足重叠偶极子模型的条件下,实现对仪器偏置的高精度实时校正.模型仿真结果表明,在30m常规飞行高度下,该方法实时测得的仪器偏置精度接近于110m高空测得的精度;校正后仪器偏置的绝对误差与理论二次场的比值即相对误差小于5%,满足反演大地电导率的精度要求.该方法不仅减少了飞行的工作量,降低了飞行成本和飞行难度,而且可更加精确地获得测线飞行过程中仪器偏置的非线性变化值,提高航测数据的观测精度.     

Prediction of response spectra via real-time earthquake measurements

The development and implementation of an earthquake early warning system (EEWS), both in regional or on-site configurations can help to mitigate the losses due to the occurrence of moderate-to-large earthquakes in densely populated and/or industrialized areas. The capability of an EEWS to provide real-time estimates of source parameters (location and magnitude) can be used to take some countermeasures during the earthquake occurrence and before the arriving of the most destructive waves at the site of interest. However, some critical issues are peculiar of EEWS and need further investigation: (1) the uncertainties on earthquake magnitude and location estimates based on the measurements of some observed quantities in the very early portion of the recorded signals; (2) the selection of the most appropriate parameter to be used to predict the ground motion amplitude both in near- and far-source ranges; (3) the use of the estimates provided by the EEWS for structural engineering and risk mitigation applications.In the present study, the issues above are discussed using the Campania–Lucania region (Southern Apennines) in Italy, as test-site area. In this region a prototype system for earthquake early warning, and more generally for seismic alert management, is under development. The system is based on a dense, wide dynamic accelerometric network deployed in the area where the moderate-to-large earthquake causative fault systems are located.The uncertainty analysis is performed through a real-time probabilistic seismic hazard analysis by using two different approaches. The first is the Bayesian approach that implicitly integrate both the time evolving estimate of earthquake parameters, the probability density functions and the variability of ground motion propagation providing the most complete information. The second is a classical point estimate approach which does not account for the probability density function of the magnitude and only uses the average of the estimates performed at each seismic station.Both the approaches are applied to two main towns located in the area of interest, Napoli and Avellino, for which a missed and false alarm analysis is presented by means of a scenario earthquake: an M 7.0 seismic event located at the centre of the seismic network.Concerning the ground motion prediction, attention is focused on the response spectra as the most appropriate function to characterize the ground motion for earthquake engineering applications of EEWS.     

Hypocenter location of mixed events in real-time processing systems

Automatic pickings in earthquake real-time monitoring systems often contain noise bursts and/or phases of different event(s) occurring almost simultaneously. Typically, a locator uses these picks as P and S waves arrival times coming from a single event and, therefore, should be complemented by a distinctive phase association logic. The method we propose manages to automatically associate data related to different events and eliminates the influence of spoiled data from single events. The method is based on “network beamforming”, a robust and stable algorithm, which utilizes a hypocenter grid search for the stack maximum of a set of complex exponents applied to the P phase readings. The algorithm separates the residual outliers and then uses them for location. If successful, a hypocenter is established for the interfering event. The solutions obtained are overall robust and independent from the estimate of origin times. The preliminary epicenter for the grid search is provided by the intersection of perpendicular bisectors in the modified “arrival order algorithm” or by the modified “T^now” algorithm, which uses non-arrival information. We applied this method to automatic first arrival phase readings of 915 events registered by the Hi-net Japan seismic network and our results are statistically promising. Here, we present two interesting and complicated examples.     

厚壳贻贝(Mytilus coruscus)金属硫蛋白MT-10:cDNA克隆、结构分析及铜离子胁迫下的表达

近年来我国沿海海域重金属污染日趋严重,软体动物尤其是双壳贝类响应重金属污染胁迫的研究发展成为目前的研究热点。金属硫蛋白是一类广泛存在于生物体内富含半胱氨酸(Cys)的小分子蛋白,具有结合金属离子的能力,能有效调节生物体细胞内的金属离子平衡,具有清除自由基,重金属解毒的功能,关于厚壳贻贝MT蛋白的研究尚未见报道。本文首次利用RT-PCR和c DNA末端快速扩增技术克隆了厚壳贻贝MT-10蛋白c DNA全长(Mc-MT-10),其包含101bp的5′UTR区,108bp的3′UTR区以及222bp的ORF区,编码73个氨基酸。经比对发现Mc-MT-10富含Cys,且具有9个软体动物中常见的CXC结构以及CKCXXXCXCX结构,系统进化树中同科或属中MT-10聚合成一个分支,显示出高度保守性。采用Q-PCR技术分析了MT-10 mRNA在Cu~(2+)胁迫下厚壳贻贝外套膜及消化腺中表达变化,结果表明Mc-MT-10 mRNA在Cu~(2+)胁迫下两种组织中的表达量急剧升高,说明Cu~(2+)能够诱导Mc-MT-10 mRNA的表达,并且表现出诱导具有阈值性特点。研究结果可为深入研究厚壳贻贝MT蛋白的功能奠定基础,为厚壳贻贝抗重金属污染新品种选育提供数据支持。     

Fast computation of the transient motions of moving vessels in irregular ocean waves

A fast time-domain method is developed in this paper for the real-time prediction of the six degree of freedom motions of a vessel traveling in an irregular seaway in infinitely deep water. The fully coupled unsteady ship motion problem is solved by time-stepping the linearized boundary conditions on both the free surface and body surface. A velocity-based boundary integral method is then used to solve the Laplace equation at every time step for the fluid kinematics, while a scalar integral equation is solved for the total fluid pressure. The boundary integral equations are applied to both the physical fluid domain outside the body and a fictitious fluid region inside the body, enabling use of the fast Fourier transform method to evaluate the free surface integrals. The computational efficiency of the scheme is further improved through use of the method of images to eliminate source singularities on the free surface while retaining vortex/dipole singularities that decay more rapidly in space. The resulting numerical algorithm runs 2–3 times faster than real time on a standard desktop computer. Numerical predictions are compared to prior published results for the transient motions of a hemisphere and laboratory measurements of the motions of a free running vessel in oblique waves with good agreement.     

The QPO phenomenon

Aspects of QPOs in the X‐ray flux of low mass X‐ray binaries are reviewed, with particular attention to the comparison between the QPOs in black holes and those in neutron stars. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Peculiarities of processing VLBI observations of pulsars

We determined the features of pulsars that were disregarded in standard amplitude-calibration procedures for VLBI observations. We suggest additional amplitude-calibration methods. These methods were used to process the VLBI observations of the pulsar PSR B0329+54 carried out with the HALCA ground—space interferometer. Data from the space radio telescope are corrected for a nonuniform reception band. We estimated the diameter of the scattering disk for this pulsar at a frequency of 1600 MHz: \( < 1\mathop .\limits^{''} 8 \times 10^{ - 3}\).     

The X-ray burst detected in 1991 from the transient source SAX J1747.0-2853

We present the observations of an intense X-ray burst from the recurrent transient source SAXJ1747.0-2853 located near the Galactic center. The burst was detected by the ART-P telescope onboard the Granat observatory on October 20, 1991. The burst time profile exhibits features that clearly point to an increase in the photospheric radius of the neutron star at the burst onset. This increase is attributable to an expansion of its outer hydrogen-rich atmospheric envelope under radiation pressure. After hydrogen-envelope outflow and photospheric contraction, the radiation flux emerging from deep within the star continued to rise, which several seconds later led to a recurrent, weaker photospheric expansion attributable to the outflow of the outer helium-envelope layers. Based on the described picture, we determined the distance to the source, d=7.9±0.4 kpc. No radiation was detected by the ART-P telescope from the source in quiescence. Actually, the source itself was discovered only seven years later by the BeppoSAX satellite during its X-ray activity.     

High-precision timing of PSR J1600−3053

We present the results of a high-precision timing campaign directed at the binary millisecond pulsar J1600−3053. Submicrosecond pulsar timing has long been the domain of bright, low dispersion measure millisecond pulsars or large diameter telescopes. This experiment, conducted using the Parkes radio telescope in New South Wales, Australia, and utilizing the latest baseband recording hardware, has allowed this pulsar, although distant and faint, to present residuals to a model of its spin behaviour of 650 ns over a period of more than 2 yr. We have also constrained the orbital inclination via Shapiro delay to be between 59° and 70° to 95 per cent confidence and obtained a scintillation velocity measurement indicating a transverse velocity less than 84 km s⁻¹. This pulsar is demonstrating remarkable stability comparable to, and in most cases improving upon, the very best long-term pulsar timing experiments. If this stability is maintained, the current limits on the energy density of the stochastic gravitational wave background will be reached in four more years.