首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
随着页岩气开采、地热能源开采、CO2封存和水库蓄水等工业活动的开展,由此产生的诱发地震危害问题日益严重。为了对诱发地震进行有效监测和管控,地震学家设计了"交通灯"系统,当地震活动达到一定阈值时,相关人员可根据系统警报及时采取应对措施,减轻地震灾害。本文调研了近年来多个国家的"交通灯"系统研究和应用进展,介绍了其发展历史、设置原理及应用案例。"交通灯"系统的发展和完善应综合考虑震级、震动强度、b值和断层分布,以及建筑结构、城镇距离、人口密度、公众反应等社会因素。同时,在诱发地震频发区域,应建设密集的专用监测台网,提高微震检测和定位能力,并引入模板匹配和人工智能等新的自动处理方法,及时产出高完备性、高精度的微震目录,对微震的时空演化进行有效监测,形成时效性较强的"交通灯"系统。  相似文献   

2.
A layeredP- andS-wave velocity model is obtained for the Friuli seismic area using the arrival time data ofP- andS-waves from local earthquakes. A damped least-squares method is applied in the inversion.The data used are 994P-wave arrival times for 177 events which have epicenters in the region covered by the Friuli seismic network operated by Osservatorio Geofisico sperimentale (OGS) di Trieste, which are jointly inverted for the earthquake hypocenters andP-wave velocity model. TheS-wave velocity model is estimated on the basis of 978S-wave arrival times and the hypocenters obtained from theP-wave arrival time inversion. We also applied an approach thatP- andS-wave arrival time data are jointly used in the inversion (Roecker, 1982). The results show thatS-wave velocity structures obtained from the two methods are quite consistent, butP-wave velocity structures have obvious differences. This is apparent becauseP-waves are more sensitive to the hypocentral location thanS-waves, and the reading errors ofS-wave arrival times, which are much larger than those ofP-waves, bring large location errors in the joint inversion ofP- andS-wave arrival time. The synthetic data tests indicated that when the reading errors ofS-wave arrivals are larger than four times that ofP-wave arrivals, the method proposed in this paper seems more valid thanP- andS-wave data joint inversion. Most of the relocated events occurred in the depth range between 7 and 11 km, just above the biggest jump in velocity. This jump might be related to the detachment line hypothesized byCarulli et al. (1982). From the invertedP- andS-wave velocities, we obtain an average value 1.82 forV p /V s in the first 16 km depth.  相似文献   

3.
The paper discusses microseismic monitoring during oil well stimulation by hydraulic fracturing, an emergent technology used for hydraulic fracturing layer control. The passive monitoring is a new widely developing technology of HFL control. The main factor affecting the results of passive seismic monitoring is the event location accuracy. The passive monitoring acquisition system utilizes one three components seismic probe deployed into the observation well. To evaluate the location accuracy of induced events for one observation well we applied traditional kinematic approach based on picking of earthquake P- and S-waves arrivals. The influence of geometric parameters of geophones location in a borehole, their quantity, picking errors of waves arrivals on the accuracy of microearthquakes location is studied.  相似文献   

4.
Scaling relations for seismic events induced by mining   总被引：1，自引：0，他引：1  
The values of seismic moment andS-wave corner frequency from 1575 seismic events induced in South African, Canadian, Polish, and German underground mines were collected to study their scaling relations. The values ofP-wave corner frequency from 649 events were also available. Seismic moments of these events range from 5*103 to 2*1015 N·m (moment magnitude is from –3.6 to 4.1), theS-wave corner frequency ranges from 0.7 to 4438 Hz, and theP-wave corner frequency is between 5 and 4010 Hz. The slope of a regression line between the logarithm ofS- andP-wave corner frequencies is equal to one, and the corner frequencies ofP waves are higher than those ofS waves on the average by about 25 percent. In studies of large and moderate earthquakes it has been found that stress drop is approximately independent of the seismic moment, which means that seismic moment is inversely proportional to the third power of corner frequency. Such a behavior was confirmed for most of the data considered here. A breakdown in the similarity betwen large and small events seems to occur for the events with moment magnitude below –2.5. The average values of seismic moment referred to the same range of corner frequency, however, are vastly different in various mining areas.  相似文献   

