首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 156 毫秒
1.
Viscoelastic modelling reveals that the interaction of compressional-wave velocity Cp, compressional-wave quality factor Qp, shear-wave velocity Cs, shear-wave quality factor Qs and Poisson's ratio as a function of time intercept τ and ray parameter p, is complicated; however, distinct, potentially diagnostic behaviours are seen for different combinations of viscoelastic parameters. Synthetic seismograms for three viscoelastic reservoir models show that variations in the Poisson's ratio produce visible differences when compared to the corresponding elastic synthetic seismograms; these differences are attributable to interaction of the elastic parameters with Qp and Qs. When the P-wave acoustic impedance contrast is small, viscoelastic effects become more apparent and more useful for interpretation purposes. The corresponding amplitude and net phase spectra reveal significant differences between the elastic and the viscoelastic responses. When P-wave reflectivities are large, they tend to dominate the total response and to mask the Q reflectivity effects. The attenuation effects are manifested as an amplitude decay that increases with both time and ray parameter. The sensitivity of the computed seismic responses for various combinations of viscoelastic parameters suggests the opportunity for diagnostic interpretation of τ-p seismic data. The interpretation of the viscoelastic parameters can permit a better understanding of the rock types and pore fluid distribution existing in the subsurface.  相似文献   

2.
The use of relaxation mechanisms has recently made it possible to simulate viscoelastic (Q) effects accurately in time-domain numerical computations of seismic responses. As a result, seismograms may now be synthesized for models with arbitrary spatial variations in compressional- and shear-wave quality factors (Q9, and Qs, as well as in density (ρ) and compressional- and shear-wave velocities (Vp, and Vs). Reflections produced by Q contrasts alone may have amplitudes as large as those produced by velocity contrasts. Q effects, including their interaction with Vp, Vs and p, contribute significantly to the seismic response of reservoirs. For band-limited data at typical seismic frequencies, the effects of Q on reflectivity and attenuation are more visible than those on dispersion. Synthetic examples include practical applications to reservoir exploration, evaluation and monitoring. Q effects are clearly visible in both surface and offset vertical seismic profile data. Thus, AVO analyses that neglect Q may produce erroneous conclusions.  相似文献   

3.
Characterization of shallow structures was performed by using different approaches analysing both P- and S-wave seismic data with different resolution. The refraction tomography provided P and S velocity models of the first 80 m, while the reflection seismic processing gives a reasonable stacking velocity field until 300 m depth for both P- and S-wave data. So, we estimated the Vp/Vs ratio and an empirical relationship between the two velocities. We characterised the shallow layers using tomographic velocity models and the deeper layers using seismic images with different resolution. The seismic images were obtained by conventional CMP reflection seismic processing and by a novel multi-refractor imaging technique.  相似文献   

4.
A new set of three-dimensional velocity models beneath Mt. Etna volcano is derived in the present work. We have used P- and S-wave arrivals from local earthquakes recorded at permanent and temporary seismic networks installed since 1980. A set of 1249 earthquakes recorded at more than four seismic stations was selected for traveltime inversion. The velocity models obtained by using different data selection criteria and parametrization display similar basic features, showing a high P-wave velocity at shallow depth in the SE quadrant, in close connection with a high gravimetric Bouguer anomaly. This area shares a low Vp/Vs ratio. High P-wave velocities and high Vp/Vs ratios are obtained along the central conduits, suggesting the presence of dense, intrusive magmatic bodies extending to a depth of about 20 km. The central intrusive core is surrounded by lower P-wave velocities. The relocated earthquake hypocenters also display the presence of an outward dipping brittle region, away from the central conduits, surrounding a ductile zone spatially related to the high P-wave velocity anomalies located in proximity to the central craters.  相似文献   

5.
We have measured the velocities and attenuations of compressional and shear waves in 29 water-saturated samples of sandstones and shales at a confining pressure of 60 MPa and at frequencies of about 0.85 MHz. The measurements were made using a pulse echo method in which the samples (diameter 5 cm, length 1.5 cm to 2.5 cm) were placed between perspex buffer rods inside a high-pressure cell. The velocity of each seismic wave was determined from the traveltime difference of equivalent phase points (corrected for diffraction effects) of the signals reflected from the top and from the base of each sample. Attenuation was determined in a similar way by comparison of the diffraction corrected amplitudes of the signals. The attenuation data are presented as ‘quality factors’: Qp and Qs for compressional and shear waves respectively. The results show that Qs is strongly correlated with Vs, that Qp is weakly correlated with Vp, and that Qp is strongly correlated with Qs. Qp is strongly dependent on the volume percentage of the assemblage of intra-pore minerals, whether they are clays or carbonates. It is concluded that the attenuation mechanism is due to the local fluid flow arising from the differential dilation of the solid rock frame and the intra-pore mineral assemblage, which is a result of their very different elastic moduli.  相似文献   

6.
Generalized inversion of the S-wave amplitude spectra from the strong-motion network data in the East-Central Iran has been used to estimate simultaneously source parameters, site response and the S-wave attenuation (Qs). In this regard, 190 three-component records were used corresponded to 40 earthquakes with the magnitudes M3.5–M7.3. These earthquakes were recorded at 42 stations in the hypocentral distance range from 9 to 200 km. The inverse problem was solved in 20 logarithmically equally spaced points in the frequency band from 0.4 to 15 Hz. The frequency-dependent site amplification was imposed, as a constraint, on two reference site responses in order to remove the undetermined degree of freedom in the inversion and obtain a unique inverse solution. Also, a geometrical spreading factor was assumed for removing the trade-off between geometrical spreading and anelastic attenuation. Different source parameters, such as seismic moment (M0), seismic energy (Es), corner frequency (fc) and Brune stress drop (Δσ), were estimated for each event by fitting an ω2 model to the spectra obtained from the inversion. The stress drop values of earthquakes, obtained in this research, are in good agreement with those of other studies. Also average site response values were correlated to the average shear wave velocities in the uppermost 30 m, in high and low frequency bands. The peak frequencies of site amplifications, estimated by the generalized inversion method, where in good agreement with those of horizontal to vertical (H/V) spectral ratios for the S-wave portion of records. However, no perfect matching in amplitude was obtained due to the deficiencies of the H/V ratio technique. By supposing a free shape for Q factor, a frequency dependent function was found, the logarithm of which could be approximated by a linear function, Q(f)=151f0.75. The uncertainties of model parameters have been evaluated by covariance matrix of least-square fit. The residuals were also analyzed in order to assess the validity of the model. The analysis of residuals with respect to magnitude and distance indicates that they are distributed normally with approximately zero mean. The robustness of the results has been studied concerning their sensitivities to the omission of different datasets, selected randomly from original database. The results obtained here can be used in predicting ground-motion parameters applying stochastic methods.  相似文献   

7.
The pressure dependence of P- and S-wave velocities, velocity anisotropy, shear wave splitting and crack-porosity has been investigated in a number of samples from different crustal rock types for dry and wet (water saturated) conditions. At atmospheric pressure, P-wave velocities of the saturated, low-porosity rocks (< 1%) are significantly higher than in dry rocks, whereas the differences for S-wave velocities are less pronounced. The effect of intercrystalline fluids on seismic properties at increased pressure conditions is particularly reflected by the variation of the Poisson's ratio because P-wave velocities are more sensitive to fluids than S-wave velocities in the low-porosity rocks. Based on the experimental data, the respective crack-density parameter (), which is a measure of the number of flat cracks per volume unit contained within the background medium (crack-free matrix), has been calculated for dry and saturated conditions. There is a good correlation between the calculated crack-densities and crack-porosities derived from the experimentally determined volumetric strain curves. The shear wave velocity data, along with the shear wave polarisation referred to a orthogonal reference system, have been used to derive the spatial orientation of effective oriented cracks within a foliated biotite gneiss. The experimental data are in reasonable agreement with the self consistent model of O'Connell and Budiansky (1974). Taking the various lithologies into account, it is clear from the present study, that combined seismic measurements ofV p andV s , using theV p V s -ratio, may give evidence for fluids on grain boundaries and, in addition, may provide an estimate on the in-situ crack-densities.  相似文献   

8.
We present the first systematic study of attenuation derived from the S-wave coda in the frequency range 1-32 Hz for the southern part of the Netherlands and its surroundings. For this we used two methods, the codaQ (Q c) method and the Multiple Lapse Time Window (MLTW) method. In the interpretation of the results both single and multiple scattering in a half space are considered. Our aim is to validate these interpretations in our region and to try to identify theeffects of attenuation due to intrinsic absoprtion (Q i)and scattering attenuation (Q s). For this we analyzedmore than 100 3-component high-quality digital seismograms from 43 crustalevents and 23 different stations in the Netherlands, Germany and Belgium.Coda Q results show smaller Q c (=Q 0fn) values for epicentral distances shorter than 25 km (Q 0=90) compared to larger epicentral distances (Q 0=190), but similar frequency dependence (f-0.9). Interpretation of MLTW results provided a seismic albedo smaller then 0.5, suggesting that the intrinsic absorption dominates over scattering in this region. Both Q i and Q s show similar frequency dependences as Q c. These results are comparable to those obtained in other areas, but we also show that more sophisticated models are required to remove ambiguities in the interpretation. For short lapse times and shortevent-station distances we find for the simple half space model a correspondinginterpretation of both methodologies, where Q c correspondsto Q t, suggesting that a model with single scattering in ahalf space is appropriate. For long lapse times and long event station distances, however, we find that the S-wave coda is, most probably, too much influenced by crust-mantel heterogenities and more sophisticated Qinversion models using larger data sets are required for more reliable attenuation estimates.  相似文献   

9.
High-frequency (≥ 2 Hz) Multi-channel Analysis of Love Waves (MALW) provides a practical way to determine velocity of horizontally polarized shear (SH) waves for a layered earth model up to 30 m below the ground surface in many geological settings. The information used in the MALW method is phase of Love waves. Information on amplitude of Love waves is not utilized in the MALW method. In this paper we present a method that uses information on amplitude of high-frequency Love waves to estimate quality factors (Qs) of near-surface materials. Unlike Rayleigh waves, attenuation coefficients (amplitude) of Love waves are independent of quality factors for P waves and are function of quality factors of Love waves. In theory, a fewer parameters make the inversion of attenuation coefficients of Love waves more stable and reduce the degree of nonuniqueness. We discussed sensitivity of an inversion system based on a linear relationship between attenuation coefficients and dissipation factors (1/Qs). The sensitivity analysis suggested that damping and constraints to an inversion system are necessary to obtain a smooth and meaningful quality factor model when no other information is available. We used synthetic and real-world data to demonstrate feasibility of inversion of attenuation coefficients of high-frequency Love-wave data acquired with the MALW method for quality factors with a linear, damped and constrained system.  相似文献   

10.
Minimization of seismic residuals does not guarantee uniqueness of the model, and this implies ambiguities in the inversion. Amplitude vs. offset (AVO) inversion does not lead to a unique solution of single elastic interface parameters unless converted and S-wave or critical angle reflections are available. Given the ambiguity of AVO inversion, this paper discusses the interaction between AVO and velocity estimation. The number of independent parameters necessary to describe an isolated reflection with AVO behaviour and residual velocity error is determined. Statistical analysis allows the establishment of an approximate equivalence of the effects of AVO and slight velocity variations; this equivalence cannot be solved without geological a priori information (kinematic equivalence). The data are then decomposed into compound events (i.e. sequences of N interfaces that follow each other at a fixed time lag). The decomposition is obtained by extrapolating the results of the analysis from narrowband to wideband data. Compound events decomposition demonstrates that AVO inversion is ambiguous, not only in the physical parameter space (P- and S-wave velocities, and density) but also kinematically. As an example of compound event decomposition, a medium is derived. This medium is geologically implausible but is kinematically equivalent.  相似文献   

11.
We derive S-wave attenuation characteristics, earthquake source parameters and site amplification functions at seismic stations used for earthquake early warning in the Irpinia–Basilicata region, using non-parametric spectral inversion of seismograms from 49 local events with M L = 1.5–3.1. We obtain relatively low Q values (Q 0 = 28 at a frequency of 1 Hz) in conjunction with a strong frequency-dependence (close to linear). The source spectra can be satisfactorily modeled using the omega-square model, with stress drops ranging between 0.01–2 MPa, and in the narrow magnitude range available for analysis, the source spectra seem to scale self-similarly. The local magnitude M L shows a linear correlation with moment magnitude M W, however with a systematic underestimation by about 0.5-magnitude units. The results obtained in this work provide important insights into the ground-motion characteristics that are required for appropriate seismic hazard assessment and are of practical relevance for a suite of applications, such as the calibration of ground-motion prediction equations or the correction for site amplification in earthquake early warning and rapid calculation of shake-maps for seismic emergency management.  相似文献   

12.
In this paper, on the basis of intensity data from 85 earthquakes occurred in China, the method of resolving overdetermined equations by using the damped least squares method is applied to inversing for the mean value of 1 s. S wave quality factor, the spreading factor and the S-wave acceleration at the foci of earthquakes in 7 regions of China. The relation between S-wave acceleration at the foci of earthquakes and magnitude is discussed. As an example, 2-DQ S distribution in Taiwan Province is obtained by inversion. It is found that there exists some corresponding relation between this distribution and Bouguer gravitational anomaly and seismicity. Preliminary analysis and discussion on results of the inversion, and estimation of error inQ S are made. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 202–211, 1991. The English version of this paper is improved by Professor Yushou Xie.  相似文献   

13.
Source, propagation path and site conditions are the factors affecting seismic ground motion. Consequently, recordings acquired at a seismic station are formed by the convolution of these three factors. In this work S-wave acceleration Fourier spectra of earthquakes recorded at local and regional scale, by the ITSAK accelerometric network for the period 2010–2016, are modeled as a product of source, propagation path (including geometric and anelastic attenuation) and site effects. The data set consists of 136 crustal earthquakes occurred in the broader Aegean area, with magnitudes 4.2?≤?Mw?≤?6.5 and epicentral distances 20 km?≤?R?≤?350 km, recorded at 112 broadband accelerometric stations installed at sites with various geologic conditions. Based on this data set, an iterative Gauss–Newton inversion method to solve the non-linear problem and retrieve the different terms of source, propagation path and site, is applied. This method uses an initial input model trying to find the best and at the same time a stable solution for the inverted parameters, which are, moment magnitude (Mw), corner frequency (fc), attenuation quality factor (Qs?=?Qofα), slope of the geometric attenuation (1/Rγ) and site transfer function (S(f)). The initial values of the starting model can be either known from other studies or inferred within a reasonable range. Depending on the level of knowledge on these input parameters, the associated standard deviation can be adjusted (large values for unknown parameters or small values for parameters which are well constrained). Results of the analyses exhibit satisfactory agreement of estimated source parameters with those proposed by seismological centers in Greece and propagation path properties similar to the ones determined in relevant previous studies for the same region. In addition, the site transfer functions obtained by the non-linear inversion are comparable with those calculated for the same sites using either standard spectral ratio or horizontal-to-vertical spectral ration (HVSR—receiver function) techniques. The aforementioned results are encouraging for reliable earthquake source parameters, propagation path properties and site effect assessment, in areas of intermediate to high seismicity.  相似文献   

14.
A quantitative AVO algorithm suitable for media with slow lateral parameter variations is developed. The method is based on a target-oriented inversion scheme for estimation of elastic parameters in a locally horizontally stratified medium. The algorithm uses band-limited PP reflection coefficients in the τ-p domain to estimate P- and S-wave velocities, densities and layer thicknesses. To obtain these reflection coefficients, a pre-processing involving the Radon transform and multiple attenuation is necessary. Furthermore, a macromodel for the velocities above the target zone must be found prior to the inversion. Various inversion tests involving synthetic data with white Gaussian noise and modelling errors that are likely to occur in conjunction with real data have been performed. In general, the inversion algorithm is fairly robust, since it is able to reproduce the main features of the reference model: main interface locations and relative contrasts in the three unknown layer parameters are recovered. From a test combining the effect of source directivity, one thin layer and 20% white Gaussian noise, it was found that neglect of the source directivity in the inversion caused the largest errors in the estimates. This indicates that it is very important either to eliminate the source directivity in a preprocessing step, or to take the directivity into account in the present algorithm. Despite these problems it was concluded that the inversion algorithm was able to reproduce the main features of the reference model.  相似文献   

15.
An inverse problem is one in which the parameters of a model are determined from measured seismic data. Important to the solution of inverse problems is the issue of whether or not a solution exists. In this paper we show, in a constructive manner, that a solution does exist to the specific inverse problem of determining the parameters of a horizontally stratified, lossless, isotropic and homogeneous layered system that is excited by a non-normal incidence (NNI) plane wave. Mode conversion between P- and S-waves is included. We develop a seven-step layer-recursive procedure for determining all of the parameters for layer j. These parameters are P-wave and S-wave velocities and angles of incidence, density, thickness, traveltimes, and reflection- and transmission-coefficient matrices. Downward continuation of data from the top of one layer to the top of the next lower layer is an important step in our procedure, just as it is in normal incidence (NI) inversion. We show that, in order to compute all parameters of layer j, we need to (and can) compute some parameters for layer j+ 1. This is a non-causal phenomenon that seems to be necessary in NNI inversion but is not present in NI inversion.  相似文献   

16.
Simultaneous estimation of effects of source, propagation path, and local site amplification was carried out using observed strong motion records in a frequency range from 0.8 to 20 Hz for the purpose of empirical evaluation of the local site effects in different geological conditions in the northwestern part of Turkey. The analyzed data are S-wave portions of 162 accelerograms from 39 shallow events observed at 14 sites of BYTNet array. A spectral separation method was applied to the observed S-wave spectra. The solutions for source spectra, inelasticity factor of propagation path for S-waves (Q s-value), and factor of site amplification at each site were obtained in a least squares sense. In the analysis, we assumed that the factor of the site amplification at a reference site is the same as that of theoretical amplification of S-waves to the soil model whose bottom layer has an S-wave velocity around 2.15 km/s. The estimated Q s-value of the propagation path is modeled as Q s(f)?=?87.4f0.78. The estimated site amplifications are characterized into three groups. The sites in the first group belong to rock site with no dominant peaks at a frequency range of 2 to 10 Hz. The second group of hard soil sites is characterized with moderately dominant peaks at a frequency of 5 Hz. The last group for soft soil sites has common peaks at a frequency of 4 Hz with larger amplitudes than those in the hard soil group. We, then, compare the amplifications with average S-wave velocity in top 30 m of the shallow S-wave profiles and proposed linear empirical formula between them at each frequency. We, furthermore, inverted the observed amplification factors into S-wave velocity and Q s-value profiles of the deep soil over the basement.  相似文献   

17.
Delineation of the top sedimentary structure and its Qs vs. Qp relationship using the travel-time difference of direct S and converted Sp phase is key to understanding the seismic hazard of any sedimentary basin area. We constructed filtered displacement waveforms from local ETNA Episensor acceleration recordings as well as local velocity recordings of aftershocks of the 2001 Bhuj earthquake recorded by the Kachchh seismological network of the National Geophysical Research Institute (NGRI), Hyderabad, India during 2001–2004. Stations are within 15–70km of epicenters, and the resulting displacement waveforms are generally simple, displaying prominent P, Sp, and S wave pulses. Particle motion of P and S waves suggest near-vertical raypaths consistent with preliminary depth estimates. The direct S wave on the horizontal component is characterized by lower frequency content than the converted Sp phase on the vertical component. This difference in frequency content between S and Sp phases can be explained in terms of different attenuation effects for P and S waves in the unconsolidated sediments. The Sp phase is generated by S-to-P phase conversion at the base of Mesozoic sediments of the Kachchh basin. Travel-time inversion (VELEST) of 2565 P and 2380 S arrivals from 658 well located aftershocks recorded at 8–14 three-component local seismic stations led to 1 D velocity models indicated very slow sediments in the upper 0–2 km depth range (Vp: 2.92 km/s and Vs: 0.90 km/s) and an increasing trend of velocities with depth at 2–40 km depth. The estimated sediment thicknesses beneath 12 accelerograph and 6 seismograph sites from the estimated velocity model and the travel-time difference between S and converted Sp phases reaches a maximum of (1.534 ± 0.117) km beneath Bandri (near the location of 2001 Bhuj mainshock) and attains a minimum sediment thickness of (0.858 ± 0.104) km beneath Ramvav and Burudia. The spectral ratios between Sp and S from 159 three-component accelerograms have been used to study seismic wave attenuation in the Kachchh rift basin. The estimated Qs vs. Qp relations for 12 accelerograph sites vary from Qs = 0.184 Qp (at Chobari) to Qs = 0.505 Qp (at Dudhai). For stations Chobari, Chopdwa, Jahawarnagar, Vondh and Tapar, the spectral ratio slopes and hence the calculated Qs vs. Qp relations are effectively the same, and the correlation coefficients are quite high (0.91–0.93). Stations Adhoi, Manfara, New Dudhai, Dudhai and Sikara have similar Qs vs. Qp relationships to each other and also have high correlation coefficients (0.78–0.87). The spectral ratios for stations Anjar and Ramvav are small and poorly constrained, resulting in less reliable Qs vs. Qp relations. This could be due to noisy data, fewer available waveforms, or scattering due to velocity heterogeneities and/or interface irregularities.  相似文献   

18.
Shear-wave velocity is a key parameter for calibrating monitoring time-lapse 4D seismic data during CO2-EOR (Enhanced Oil Recovery) and CO2 sequestration. However, actual S-wave velocity data are lacking, especially in 4D data for CO2 sequestration because wells are closed after the CO2 injection and seismic monitoring is continued but no well log data are acquired. When CO2 is injected into a reservoir, the pressure and saturation of the reservoirs change as well as the elastic parameters of the reservoir rocks. We propose a method to predict the S-wave velocity in reservoirs at different pressures and porosities based on the Hertz–Mindlin and Gassmann equations. Because the coordination number is unknown in the Hertz–Mindlin equation, we propose a new method to predict it. Thus, we use data at different CO2 injection stages in the Gao89 well block, Shengli Oilfield. First, the sand and mud beds are separated based on the structural characteristics of the thin sand beds and then the S-wave velocity as a function of reservoir pressure and porosity is calculated. Finally, synthetic seismic seismograms are generated based on the predicted P- and S-wave velocities at different stages of CO2 injection.  相似文献   

19.
Long-range seismic sounding carried out during the last few years on the territory of the U.S.S.R. has shown a basic inhomogeneity of the uppermost mantle, as well as evidence of regularities in the distribution of its seismic parameters. The following data were used: times and apparent velocities of P- and S-waves for investigation of mantle velocities, converted waves for seismic discontinuity model studies and wave attenuation for Q-factor estimation. Strong regularities were distinguished in the distribution of average seismic velocities for the uppermost mantle, in their dependence on the age and type of geostructure and on their position relative to the central part of the continent. Old platforms and the inner part of the continent are marked by velocities under the Mohorovi?i? discontinuity of more than 8.2–8.3 km s?1, young platforms and outer parts of the continent by 8.0–8.2 km s?1, and orogenic and rift zones by 7.8–8.0 km s?1. The difference becomes more pronounced at a depth of about 100–200 km: for the old platform mantle velocities of 8.5–8.6 km s?1 are typical; beneath the orogenic and rift areas, inversion zones with velocities less than 7.8 km s?1 are observed.The converted waves show fine inhomogeneities of the crust and uppermost mantle, the presence of many discontinuities with positive and negative changes of velocity, and anisotropy of seismic waves in some of the layers. Wave attenuation allowed the determination of the Q-factor in the mantle. It varied from one region to another but a close relation between Q and P-wave velocity is the main cause of its variation.  相似文献   

20.
An approximation is developed that allows mapped 4D seismic amplitudes and time‐shifts to be related directly to the weighted linear sum of pore pressure and saturation changes. The weights in this relation are identified as key groups of parameters from a petroelastic model and include the reservoir porosity. This dependence on groups of parameters explains the inherent non‐uniqueness of this problem experienced by previous researchers. The proposed relation is of use in 4D seismic data feasibility studies and inversion and interpretation of the 4D seismic response in terms of pore pressure and water saturation changes. A further result is drawn from analysis of data from the North Sea and West Africa, which reveals that the relative interplay between the effects of pore pressure and saturation changes on the seismic data can be simplified to the control of a single, spatially variant parameter CS/CP. Combining these results with those from published literature, we find that CS/CP = 8 appears to be a generality across a range of clastic reservoirs with a similar mean porosity. Using this CS/CP value, an in situ seismic‐scale constraint for the rock stress sensitivity component of the petroelastic model is constructed considering this component carries the largest uncertainty.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号