香港与珠三角产业优化合作与经济协调发展

本文分析了香港和珠三角的产业结构现状，提出了两地今后各自的产业发展路向，在此基础上，为发挥两地的整体优势，作者提出两地必须进行产业合作并得出结论：认为两地今后应发展结构性合作。香港当前应致力于制造业的升级转型，以此带动整个经济转型和多功能国际中心地位的强化；珠三角应适时向高新技术产业发展，２０１０年前应致力于工业化向成熟阶段的推进，建设成为一个高新技术产品输出基地。在此基础上，两地采取不同策略以实     

Windowed f-k spectra of a three-dimensional wavefield excited by a point source in a two-dimensional multilayered elastic medium

When full 3-D modelling is too costly or cumbersome, computations of 3-D elastic wave propagation in laterally heterogeneous, multilayered 2-D geological structures may enhance considerably our ability to predict strong ground motion for seismological and engineering purposes. Towards this goal, we extend the method based on the combination of the thin-layer finite-element and boundary-element methods (TLFE-BEM) and calculate windowend f - k spectra of the 3-D wavefield. The windowed f - k spectra are spatially localized spectra from which the local properties of the wavefield can be extracted. The TLFE-BEM is particularly suited for calculating the complete wavefield where surface waves are dominant in multilayered media. The computations are performed in the frequency domain, providing the f - k spectra directly. From the results for the 3-D wavefield excited by a point source in a 2-D multilayered, sloped structure, it can be said that the phase velocity of the fundamental-mode Rayleigh wave in a laterally heterogeneous multilayered medium, estimated from the windowed f - k spectra, varies with the location of the point source. For the model calculated in this article, the phase velocity varies between the value for the flat layered structure of the thick-layer side and that for the structure just under the centre of the window. The exact subsurface structure just under the centre of an array in a laterally heterogeneous medium cannot be obtained if we use the f - k spectral analysis assuming a flat layered structure.     

A new class of stratified viscoelastic models by analytical techniques

Multilayer, spherically stratified, self-gravitating relaxation models with a large number of layers (more than 100) can be dealt with analytically. Relaxation processes are studied for both Heaviside surface loads and tidal forcings. Simulations of the relaxation process of a realistic earth model with an incompressible Maxwell rheology show that models containing about 30 to 40 layers have reached continuum limits on all timescales and for all harmonic degrees up to at least 150 whenever an elastic lithosphere is present, irrespective of the viscosity profile in the mantle. In particular, fine-graded stratification of the shallow layers proves to be important for high harmonic degrees in these models. The models produce correct long-time (fluid) limits. It is shown that differences in the transient behaviour of the various models are due to the applied volume-averaging procedure of the rheological parameters. Our earlier proposed hypothesis that purported shortcomings in the fundamental physics of (discrete) normal-mode theory are artificial consequences of numerical inaccuracies, theoretical misinterpretations and the use of incomplete sets of normal modes is reinforced by the results presented. We show explicitly that the models produce both continuous behaviour resulting from continuous rheological stratifications and discrete behaviour resulting from sharp density contrasts, as at the outer surface and the core-mantle boundary. The differences between volume-averaged models and fixed-boundary contrast models are outlined. Reducing many-layer models with a volume-averaging procedure before employing a normal mode analysis is both economical and highly accurate on all timescales and for all spherical harmonic degrees. The procedure minimizes the chances of missing contributing modes, while using models with more layers will not result in any substantial increase of accuracy.     

Rayleigh-wave dispersion function for a transversely isotropic layered spherical earth using the Thomson-Haskell matrix method

We consider the two coupled differential equations of the two radial functions appearing in the displacement components of spheroidal oscillations for a transversely isotropic (TI) medium in spherical coordinates. Elements of the layer matrix have been explicitly written—perhaps for the first time—to extend the use of the Thomson-Haskell matrix method to the derivation of the dispersion function of Rayleigh waves in a transversely isotropic spherical layered earth. Furthermore, an earth-flattening transformation (EFT) is found and effectively used for spheroidal oscillations. The exponential function solutions obtained for each layer give the dispersion function for TI spherical media the same form as that on a flat earth. This has been achieved by assuming that the five elastic parameters involved vary as r ^p and that the density varies as r ^p-2, where p is an arbitrary constant and r is the radial distance. A numerical illustration with p = - 2 shows that, in spite of the inhomogeneity assumed within layers, the results for spherical harmonic degree n , versus time period T , obtained here for the Primary Reference Earth Model (PREM), agree well with those obtained earlier by other authors using numerical integration or variational methods. The results for isotropic media derived here are also in agreement with previous results. The effect of transverse isotropy on phase velocity for the first two modes of Rayleigh waves in the period range 20 to 240 s is calculated and discussed for continental and oceanic models.     

Some remarks about the degree-one deformation of the Earth

The degree-one deformation of the Earth (and the induced discrepancy between the figure centre and the mass centre of the Earth) is computed using a theoretical approach (Love numbers formalism) at short timescales (where the Earth has an elastic behaviour) as well as at long timescales (where the Earth has a viscoelastic or quasi-fluid behaviour). For a Maxwell model of rheology, the degree-one relaxation modes associated with the viscoelastic Love numbers have been investigated: the M_o mode does not exist and there is only one transition mode (instead of two) generated by a viscosity discontinuity.
The translations at each interface of the incompressible layers of the earth model [surface, 670 km depth discontinuity, core-mantle boundary (CMB) and inner-core boundary (ICB)] are computed. They are elastic with an order of magnitude of about 1 mm when the excitation source is the atmospheric continental loading or a magnetic pressure acting at the CMB. They are viscoelastic when the earth is submitted to Pleistocene deglaciation, with an order of magnitude of about 1 m. In a quasi-fluid approximation (Newtonian fluid) because of the mantle density heterogeneity their order of magnitude is about 100 m (except for the ICB, which is in quasi-hydrostatic equilibrium at this timescale).     

地球内核的径向振动

研究了地球内核围绕其平衡位置的径向振动，即所谓Ｓｌｉｃｈｔｅｒ模，求出了这一现象的周期。在求解导出的常微分方程组时采用了切贝谢夫配点法，容易而成功地解决了方程在地心处的奇异性带来的困难。所得结果与国外目前公认的结果一致     

次天气尺度与中尺度正、斜压模演变揭示的武汉暴雨过程

利用滤波原理提取出大气流场中的次天气尺度和中尺度信息，再把大尺度和中尺度水平风场分别分解为正压分量（垂直平均）和斜压分量（扰动）流场。对1998年7月21～22日发生在武汉附近的强暴雨过程进行了次天气尺度与中尺度流场正、斜压分量演变特征的分析。结果表明：次天气尺度与中尺度流场正压分量的演变与此次强暴雨的酝酿、发展和消亡具有内在的联系；次天气尺度与中尺度流场高层200hPa斜压分量很强，低层850hPa正压分量很强；次天气尺度与中尺度流场斜压性占主导地位，随着暴雨的发展，中尺度流场的正压性减弱而斜压性进一步增强，而次天气尺度流场的正压性增强而斜压性减弱。以上结论对于揭示中尺度暴雨过程发生发展的本质有一定的意义。     

煤中有害微量元素的赋存状态

基于前人资料的分析总结，系统探讨了煤中Ag等26种有害微量元素的赋存状态。研究结果表明：煤中有害元素或多或少都与无机质、有机质有联系，只是联系程度不同。有些在地壳中被认为彼此关系不大的元素，在煤中却密切相关。煤中B、Be、Br等主要以有机相赋存，其它有害元素主要与矿物质有关。