5.
Interpretation ofP/S-wave successions is used in induced or passive microseismicity. It makes the location of microseismic events possible when the triangulation technique cannot be used. To improve the reliability of the method, we propose a technique that identifies theP/S-wave successions among recorded wave successions. A polarization software is used to verify the orthogonality between theP andS polarization axes. The polarization parameters are computed all along the 3-component acoustic signal. Then the algorithm detects time windows within which the signal polarization axis is perpendicular to the polarization axis of the wave in the reference time window (representative of theP wave). The technique is demonstrated for a synthetic event, and three application cases are presented. The first one corresponds to a calibration shot within which the arrivals of perpendicularly polarized waves are correctly detected in spite of their moderate amplitude. The second example presents a microseismic event recorded during gas withdrawal from an underground gas storage reservoir. The last example is chosen as a counter-example, concerning a microseismic event recorded during a hydraulic fracturing job. The detection algorithm reveals that, in this case, the wave succession does not correspond to aP/S one. This implies that such an event must not be located by the method based on the interpretation of aP/S-wave succession as no such a succession is confirmed.  相似文献   

6.
—We analyze source parameters of microseismic events (M < ? 1) associated with high flow-rate water injections in a shale formation at a depth of 220 m. Two types of events were observed several hundred impulsive events with clear P- and S-wave arrivals, and continuous emissions with peaked spectra detected well into the experiment. For a representative collection of impulsive events, an ?2 model provided satisfactory fits to displacement spectra corrected for attenuation, and average quality factors of 34 and 15 were obtained for P and S waves. P-wave first motion analysis and E S ?/E P ratios indicated the existence of a non-double-couple component in some events, particularly early in the experiment. A clear difference was observed for estimates of stress release parameters as non-double-couple events had smaller stress drops and apparent stresses. The seismic efficiency of double-couple and non-double-couple events was limited to 0.9% and 0.05% respectively, with average values being 0.25% and 0.02%. A comparison of our results with those reported for a similar magnitude range in a hard-rock formation indicates considerably smaller estimates of stress drop and apparent stress in our case while seismic efficiencies are comparable.  相似文献   

7.
Shear wave velocity modelling in crustal rock for seismic hazard analysis   总被引：2，自引：1，他引：2  
P-wave velocity data along with the thickness of sedimentary and crystalline layers within bedrock were collected from all global regions and presented in the Global Crustal Model CRUST2.0, published in 2001. This well-organised database provides invaluable potential contributions towards future seismic hazard modelling, particularly for stable continental regions (SCRs), where there is a scarcity of representative strong motion records for conventional modelling purposes. The P-wave velocity information presented in CRUST2.0 has been converted herein to S-wave velocity information. The latter is especially important for purposes of seismic hazard modelling. The value of the CRUST2.0 model has therefore been greatly enhanced by the important findings presented and further developed in this paper. By making the best use of available information on crustal conditions, the amplification behaviour of seismic waves affecting a region, an area or a site for any given earthquake scenario may be predicted. The developed methodology, which is intended for worldwide applications, has been illustrated by case studies in which model S-wave velocity profiles were developed for different geological regions within North America. The model profiles were found to be in excellent agreement with field measurements reported for each respective region.  相似文献   

8.
The thin-layer build of the Carpathian Foredeep Miocene formations and large petrophysical parameter variation cause seismic images of gas-saturated zones to be ambiguous, and the location of prospection wells on the basis of anomalous seismic record is risky. A method that assists reservoir interpretation of standard recorded seismic profiles (P waves) can be a converted wave recording (PS waves). This paper presents the results of application of a multicomponent seismic survey for the reservoir interpretation over the Chałupki Dębniańskie gas deposit, carried out for the first time in Poland by Geofizyka Kraków Ltd. for the Polish Oil and Gas Company. Seismic modeling was applied as the basic research tool, using the SeisMod program based on the finite-difference solution of the acoustic wave equation and equations of motion. Seismogeological models for P waves were developed using Acoustic Logs; S-wave model (records only from part of the well) was developed on the basis of theoretical curves calculated by means of the Estymacja program calibrated with average S-velocities, calculated by correlation of recorded P and PS wavefields with 1D modeling. The conformity between theoretical and recorded wavefields makes it possible to apply the criteria established on the basis of modeling for reservoir interpretation. Direct hydrocarbon indicators (bright spots, phase change, time sag) unambiguously identify gas-prone layers within the ChD-2 prospect. A partial range of the indicators observed in the SW part of the studied profile (bright spot that covers a single, anticlinally bent seismic horizon) points to saturation of the horizon. The proposed location is confirmed by criteria determined for converted waves (continuous seismic horizons with constant, high amplitude) despite poorer agreement between theoretical and recorded wavefields.  相似文献   

