期刊文献+
共找到14篇文章
< 1 >
每页显示 20 50 100
An improved rotated staggered grid fi nite difference scheme in coal seam
1
作者 Li Qin Ma Sui-Bo +1 位作者 Zhao Bin Zhang Wei 《Applied Geophysics》 SCIE CSCD 2020年第5期890-898,905,共10页
A coal seam is thin compared to the wavelength of seismic waves and usually shows strong anisotropy.It may form special geological formations,such as goafs and collapses,in coal mines.The existence of these formations... A coal seam is thin compared to the wavelength of seismic waves and usually shows strong anisotropy.It may form special geological formations,such as goafs and collapses,in coal mines.The existence of these formations may lead to instability in numerical simulations of the goaf area in a coal seam.The calculation speed of simulations is always a factor that restricts the development of simulation techniques.To improve the accuracy and effi ciency of seismic numerical simulations of goaf areas,an improved vacuum method has been incorporated into a rotated staggered grid scheme and calculations implemented by combining parallel computing and task parallelism.This ensures that the proposed numerical simulation method can be utilized in a geological model with large differences in elastic parameters among layers and improve the performance of a parallel application by enabling the full use of processor resources to expedite the calculations.We set up anisotropic coal seam models and then analyze numerically the characteristics of synthetic seismograms and snapshots of diff erent goaf areas with or without collapse.The results show that the proposed method can accurately simulate the goaf area and the calculation method can run with a high speed and parallel efficiency.The research will further advance the technology of anisotropic seismic exploration in coal fi elds,provide data for seismic inversion and build a theoretical support for coal mine disaster prediction. 展开更多
关键词 coal seam ANISOTROPY fi nite diff erence goaf area SIMULATION
在线阅读 下载PDF
Time difference correlation between earthquake lights and seismic ground accelerations
2
作者 Juan Antonio Lira Maurizio Mulas 《Applied Geophysics》 SCIE CSCD 2021年第1期9-16,128,共9页
Although earthquake lights have been known since ancient times,it has not been easy to study them.It was not until the 60s that the first photographs of them were taken.During the Peruvian earthquake in 2007,it was po... Although earthquake lights have been known since ancient times,it has not been easy to study them.It was not until the 60s that the first photographs of them were taken.During the Peruvian earthquake in 2007,it was possible to obtain the fi rst fi lm recording on earthquake lights.Likewise,during the earthquakes in Ecuador in 2016 and in Mexico in 2017,two fi lms of the earthquake lights were recorded.These fi lm recordings have helped in the study of earthquake lights,both for their objectivity and for their informational content.Several causal mechanisms have been proposed to explain earthquake lights:piezoelectricity,radon emanation,fluid diffusion,friction-vaporization,positive holes and dipole currents,among others.In this work a time difference correlation between earthquake lights and seismic ground accelerations was found and we use both seismic data and fi lm recordings of earthquake lights to explain its origin.In the discussion section it is suggested that fracturing of rocks manifest itself to some extent in the form of static electricity producing earthquake lights through induction The induction model proposed is new and it can explain the formation of EQL,even if the earth’s crust has layers of large electrical resistivity.The model also explains the formation of seismic lights without the need for special conditions on the earth’s surface or in the atmosphere.A better understanding of the earthquake lights generation process can improve our understanding of seismicity and help in the prediction of earthquakes. 展开更多
关键词 time diff erence correlation earthquake lights seismic ground accelerations rock fracturing
在线阅读 下载PDF
Blood-brain barrier and cerebral blood flow: Age differences in hemorrhagic stroke
3
作者 Semyachkina-Glushkovskaya Oxana Sindeeva Olga +7 位作者 Abdurashitov Arkady Sindeev Sergey Zinchenko Ekaterina Gekaluk Artem Ulanova Maria Mohanad Kassim Yankovskaya Ludmila Tuchin Valery 《Journal of Innovative Optical Health Sciences》 SCIE EI CAS 2015年第6期101-110,共10页
Neonatal stroke is similar to the stroke that occurs in adults and produces a significant morbidity and long-term neurologic and cognitive deficits.There are important differences in the factors,clinical events and ou... Neonatal stroke is similar to the stroke that occurs in adults and produces a significant morbidity and long-term neurologic and cognitive deficits.There are important differences in the factors,clinical events and outcomes associated with the stroke in infants and adults.However,mechanisms underlying age differences in the stroke development remain largely unknown.Therefore,treatment guidelines for neonatal stroke must extrapolate from the adult data that is often not suitable for children.The new information about differences between neonatal and adult stroke is essential for identification of significant areas for future treatment and effective prevention of neonatal stroke.Here,we studied the development of stress-induced hemorrhagic stroke and possible mechanisms underlying these processes in newborn and adult rats.Using histological methods and magnetic resonance imaging,we found age differences in the type of intracranial hemorrhages.Newborn rats demonstrated small superficial bleedings in the cortex while adult rats had more severe deep bleedings in the cerebellum.Using Doppler optical coherent tomography,we found higher stress-reactivity of the sagittal sinus to deleterious effects of stress in newborn vs.adult rats suggesting that the cerebral veins are more vulnerable to negative stress factors in neonatal vs.adult brain in rats.However,adult but not newborn rats demonstrated the stroke-induced breakdown of blood brain barrier(BBB)permeability.The one of possible mechanisms underlying the higher resistance to stress-related stroke injures of cerebral vessels in newborn rats compared with adult animals is the greater expression of two main tight junction proteins of BBB(occludin and claudin-5)in neonatal vs.mature brain in rats. 展开更多
关键词 STROKE age di®erences cerebral blood flow brain blood barrier
原文传递
A Finite-Difference Approximation for the Oneand Two-Dimensional Tempered Fractional Laplacian
4
作者 Yaoqiang Yan Weihua Deng Daxin Nie 《Communications on Applied Mathematics and Computation》 2020年第1期129-145,共17页
This paper provides a finite-difference discretization for the one-and two-dimensional tempered fractional Laplacian and solves the tempered fractional Poisson equation with homogeneous Dirichlet boundary conditions.T... This paper provides a finite-difference discretization for the one-and two-dimensional tempered fractional Laplacian and solves the tempered fractional Poisson equation with homogeneous Dirichlet boundary conditions.The main ideas are to,respectively,use linear and quadratic interpolations to approximate the singularity and non-singularity of the one-dimensional tempered fractional Laplacian and bilinear and biquadratic interpolations to the two-dimensional tempered fractional Laplacian.Then,we give the truncation errors and prove the convergence.Numerical experiments verify the convergence rates of the order O(h^2−2s). 展开更多
关键词 Tempered fractional Laplacian Finite-diff erence scheme Linear and quadratic interpolations Bilinear and biquadratic interpolations Convergence rates
在线阅读 下载PDF
Mathematical stencil and its application in finite difference approximation to the Poisson equation 被引量:4
5
作者 FENG Hui ZHANG Baoin LIU Yang 《Science China Mathematics》 SCIE 2005年第10期1421-1429,共9页
The concept of mathematical stencil and the strategy of stencil elimination for solving the finite difference equation is presented,and then a new type of the iteration algorithm is established for the Poisson equatio... The concept of mathematical stencil and the strategy of stencil elimination for solving the finite difference equation is presented,and then a new type of the iteration algorithm is established for the Poisson equation.The new algorithm has not only the obvious property of parallelism,but also faster convergence rate than that of the classical Jacobi iteration.Numerical experiments show that the time for the new algorithm is less than that of Jacobi and Gauss-Seidel methods to obtain the same precision,and the computational velocity increases obviously when the new iterative method,instead of Jacobi method,is applied to polish operation in multi-grid method,furthermore,the polynomial acceleration method is still applicable to the new iterative method. 展开更多
关键词 mathemati cal stencil stencil elimination Poisson equation fi nite dif erence iterative algorithm PARALLELISM
原文传递
Aerodynamic interaction between forewing and hindwing of a hovering dragonfly 被引量:5
6
作者 Zheng Hu Xin-Yan Deng 《Acta Mechanica Sinica》 SCIE EI CAS CSCD 2014年第6期787-799,共13页
The phase change between the forewing and hindwing is a distinct feature that sets dragonfly apart from other insects.In this paper,we investigated the aerodynamic effects of varying forewing-hindwing phase di ff eren... The phase change between the forewing and hindwing is a distinct feature that sets dragonfly apart from other insects.In this paper,we investigated the aerodynamic effects of varying forewing-hindwing phase di ff erence with a60 inclined stroke plane during hovering flight.Force measurements on a pair of mechanical wing models showed that in-phase flight enhanced the forewing lift by 17%and the hindwing lift was reduced at most phase differences.The total lift of both wings was also reduced at most phase di ff erences and only increased at a phase range around in-phase.The results may explain the commonly observed behavior of the dragonfly where 0 is employed in acceleration.We further investigated the wing-wing interaction mechanism using the digital particle image velocimetry(PIV)system,and found that the forewing generated a downwash flow which is responsible for the lift reduction on the hindwing.On the other hand,an upwash flow resulted from the leading edge vortex of the hindwing helps to enhance lift on the forewing.The results suggest that the dragonflies alter the phase di ff erences to control timing of the occurrence of flow interactions to achieve certain aerodynamic effects. 展开更多
关键词 Dragonfly aerodynamics Forewing-hindwing interaction Phase di ff erence Interactive flow
在线阅读 下载PDF
Study on the crustal stress fi eld of the Tengchong volcanic area using composite focal mechanism method 被引量:3
7
作者 Sheng Shu-zhong Wan Yong-ge +4 位作者 Jiang Chang-sheng Wang Xiao-shan Liang Shan-shan Xiao Gen-ru Hu Xiao-hui 《Applied Geophysics》 SCIE CSCD 2021年第2期238-251,274,共15页
We calculated the crustal stress field using the composite focal mechanism method based on the P-wave initial motion polarity data of the Tengchong volcanic area from January 2011 to April 2019 obtained from the Bulle... We calculated the crustal stress field using the composite focal mechanism method based on the P-wave initial motion polarity data of the Tengchong volcanic area from January 2011 to April 2019 obtained from the Bulletin of Seismological Observations of Chinese stations.The magnitude range of earthquakes used in this study is 0–4,and their magnitudes are mainly approximately 1.0.To investigate the infl uence of the source location on the stress fi eld and obtain reliable stress fi elds of the study area,we applied the double-diff erence algorithm to relocate the seismic events,obtaining more accurate and reliable relative positions of seismic events with a clearer seismic belt.On the basis of relocation results,the study on the stress fi eld along the fault zone was conducted,and the infl uence of seismic event position on the stress fi eld was analyzed.Results show that,fi rst,the current stress regime in the shallow crust of the Tengchong volcanic area is strike-slip faulting,the orientation of the principal compressive stress axis is NE–SW,the orientation of the principal extension stress axis is SE–NW,the principal compressive and extension stress axes are nearly horizontal,and the dip angle of intermediate principal stress axis is relatively large.This reflects that the volcanic and seismic activities in the Tengchong volcanic area are mainly controlled by the collision and squeezing eff ect of the Indian–Eurasian plate.It also refl ects that the current tensile action caused by deep magma activity has little infl uence on the shallow crustal stress field.Second,the stress field along fault zones reveals that there exist local stress fi elds,such as the thrust stress regime at the strike-slip fault terminal area,which is consistent with the compressional area at the intersection of conjugate strike-slip faults indicated by previous study.Third,the stress fi eld results are consistent,regardless of using the original location in the bulletin or the relocated location,indicating that the infl uence of the event location error can be neglected when there are suffi cient data and refl ecting the stability of the composite focal mechanism method.The findings can serve as a reference for investigating geological structure movement,seismic activities,and volcanic activities in the Tengchong volcanic area. 展开更多
关键词 Tengchong volcanic area Composite focal mechanism Stress fi eld Double-diff erence location Conjugate Fault
在线阅读 下载PDF
Smart Semi-active PID-ACO control strategy for tower vibration reduction in Wind Turbines with MR damper 被引量:2
8
作者 Mahmudur Rahman Zhi Chao Ong +1 位作者 Wen Tong Chong Sabariah Julai 《Earthquake Engineering and Engineering Vibration》 SCIE EI CSCD 2019年第4期887-902,共16页
Wind turbine technology is well known around the globe as an eco-friendly and eff ective renewable power source. However, this technology often faces reliability problems due to structural vibration. This study propos... Wind turbine technology is well known around the globe as an eco-friendly and eff ective renewable power source. However, this technology often faces reliability problems due to structural vibration. This study proposes a smart semi-active vibration control system using Magnetorheological (MR) dampers where feedback controllers are optimized with nature-inspired algorithms. Proportional integral derivative (PID) and Proportional integral (PI) controllers are designed to achieve the optimal desired force and current input for MR the damper. PID control parameters are optimized using an Ant colony optimization (ACO) algorithm. The eff ectiveness of the ACO algorithm is validated by comparing its performance with Ziegler-Nichols (Z-N) and particle swarm optimization (PSO). The placement of the MR damper on the tower is also investigated to ensure structural balance and optimal desired force from the MR damper. The simulation results show that the proposed semi-active PID-ACO control strategy can signifi cantly reduce vibration on the wind turbine tower under diff erent frequencies (i.e., 67%, 73%, 79% and 34.4% at 2 Hz, 3 Hz, 4.6 Hz and 6 Hz, respectively) and amplitudes (i.e. 50%, 58% and 67% for 50 N, 80 N, and 100 N, respectively). In this study, the simulation model is validated with an experimental study in terms of natural frequency, mode shape and uncontrolled response at the 1st mode. The proposed PID-ACO control strategy and optimal MR damper position is also implemented on a lab-scaled wind turbine tower model. The results show that the vibration reduction rate is 66% and 73% in the experimental and simulation study, respectively, at the 1st mode. 展开更多
关键词 ant COLONY optimization fi nite diff erence method MAGNETORHEOLOGICAL damper PID-ACO controller semiactive control system identifi cation wind turbine TOWER
在线阅读 下载PDF
Three-dimensional high-resolution velocity structure imaging and seismicity study of Yangbi Ms6.4 earthquake 被引量:2
9
作者 Ma Yong Zhang Hai-Jiang +1 位作者 Gao Lei Bi Jin-Meng 《Applied Geophysics》 SCIE CSCD 2021年第4期579-591,595,共14页
In this study,based on the body wave arrival data of 5506 earthquakes recorded by 32 fi xed stations and 94 temporary stations in Yangbi and surrounding areas,the source parameters of Yangbi Ms6.4 earthquake sequence ... In this study,based on the body wave arrival data of 5506 earthquakes recorded by 32 fi xed stations and 94 temporary stations in Yangbi and surrounding areas,the source parameters of Yangbi Ms6.4 earthquake sequence and three-dimensional(3-D)fi ne Vp,Vs,and Vp/Vs were inverted by using the consistency-constrained double-diff erence tomography method.The results showed that the focal depth after relocation was mostly in the range of 3–10 km,evidently nearly horizontally distributed,and concentrated in the weak area of the high-velocity body or at the side of the high-low-velocity body transition zone toward the high-velocity body,showing a good corresponding relationship with the velocity structure.The velocity structure in the Yangbi area has remarkably uneven characteristics.The seismic activity area is dominated by high-velocity bodies prone to brittle fracture near the surface.As the depth increases,low-velocity anomalies appear.A signifi cant diff erence was observed in the wave velocity ratio between the upper and lower sides of the seismically dense strip.Based on the focal mechanism of the Yangbi Ms6.4 earthquake and the fine 3-D velocity structure,this article concludes that the Yangbi Ms6.4 earthquake was caused by a strong regional tectonic stress concentrated in the relatively weak area by hard high-velocity bodies on the northwest sides.The Ms5.6 foreshock broke the inherent balance of regional stress and promoted the occurrence of the Yangbi Ms6.4 mainshock.Afterward,the stress was adjusted to a new equilibrium state through a large number of aftershocks,forming a foreshock–mainshock–aftershock type of seismic activity model.Based on the activity law of the Yangbi Ms6.4 earthquake sequence and characteristics of the 3D velocity structure distribution,this paper speculates that the seismogenic structure of the Yangbi earthquake was possibly a northwest strike-slip buried fault with a depth of 3–10 km on the southwest side of the Weixi–Qiaohou fault. 展开更多
关键词 Tangshan fault double-diff erence tomography velocity structure seismic activity Vp/Vs inversion
在线阅读 下载PDF
Q estimation based on logarithmic spectral areas with different high and low frequencies 被引量:2
10
作者 An Yong Wang Xin-You 《Applied Geophysics》 SCIE CSCD 2021年第1期75-84,129,130,共12页
When a seismic wave propagates through subsurface viscoelastic media,the formation absorbs the high-frequency energy of the seismic wave more strongly than the low-frequency energy.As the depth and the off set increas... When a seismic wave propagates through subsurface viscoelastic media,the formation absorbs the high-frequency energy of the seismic wave more strongly than the low-frequency energy.As the depth and the off set increase,the diff erence between the logarithmic spectral areas with high and low frequencies gradually increases.Based on this seismic wave characteristic,we have developed a novel Q-estimation method based on logarithmic spectral area diff erence of high and low frequencies(referred to as the LSAD_LH method).In this paper,we derive the theoretical relationship between the Q value and difference of logarithmic spectral areas with high and low frequencies and prove the applicability of the LSAD_LH method using a single-layer medium numerical model.To verify the sensitivity of the LSAD_LH method to bandwidth selection and noise,we compare the LSAD_LH method with two credible methods—the logarithmic spectral ratio(LSR)and logarithmic spectral area diff erence(LSAD)methods using a synthetic model containing the random noise.The results demonstrate that the LSAD_LH method is not highly dependent on bandwidth,and in terms of noise immunity,it is signifi cantly better than the LSR method and has the same advantages as the LSAD method.To further highlight the advantages of the LSAD_LH method,we apply the LSAD_LH and LSAD methods to the vertical seismic profi ling(VSP)numerical simulation of the multilayer media and the fi eld zero-off set VSP data.The application of the two cases proves the applicability of the LSAD_LH method and the accuracy of the high Q-value estimation relative to the LSAD method.Moreover,via the transmission coeffi cient,the LSAD_LH method overcomes the weakness of the LSAD method. 展开更多
关键词 Attenuation of high-and low-frequency components Quality factor Logarithmic spectral area diff erence
在线阅读 下载PDF
Multi-temporal NDVI analysis using UAV images of tree crowns in a northern Mexican pine-oak forest 被引量:1
11
作者 JoséLuis Gallardo-Salazar Marcela Rosas-Chavoya +4 位作者 Marín Pompa-García Pablito Marcelo López-Serrano Emily García-Montiel Arnulfo Meléndez-Soto Sergio Iván Jiménez-Jiménez 《Journal of Forestry Research》 SCIE CAS CSCD 2023年第6期1855-1867,共13页
The use of unmanned aerial vehicles(UAV)for forest monitoring has grown significantly in recent years,providing information with high spatial resolution and temporal versatility.UAV with multispectral sensors allow th... The use of unmanned aerial vehicles(UAV)for forest monitoring has grown significantly in recent years,providing information with high spatial resolution and temporal versatility.UAV with multispectral sensors allow the use of indexes such as the normalized difference vegetation index(NDVI),which determines the vigor,physiological stress and photo synthetic activity of vegetation.This study aimed to analyze the spectral responses and variations of NDVI in tree crowns,as well as their correlation with climatic factors over the course of one year.The study area encompassed a 1.6-ha site in Durango,Mexico,where Pinus cembroides,Pinus engelmannii,and Quercus grisea coexist.Multispectral images were acquired with UAV and information on meteorological variables was obtained from NASA/POWER database.An ANOVA explored possible differences in NDVI among the three species.Pearson correlation was performed to identify the linear relationship between NDVI and meteorological variables.Significant differences in NDVI values were found at the genus level(Pinus and Quercus),possibly related to the physiological features of the species and their phenology.Quercus grisea had the lowest NDVI values throughout the year which may be attributed to its sensitivity to relative humidity and temperatures.Although the use of UAV with a multispectral sensor for NDVI monitoring allowed genera differentiation,in more complex forest analyses hyperspectral and LiDAR sensors should be integrated,as well other vegetation indexes be considered. 展开更多
关键词 Multispectral images Normalized diff erence Vegetation index PHENOLOGY Unmanned aerial vehicles Multitemporal analysis
在线阅读 下载PDF
Phylogenetic relationship of Picea mongolica with other Picea species in the same area based on chloroplast gene variations 被引量:1
12
作者 Hao Di Jingjing Ma +3 位作者 Kuaikuai He Fangxu Han Yue Li Shihui Niu 《Journal of Forestry Research》 SCIE CAS CSCD 2021年第1期297-305,共9页
Picea mongolica is a conifer with a limited distribution in China and its taxonomic status is controversial.In order to explore genetic diff erences between P.mongolica and other nearby Picea species and to investigat... Picea mongolica is a conifer with a limited distribution in China and its taxonomic status is controversial.In order to explore genetic diff erences between P.mongolica and other nearby Picea species and to investigate its taxonomic status,phylogenetic relationships were analyzed between P.mongolica and Picea koraiensis,Picea meyeri and Picea wilsonii by three chloroplast gene sequences matK,chlB and atpA.The length of joint chloroplast sequence is 2379 bp.The fi fteen haplotypes were identifi ed by haplotype network analysis,among which two were major haplotypes and nine were unique.In addition,the genetic diversity of the sample collection species was inferred.Based on the haplotype network and Neighbor Joining phylogenetic tree analysis,P.mongolica was located on the basal clade of the phylogenetic tree which had more primitive taxa,and the genetic diversity of P.mongolica was highest.The signifi cant diff erences between P.mongolica and these other Picea species were identifi ed by this research. 展开更多
关键词 CPDNA Genetic diff erence PHYLOGENY Picea mongolica Taxonomic status
在线阅读 下载PDF
Role of Microalloyed Y and Gd in Improving the Corrosion Resistance of Rolled Mg-3Al-1Zn Alloy
13
作者 Panpan Wang Haitao Jiang +2 位作者 Yujiao Wang Yun Zhang Jingchao Tao 《Acta Metallurgica Sinica(English Letters)》 SCIE EI CAS CSCD 2022年第6期941-947,共7页
The corrosion performances of rolled AZ31 alloy with addition of Y and Gd were comparatively investigated. The corrosion rates of AZ31 alloy modified with Y and Gd were 3.3 and 3.7 mm/y immersing for 24 h in 3.5 wt% N... The corrosion performances of rolled AZ31 alloy with addition of Y and Gd were comparatively investigated. The corrosion rates of AZ31 alloy modified with Y and Gd were 3.3 and 3.7 mm/y immersing for 24 h in 3.5 wt% NaCl solution, respectively, which were much lower than that of AZ31 alloy(23.6 mm/y). The remarkable improvement in corrosion performance by adding Y or Gd was ascribed to preferentially form less noble Al_(2)Y or Al_(6)GdMn_(6) and more protective corrosion product layer, resulting in the reduction in H_(2) evolution rate and the enhancement of passivity. 展开更多
关键词 AZ31 alloy Second phases Volta potential diff erence Y GD
原文传递
3D prestack reverse time migration of ground penetrating radar data based on the normalized correlation imaging condition
14
作者 Wang Hong-Hua Gong Jun-bo +4 位作者 Zhang Zhi Xiong Bin Lv Yu-zeng Feng De-shan Dai Qian-wei 《Applied Geophysics》 SCIE CSCD 2020年第5期709-718,901,共11页
The reverse time migration(RTM)of ground penetrating radar(GPR)is usually implemented in its two-dimensional(2D)form,due to huge computational cost.However,2D RTM algorithm is difficult to focus the scattering signal ... The reverse time migration(RTM)of ground penetrating radar(GPR)is usually implemented in its two-dimensional(2D)form,due to huge computational cost.However,2D RTM algorithm is difficult to focus the scattering signal and produce a high precision subsurface image when the object is buried in a complicated subsurface environment.To better handle the multi-off set GPR data,we propose a three-dimensional(3D)prestack RTM algorithm.The high-order fi nite diff erence time domian(FDTD)method,with the accuracy of eighth-order in space and second-order in time,is applied to simulate the forward and backward extrapolation electromagnetic fi elds.In addition,we use the normalized correlation imaging condition to obtain pre-stack RTM result and the Laplace fi lter to suppress the low frequency noise generated during the correlation process.The numerical test of 3D simulated GPR data demonstrated that 3D RTM image shows excellent coincidence with the true model.Compared with 2D RTM image,the 3D RTM image can more clearly and accurately refl ect the 3D spatial distribution of the target,and the resolution of the imaging results is far better.Furthermore,the application of observed GPR data further validates the eff ectiveness of the proposed 3D GPR RTM algorithm,and its fi nal image can more reliably guide the subsequent interpretation. 展开更多
关键词 Ground Penetrating Radar(GPR) 3D Reverse Time Migration(RTM) Finite Diff erence Time Domain(FDTD) Normalized correlation imaging condition
在线阅读 下载PDF
上一页 1 下一页 到第
使用帮助 返回顶部