9.
—?The aim of our study consists of analyzing potentially non-double-couple seismic events recorded at regional distances. In order to define the nature of the seismic source, a moment tensor inversion is carried out as this method is general enough not to initially constrain the source mechanism. In this paper we present an application to a seismic event induced by a mine collapse which occurred near the town of Halle in Germany. Because of its induced nature, many parameters such as the location and geometry of this seismic source are known. This information allows us to test the influence of inadequate propagation modeling on the moment tensor obtained from the inversion. Green's functions have been computed with the reflectivity method in a flat layered medium, using the European model EurID (Du et? al., 1998; Dufumier et al., 1997). From the inversion of P-wave seismograms recorded by the German Regional Seismic Network will, we obtained a source time function which can be decomposed into two subevents. The first one has a large isotropic part and a deviatoric mechanism with near vertical nodal planes. No volume change is observed for the second subevent, but a deviatoric component opposite of the first one. The addition of S-waves does not change the results of the inversion which are stable. Surface waves were not used because of their poor dispersion curves. Based on the moment tensor obtained from these inversions, the physical process at the source is compatible with a large cavity collapse.  相似文献   

10.
The attenuation of amplitude is seen in seismic waves which pass through the central region of the Aso caldera, in Kyushu, Japan. It is also recognized from spectral analysis of seismic waves that the higher frequencies of the P-wave are reduced in the waves which pass through the central region of the caldera. It is shown that the relative attenuation increases remarkably for the frequency range of 5 to 10 Hz. The specific attenuation factor Q of the P-wave train is about 100. From the surface projection of the ray paths with low Q values through the Aso caldera to each station, the attenuating region is located beneath the center of the caldera, extending to the north of the central cones. In conjunction with the low Q value of the P-wave and the decreases of S-wave amplitudes, the relative P-wave residual times have comparatively large values for seismic waves passing through the central region beneath the caldera. In order to attempt to provide additional information on the depth configuration of the attenuating material, the ray paths of P-wave's first arrivals are located in three-dimensional space. It indicates that the low-velocity material is located beneath the center of the caldera at depths of about 6 to 9 km. However, lowvelocity anomalies above the depth of 6 km and below the depth of 15 km were not able to be detected, because most of the available seismic ray paths had crossed the caldera at depths of about 6 to 15 km. Furthermore, the relative residual times have numerous errors resulting from incorrect hypocenter locations, origin times, inhomogeneities in the structure and uncertainty of the velocity structure. At shallow depths in the Aso caldera, refraction or reflection studies are required for an accurate estimate of the structure and more detailed properties of the attenuating material.  相似文献   

11.
We demonstrate how multiples, generated at the interfaces of plane parallel beds, modify the propagation characteristics of an originally coherent seismic wave. For waves propagating at an angle to the bedding plane we find that theSV andP-waves couple so that neither is a pure mode. TheSH-wave, while modified in its propagation characteristics by multiples, remains a pure mode. The coupling ofSV-multiples into the quasi-P-mode appears weaker than the coupling ofP-wave multiples into the quasi-SV mode; at least this is so for the two simple cases of (a) density fluctuations only and (b) correlatedV p andV s fluctuations which conserve Poisson's ratio.We also find that the coupling is sensitive to both the angle of propagation and frequency. In addition there is a cut-off angle forP-wave multiples influencing the quasi-SV mode. Propagation angles larger than the cut-off permit theP-multiples to modify the phase of the quasi-SV mode, but not its effective attenuation. No such cut-off effect is found for SV-multiples influencing the quasi-P mode, whose angle-dependent and frequency-dependent phase distortion and effective attenuation are influenced both byP-wave multiples andSV-multiples.In view of the mathematical complexity of the expressions describing the phase, and effective attenuation of modes when allowance is made forP-andS-wave multiples, we strongly advocate numerical coding of the major mathematical formulae. By so doing a systematic study can be undertaken of the frequency and offset dependence of seismic waves as a function of seismic source input and power spectral behavior of the fluctuations in density and elastic constants of beds. It is our opinion that the full mathematical expressions are too involved to permit an analytic, systematic investigation to be given of the phase and attenuation of seismic waves with any degree of sophistication or generality.  相似文献   

12.
The source parameters, moment, stress drop and source dimension are estimated for 61 events from the January 1975 Brawley earthquake swarm. Earthquakes studied range in local magnitude from 1.0 to 4.7. Stress drops range from 1 to 636 bars and increase with source depth. It is estimated that the sedimentary structure of the Imperial Valley amplifies shear waves by a factor of 2 to 3 in addition to the free surface amplification of 2. Estimates of moment from 10 sec surface waves are 4 to 6 times larger than the moment estimated from the relatively flat part of the local body wave spectrum at 1 sec. This may be due to after-slip on the fault, a long thin fault, or partial stress drop. It is shown that the experimentally determined ratio of stress drop to apparent stress should be approximately 4.0 when spectrum integration is used to obtainS-wave energy and theP-wave energy is 1/3 theS-wave energy.  相似文献   

13.
Seismic anisotropy which is common in shale and fractured rocks will cause travel-time and amplitude discrepancy in different propagation directions. For microseismic monitoring which is often implemented in shale or fractured rocks, seismic anisotropy needs to be carefully accounted for in source location and mechanism determination. We have developed an efficient finite-difference full waveform modeling tool with an arbitrary moment tensor source. The modeling tool is suitable for simulating wave propagation in anisotropic media for microseismic monitoring. As both dislocation and non-double-couple source are often observed in microseismic monitoring, an arbitrary moment tensor source is implemented in our forward modeling tool. The increments of shear stress are equally distributed on the staggered grid to implement an accurate and symmetric moment tensor source. Our modeling tool provides an efficient way to obtain the Green’s function in anisotropic media, which is the key of anisotropic moment tensor inversion and source mechanism characterization in microseismic monitoring. In our research, wavefields in anisotropic media have been carefully simulated and analyzed in both surface array and downhole array. The variation characteristics of travel-time and amplitude of direct P- and S-wave in vertical transverse isotropic media and horizontal transverse isotropic media are distinct, thus providing a feasible way to distinguish and identify the anisotropic type of the subsurface. Analyzing the travel-times and amplitudes of the microseismic data is a feasible way to estimate the orientation and density of the induced cracks in hydraulic fracturing. Our anisotropic modeling tool can be used to generate and analyze microseismic full wavefield with full moment tensor source in anisotropic media, which can help promote the anisotropic interpretation and inversion of field data.  相似文献   

14.
—?Microseismic monitoring systems are generally installed in areas of induced seismicity caused by human activity. Induced seismicity results from changes in the state of stress which may occur as a result of excavation within the rock mass in mining (i.e., rockbursts), and changes in hydrostatic pressures and rock temperatures (e.g., during fluid injection or extraction) in oil exploitation, dam construction or fluid disposal. Microseismic monitoring systems determine event locations and important source parameters such as attenuation, seismic moment, source radius, static stress drop, peak particle velocity and seismic energy. An essential part of the operation of a microseismic monitoring system is the reliable detection of microseismic events. In the absence of reliable, automated picking techniques, operators rely upon manual picking. This is time-consuming, costly and, in the presence of background noise, very prone to error. The techniques described in this paper not only permit the reliable identification of events in cluttered signal environments they have also enabled the authors to develop reliable automated event picking procedures. This opens the way to use microseismic monitoring as a cost-effective production/operations procedure. It has been the experience of the authors that in certain noisy environments, the seismic monitoring system may trigger on and subsequently acquire substantial quantities of erroneous data, due to the high energy content of the ambient noise. Digital filtering techniques need to be applied on the microseismic data so that the ambient noise is removed and event detection simplified. The monitoring of seismic acoustic emissions is a continuous, real-time process and it is desirable to implement digital filters which can also be designed in the time domain and in real-time such as the Kalman Filter. This paper presents a real-time Kalman Filter which removes the statistically describable background noise from the recorded seismic traces.  相似文献   

15.
v--vRegional crustal waveguide calibration is essential to the retrieval of source parameters and the location of smaller (M < 4.8) seismic events. This path calibration of regional seismic phases is strongly dependent on the accuracy of hypocentral locations of calibration (or master) events. This information can be difficult to obtain, especially for smaller events. Generally, explosion or quarry blast generated travel-time data with known locations and origin times are useful for developing the path calibration parameters, but in many regions such data sets are scanty or do not exist. We present a method which is useful for regional path calibration independent of such data, i.e. with earthquakes, which is applicable for events down to Mw = 4 and which has successfully been applied in India, central Asia, western Mediterranean, North Africa, Tibet and the former Soviet Union. These studies suggest that reliably determining depth is essential to establishing accurate epicentral location and origin time for events. We find that the error in source depth does not necessarily trade-off only with the origin time for events with poor azimuthal coverage, but with the horizontal location as well, thus resulting in poor epicentral locations. For example, hypocenters for some events in central Asia were found to move from their fixed-depth locations by about 20 km. Such errors in location and depth will propagate into path calibration parameters, particularly with respect to travel times. The modeling of teleseismic depth phases (pP, sP) yields accurate depths for earthquakes down to magnitude Mw = 4.7. This Mw threshold can be lowered to four if regional seismograms are used in conjunction with a calibrated velocity structure model to determine depth, with the relative amplitude of the Pnl waves to the surface waves and the interaction of regional sPmP and pPmP phases being good indicators of event depths. We also found that for deep events a seismic phase which follows an S-wave path to the surface and becomes critical, developing a head wave by S to P conversion is also indicative of depth. The detailed characteristic of this phase is controlled by the crustal waveguide. The key to calibrating regionalized crustal velocity structure is to determine depths for a set of master events by applying the above methods and then by modeling characteristic features that are recorded on the regional waveforms. The regionalization scheme can also incorporate mixed-path crustal waveguide models for cases in which seismic waves traverse two or more distinctly different crustal structures. We also demonstrate that once depths are established, we need only two-stations travel-time data to obtain reliable epicentral locations using a new adaptive grid-search technique which yields locations similar to those determined using travel-time data from local seismic networks with better azimuthal coverage.  相似文献   

16.
v—vRegional seismic discriminants for the May 11, 1998 Indian underground nuclear test(s) and earthquakes recorded at station NIL (Nilore, Pakistan) provide new data to test strategies that can be used to monitor the Comprehensive Nuclear-Test-Ban Treaty (CTBT). Three categories of regional discriminants (ratios of P- and/or S-wave energy) were measured on short-period (0.5–6uHz) seismograms: P/S amplitude ratios (phase ratios) measured in the same frequency band, P- and S-wave spectral ratios (i.e., low frequency to high frequency for the same phase) and P/S cross-spectral ratios (i.e., low frequency S-wave to high frequency P-wave). The P/S amplitude ratios show good separation of the Indian nuclear test and regional earthquakes for Pn/Lg and Pn/Sn, however Pg/Lg does not discriminate as well. Pn/Lg and Pn/Sn discriminate well at frequencies as low as 0.5–2uHz, especially after accounting for path effects. This observation differs from previous studies that report poor separation of earthquakes and explosions at lower frequencies. The P/S amplitude ratios do not show any magnitude dependence, suggesting that forming the ratios in a fixed frequency band cancels the effects of source size-corner frequency scaling. Spatial variability of the observed discriminants arises from variations in crustal waveguide and/or attenuation structure (path propagation effects). Grouping amplitude ratios for earthquakes with paths similar to the Indian test greatly improves discrimination. Removing distance trends does not generally improve discrimination. Accounting for path effects with Bayesian kriging significantly improves discrimination. Spectral ratios (e.g., Pn [0.5–1uHz]/Pn [4–6uHz]) and cross-spectral ratios (e.g., Lg [1–2uHz]/Pn [4–6uHz]) show distance and magnitude dependence. We developed a technique for simultaneously removing the effects of distance and source size-corner frequency scaling on amplitude and spectral ratios. The technique uses a gridsearch to find several parameters that characterize the observed distance and magnitude dependence. Discrimination of the Indian test improved dramatically after the distance and magnitude trends were removed from the spectral and cross-spectral ratio data.  相似文献   

17.
张达  戴锐  曾志毅  冀虎  石雅倩  常莹  韩鹏 《中国地震》2021,37(2):332-348
深井高应力、高岩压诱发的冒顶、片帮、塌方、岩爆等地压灾害严重影响了矿山安全生产.随着高精度、高采样率传感器的研发和计算机技术的飞速发展,微震高精度实时监测技术不断完善,微震监测已成为矿山地压灾害监测预警的重要手段.通过设计、研发完全自主知识产权的BSN(BGRIMM Seismic Network)微震监测系统,实现了...  相似文献   

18.
Heavily populated by Beijing and Tianjin cities, Bohai basin is a seismically active Cenozoic basin suffering from huge lost by devastating earthquakes, such as Tangshan earthquake. The attenuation (QP and QS) of the surficial Quaternary sediment has not been studied at natural seismic frequency (1?10 Hz), which is crucial to earthquake hazards study. Borehole seismic records of micro earthquake provide us a good way to study the velocity and attenuation of the surficial structure (0?500 m). We found that there are two pulses well separated with simple waveforms on borehole seismic records from the 2006 MW4.9 Wen'an earthquake sequence. Then we performed waveform modeling with generalized ray theory (GRT) to confirm that the two pulses are direct wave and surface reflected wave, and found that the average vP and vS of the top 300 m in this region are about 1.8 km/s and 0.42 km/s, leading to high vP/vS ratio of 4.3. We also modeled surface reflected wave with propagating matrix method to constrain QS and the near surface velocity structure. Our modeling indicates that QS is at least 30, or probably up to 100, much larger than the typically assumed extremely low Q (~10), but consistent with QS modeling in Mississippi embayment. Also, the velocity gradient just beneath the free surface (0?50 m) is very large and velocity increases gradually at larger depth. Our modeling demonstrates the importance of borehole seismic records in resolving shallow velocity and attenuation structure, and hence may help in earthquake hazard simulation.  相似文献   

19.
From August 2016 to July 2017, a passive seismic survey was conducted in South Western Iran as a part of a pilot project aimed to improve the imaging in geologically complex areas. Passive seismic methods have shown to be a useful tool to infer the physical properties of the underground geological structures where traditional hydrocarbon exploration methods are challenging. For this purpose, a dense passive seismic network consisting of 119 three-component borehole seismic stations was deployed over an area of 400 km2 around the city of Dehdasht. This paper focuses on the details of the network design, which was devoted to high-resolution seismological applications, including local earthquake tomography and seismic attenuation imaging. In this regard, we describe the instrument types and the station installation procedures used to obtain high-quality data that were used to retrieve three-dimensional models of P- and S-wave velocity and P-wave attenuation in the area using tomographic inversion techniques. We also assess the network performance in terms of the seismic ambient noise levels recorded at each station site, and we revise the horizontal orientation of the sensors using surface waves from teleseismic earthquakes.  相似文献   

20.
—?In this paper, we use data from seismic stations operated by NORSAR, the Kola Regional Seismological Centre (KRSC) and IRIS to study the characteristics of regional phases in the European Arctic, with emphasis on the P/S ratio discriminant. While the detection and location capability of the regional station network is outstanding, source classification of small seismic events has proved very difficult. For example, the m b ?=?3.5 seismic event near Novaya Zemlya on 16 August, 1997 has been the subject of extensive analysis in order to locate it reliably and to classify the source type. We consider the application of the P/S discriminant in the context of this event and other events observed at regional distances in the European Arctic. We show that the P/S ratios of Novaya Zemlya nuclear explosions measured in the 1–3?Hz filter band scale with magnitude, indicating a need for caution and further research when applying P/S discriminants. Using mainly data from the large NORSAR array, we note that observed P/S amplitude ratios in the European Arctic show large variability for the same source type and similar propagation paths, even when considering closely spaced observation points. This effect is most pronounced at far regional distances and relatively low frequencies (typically 1–3?Hz), but it is also significant on closer recordings (around 10 degrees) and at higher frequencies (up to about 8?Hz). Our conclusion from this study is that the P/S ratio at high frequencies (e.g., 6–8?Hz) shows promise as a discriminant between low-magnitude earthquakes and explosions in the European Arctic, but its application will require further research, including extensive regional calibration and detailed station-source corrections. Such research should also focus on combining the P/S ratio with other short-period discriminants, such as complexity and spectral ratios.  相似文献   

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

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