80 results about "P-wave" patented technology

The invention discloses a multi-wave chromatography and reflection combined detection method and device for a thin-covering-layer urban active fault and a medium. The method comprises the steps that longitudinal wave seismic data and transverse wave seismic data are collected; recording vertical component data based on the longitudinal wave seismic data for first arrival travel time pickup and tomography inversion to obtain a longitudinal wave velocity structure model Vp; extracting Rayleigh surface wave frequency dispersion through the vertical component data and carrying out multichannel surface wave analysis inversion to obtain a transverse wave velocity structure model Vs; calculating a longitudinal and transverse wave velocity ratio parameter structure model Vp / Vs, and calculating to obtain a Poisson's ratio parameter structure model v; obtaining a longitudinal wave (P wave) reflection profile and a transverse wave (S wave) reflection profile; identifying and summarizing into an interpretation fracture set Fc, and recording key elements; and by taking the interpretation fracture set Fc as a structural skeleton, completing unified registration and superposition display on Vp, Vs, Vp / Vs and v, a longitudinal wave reflection profile and a transverse wave reflection profile, and outputting a shallow underground structure comprehensive interpretation map of the thin covering layer region.

The invention discloses an ultra-deep reflection imaging method based on passive source seismic data and a related device, and relates to the technical field of earth depth structure imaging, and the method comprises the following steps: firstly, for original continuous seismic waveform data, intercepting seismic waveform data with a first arrival of a P wave and a non-arrival of an S wave, and carrying out the preprocessing of the data, so as to obtain a far earthquake event data segment; the method comprises the following steps: firstly, analyzing and identifying a far earthquake event with a high signal-to-noise ratio by using a time-frequency-wavenumber domain, removing data which are greatly influenced by noise, and then, reconstructing a reflection wave field of a high-quality far earthquake data fragment by using a multi-dimensional deconvolution technology, and carrying out convolution on a seismic source wavelet to obtain ultra-deep reflection wave field data. And finally, performing reflection imaging processing in combination with velocity analysis and superposition imaging technologies, optimizing the analysis capability of a deep reflection interface, and finally obtaining an earth deep reflection imaging profile, so that a continuous structure interface extending from the shallow earth crust to the deep mantle can be clearly revealed, and coupling characteristics between shallow and deep earth structures can be effectively captured.

The disclosure provides a method and device for velocity tomography imaging of seabed shallow media, and electronic equipment. The method includes performing multichannel analyses of Guided-P and Scholte waves on seabed multi-component seismic gathers to determine measured multi-order dispersion curves of the two waves; based on a first seabed shallow media model and a theoretical dispersion equation of Guided-P wave, performing joint inversion of multi-order dispersion curves of Guided-P wave to iteratively update the P-wave velocity of seabed media; based on a second seabed shallow media model and a theoretical dispersion equation of Scholte wave, performing joint inversion of multi-order dispersion curves of Scholte wave to iteratively update the S-wave velocity of seabed media under the constraint of the P-wave velocity determined by the Guided-P wave dispersion inversion; and performing tomography imaging of seabed velocity structures with the inverted P-wave and S-wave velocities along a survey line.

According to the atrial fibrillation recognition method based on multi-feature fusion, multiple feature extraction task branches are constructed in the same deep learning model, the RR interval features, the P-wave features and the electrocardio waveform morphological features are extracted synchronously and parallelly, and the problems of feature loss and time sequence mismatching caused by multi-step processing are avoided; the model captures RR interval irregularity in a time domain, captures P wave missing and baseline f wave features in a morphological domain, forms a rhythm-morphological-frequency spectrum triple feature system, realizes more stable identification of atrial fibrillation, can effectively distinguish similar expressions of atrial fibrillation and ventricular premature, atrial premature, sinus arrhythmia and the like based on a coordination mechanism of appetite and morphological domain, and improves the accuracy of atrial fibrillation identification. And the identification rate of the atrial fibrillation signal in the electrocardiogram data can be effectively improved.

The invention discloses a multi-stage grouting reinforcement method for water-rich broken surrounding rock, which comprises the following steps: grading broken rock mass according to the breaking degree and the water content condition of disturbed rock mass and the result of a rock mass water content correction P wave test tested by a resistivity method, flexibly adjusting the type of grouting materials according to grading, and matching different grouting processes. While the grouting material is effectively utilized and the engineering budget is reduced, the reinforcing and plugging quality of multiple areas of the fractured rock mass can be guaranteed, the safe operation of the underground engineering is guaranteed, the reinforcing and plugging quality of the fractured rock mass is guaranteed, and the safe operation of the underground engineering is guaranteed. Rock masses with different crushing degrees can obtain the best reinforcing effect through the graded grouting strategy, the overall reinforcing strength is improved, and the anti-permeability performance and the reinforcing strength of a reinforcing body are greatly improved through application of the high polymer grouting material.

The invention provides a micro-seismic P-wave first arrival pickup method based on deep learning and an adaptive time window, and the method comprises the steps: (1) collecting and preprocessing a micro-seismic signal, making a data set, and dividing the data set into a training set and a test set; (2) building an MSC-NET deep learning neural network, taking a data set as input for training, and adjusting model parameters to achieve an optimal model effect; (3) performing a micro-seismic signal classification task by applying the trained deep learning model, and dividing a large number of micro-seismic signals into an effective available class and an ineffective unavailable class; and (4) finally, carrying out P-wave first arrival time pickup on the effective micro-seismic signals based on an energy ratio algorithm of an adaptive time window to obtain a relatively reliable pickup result. According to the method, the deep learning technology and the improved traditional signal processing technology are combined, the reliability and accuracy of the automatic pickup result of the micro-seismic signal P wave first arrival time are greatly improved, and meanwhile the method has the advantages of being high in robustness and applicability and the like.

The invention discloses a reservoir water fluctuation propagation direct numerical simulation method and device in gravity dam anti-seismic analysis, a medium and equipment. According to the method, fluctuation propagation in reservoir water under combined incidence of P waves and SV waves is considered, a reservoir water-foundation site free field is analyzed and calculated on the basis of the fluctuation theory and the potential function theory, then the site free field is converted into node counterforce on fluid and solid boundaries on the basis of a fluid and solid power artificial boundary substructure model, and the reservoir water-foundation site free field is obtained. Applying node counterforce on the reservoir water-foundation artificial boundary, and directly performing numerical simulation on reservoir water-foundation site free field response; based on a static-dynamic artificial boundary conversion method, considering the influence of the static effect of the gravity dam-reservoir water-foundation system on the dynamic effect of the system, and constructing a gravity dam seismic oscillation input model for directly simulating reservoir water fluctuation propagation; according to the method, the anti-seismic analysis precision of the gravity dam is improved, and the limitations of single-wave incidence during reservoir water fluctuation propagation simulation, failure to consider the static effect and the like are overcome.

The invention discloses a magnetocardiogram P-wave detection method based on wavelet analysis and a local noise self-adaptive threshold, which is characterized in that accurate and robust detection of a P wave is realized by searching on a characteristic scale which can highlight the P wave, a self-adaptive threshold confirmation mechanism based on local noise estimation is introduced, and a self-adaptive threshold confirmation mechanism is established for a time domain search window of each cardiac cycle. According to the method, the internal noise level is dynamically evaluated, and an estimated value capable of accurately reflecting the local noise condition of the current heart beat can be obtained by carrying out statistical analysis on wavelet coefficients in a characteristic scale in an area where P waves are most likely to appear.

The invention relates to the technical field of civil engineering, in particular to a method for calculating the seismic response of a pile group under the action of seismic P waves, a pile-soil coupling model is established based on a continuous elastic dynamics theory, and pile-pile dynamic interaction factors are calculated by decoupling displacement of a free field and a scattered field of a soil body and combining coordinate transformation, so that the seismic response of the pile group is calculated. The displacement response of a pile group system is solved according to the superposition principle, the influence of pile spacing and soil mass parameters on the anti-seismic performance of the pile group is quantitatively evaluated by adopting a method of combining frequency domain amplitude-frequency characteristics and time domain PGA attenuation analysis, and theoretical support is provided for pile foundation optimization design.

The present application relates to oil and gas geophysical exploration technology field, more particularly, it relates to a kind of P-S wave data matching method based on wavelet compression, it includes the following steps: S1: PS wave seismic data is converted to compress to PP time domain;Extract seismic PP wavelet, PS wavelet, PP wavelet spectrum and PS wavelet spectrum;S2: introduce frequency variable scale conversion factor, utilize Fourier transform method to carry out wavelet compression to PS wavelet;S3: obtain filter according to PS wavelet and new PS wavelet;S4: utilize filter to carry out filter processing to PS wave seismic data after time domain conversion;S5: carry out waveform matching.The present application is based on the pre-matching in time domain, introduce frequency variable scale conversion factor, utilize the Fourier scale conversion method of frequency variable to carry out compression to PS wavelet.Then filter is designed using regularization method to carry out filter processing to original PS wave seismic record, can improve the similarity of PP, PS wave both seismic record waveform, improve the precision of final matching.

The invention discloses an inversion method based on three-dimensional reflection P wave monoclinic medium anisotropy parameters. The invention provides the three-dimensional reflection seismic wave travel time inversion method which is higher in stability, less in dependence on prior information and more sufficient in excavation of underground medium information. In the inversion process, firstly, dependence on prior information such as stratum thickness and vertical propagation speed is reduced, and meanwhile, more anisotropic parameters and relatively accurate stratum thickness can be obtained through inversion. Therefore, three-dimensional reflection seismic wave data can be explored more fully, richer underground medium information is provided, and a foundation is laid for a subsequent oil-gas exploration process, especially an unconventional oil-gas exploration process. A geophysical data processor can process the multi-azimuth single-layer horizontal plane reflection P wave travel time according to the method, anisotropy information of a bottom layer is obtained, or a multi-layer reflection P wave travel time and converted wave travel time inversion method is further developed according to the method provided by the invention.

The application relates to the technical field of earthquake prediction, and particularly provides a lightweight on-site ground motion parameter prediction model construction method and an earthquake early warning system. The lightweight on-site ground motion parameter prediction model construction method comprises the following steps: training an initial prediction model according to acceleration waveform data when a P wave of an earthquake arrives, to obtain a pre-training prediction model; taking the pre-training prediction model as a teacher model, and transferring the knowledge of the teacher model to a student model through distillation; and constructing a lightweight on-site ground motion parameter prediction model based on the student model. The method can reduce the complexity and calculation amount of the student model by transferring the knowledge in the pre-training prediction model to the student model through distillation; and the lightweight on-site ground motion parameter prediction model can be constructed based on the student model, so that the demand of the model for hardware resources can be reduced while the on-site ground motion parameter prediction efficiency is ensured, and the deployable scene of the model is expanded.

The invention discloses an electrocardiogram ventricular escape real-time judgment method and system, relates to the technical field of electrocardiograms, solves the problem of lack of multi-dimensional analysis and dynamic characteristic verification of electrocardiogram waveforms, and aims at solving the problem that a time period in which waveform abnormity really exists is focused through continuous analysis of fluctuation points of a to-be-recognized section, so that the accuracy of the electrocardiogram ventricular escape real-time judgment is improved. A symmetric space vector range of a normal P wave is used as a reference, whether a to-be-analyzed waveform contains the P wave is visually judged through matching of the form and the position of a to-be-verified vector, missed judgment caused by low amplitude of the P wave or overlapping of the P wave in a T wave is avoided, and the accuracy of the judgment is improved. A complete analysis chain from electrocardiosignals to clinical diagnosis is constructed through a four-layer logic architecture of pacing point positioning, feature waveform matching, abnormal time capturing and escape pacing mechanism verification, high-precision recognition of ventricular escape pacing is achieved, and the accuracy of ventricular escape pacing recognition is improved through deep fusion of a physiological mechanism and an algorithm. And an interpretable and traceable technical path is provided for automatic electrocardio diagnosis.

The utility model relates to the technical field of geothermal wells, in particular to a damping ring for damping of a geothermal well, which comprises a circular annular shell, a plurality of damping mechanisms are arranged in the annular shell, and the damping mechanisms are uniformly distributed in the annular shell and are symmetrically arranged in pairs along the radial direction of the annular shell; the damping mechanism comprises a plurality of longitudinal damping assemblies which are longitudinally arranged and a plurality of transverse damping assemblies which are horizontally arranged, the longitudinal damping assemblies are used for conducting damping and energy dissipation on vibration generated by P waves in seismic waves, and the transverse damping assemblies are used for conducting damping and energy dissipation on vibration generated by S waves in the seismic waves. Vibration caused by P waves and S waves in an earthquake can be subjected to vibration reduction and energy dissipation treatment, and the anti-seismic performance of the geothermal well is effectively improved.

A smart ring (104) for detecting atrial fibrillation using photoplethysmography is disclosed The smart ring (104) includes a PPG sensor for gathering PPG signals corresponding to the blood flow data of the user (102). A user device (106) receives the blood flow data and extracts a plurality of features from the blood flow data including time between heartbeats (RR intervals), variability in RR intervals, absence of P waves, presence of rapid and irregular heartbeat, SpO2, amplitude and width of the peaks, harmonic components of PPG, pulse wave velocity, etc. A weightage is assigned to each of the features. Further, the user device (106) implements a machine learning model (166) to detect a probability of occurrence of AFib based on the blood flow data. The user device (106) may further notify the user (102) about the probability of occurrence of AFib.

Techniques are disclosed for segmenting electrocardiogram (ECG) signals. In one example, a method to segment an electrocardiogram (ECG) signal may include detecting consecutive heartbeats in an ECG signal. The method also includes segmenting the ECG signal into multiple ECG segments surrounding the detected consecutive heartbeats and generating an ECG data set by joining consecutive ECG segments. The generated the ECG data set represents the detected heartbeats. In some such examples, each ECG segment is of a duration to include a QRS complex, a P wave, and a T wave.

The invention discloses a test device and method for simulating a coupling load of local explosion and four-way hysteresis on a component or a node. Comprising a vertical axial force loading assembly, a transverse pneumatic impact assembly, a double-spherical-hinge MTS actuator, a clamping assembly and a test component arranged in the vertical direction, and the clamping assembly comprises counter-force frames located on the two sides and the top of the test component correspondingly. A conventional axial load within the service period of a component can be accurately applied through the vertical axial force loading assembly, the double-spherical-hinge MTS actuator can alternately simulate P waves and Love surface waves after the P waves in an earthquake, and the transverse pneumatic impact assembly can apply a local explosion load at any time period under the action of a vertical axial force and a hysteretic load. According to the test device, the test working condition of coupling of local explosion and four-way hysteresis load on a test piece or a node can be truly simulated.

The application discloses a kind of based on transmission tower anti-seismic technology transmission tower anti-seismic device.The present application includes work platform and its upper surface swing component, the bottom end of the swing component is fixedly connected with bearing component, the top of the swing component is fixedly connected with pedestal;Including gravity center component, the gravity center component is located at the bottom of work platform and is fixedly connected with the lower surface of work platform;Including support component, including electric energy conversion module, calibration module and electromagnetic induction module.The present application is reduced by setting swing component, the degree of transmission tower gravity center deviation, improve the service life of transmission tower;By setting gravity center component, by reducing the gravity center of transmission tower when earthquake, to further control the degree of gravity center deviation of transmission tower;By setting support component, the impact caused by P wave to transmission tower in earthquake can be effectively absorbed and converted into the support force applied to transmission tower, so as to improve the anti-seismic strength of transmission tower.

A computer device for processing cardiac signals comprises: - a data storage (114) arranged to receive input data sets each comprising a plurality of P wave segments each associated with an electrocardiogram track and with an acquisition time window, and a plurality of coronary sinus signals associated with the same acquisition time window and having one or more activation sequence, - an input module comprising a plurality of classifiers, each associated with a respective atrial tachycardia class and arranged to receive an input data set, and to return a value indicating the probability that the input data set indicates an atrial tachycardia which belongs to the associated atrial tachycardia class, each classifier being further associated with a respective confusion matrix, said atrial tachycardia classes being arranged into a hierarchy made of layers each comprising one more node, each node of a given layer being connected to at least one other node in the hierarchy, only part in the layer immediately above and / or in the layer immediately below said given layer in said hierarchy, - a Bayesian engine comprising a prevalence matrix for pair of classifiers which atrial tachycardia classes are connected in said atrial tachycardia class hierarchy, arranged to receive a set of assessments comprising an atrial tachycardia class probability for each classifier, and to process said set of assessments according to the confusion matrices associated with the classifiers and according to said prevalence matrices in order to determine an updated set of assessments which maximizes the expression (I), where Z i are the updated set of assessments values, Y i are the set of assessments, YC i represents the states of the child nodes of the classifier associated with Y i in said atrial tachycardia class hierarchy, and N is the total number of classifiers, said device being arranged to receive an input data set, to process it with the input module, to derive therefrom a set of assessments provided to said Bayesian engine, and to return one or more atrial tachycardia information derived from said updated set of assessments.

To provide an ultra-high-speed radio communication method that can contribute to achieving higher-speed radio wave communication, which is different from a conventional communication method using a radio wave signal in which it has been considered that radio wave communication can be speeded up by shortening wavelengths of transverse waves because radio waves have been regarded as the transverse waves and which has achieved speed-up of the radio wave communication up to now, where the speed-up of communication enables more information to be communicated in a specific time, which also can enables communication volumes to increase.SOLUTION: Wave energy (a wave phenomenon) such as radio waves is divided into transverse waves that oscillate in a direction perpendicular to a transmitting direction of radio waves and longitudinal waves that oscillate in a direction parallel to the transmitting direction of the waves. The longitudinal waves are higher in transmitting speed than the transverse waves. (The same is true of longitudinal waves (P waves) and transverse waves (S waves) in an earthquake waves.) Although the radio waves have been regarded as transverse waves (limited to transverse waves) up to now, radio waves of longitudinal waves also can be generated. Therefore, higher-speed radio wave communication can be achieved by using the radio waves of the longitudinal waves.SELECTED DRAWING: None

The invention relates to a seismic source parameter prediction method based on physical constraint and array geometric adaptation, which comprises the following steps of: determining input, output and initial prediction models, sequentially obtaining station characteristics and predicted seismic source parameters, comparing the predicted seismic source parameters with a seismic source true value to obtain regression loss (if the seismic source true value does not exist, the regression loss is not recorded); obtaining phase loss according to the predicted seismic source parameter and the seismic phase label, and then obtaining total loss; carrying out a training process on the model according to the previous steps, and storing an optimal model; inputting related data into the optimal model, and obtaining predicted seismic source parameters; the phase loss only needs the arrival time data of the P wave and the S wave, and is the signal feature which is the most easily obtained signal feature, so that the application range is relatively wide; phase loss takes the physical law of seismic wave propagation as a hard constraint to be embedded into a model optimization process, so that the interpretability of a prediction result is improved; therefore, the method is wide in application range and high in interpretability when used for predicting the seismic source parameters.

A control method for an intracardiac defibrillation catheter system is provided for performing a safer treatment by shortening the time between switching the mode to perform defibrillation and synchronization with R wave. The system includes a catheter, electrode part, power supply, input part inputting an electrocardiogram, and changeover part connected to the power supply, for switching between a first mode and second mode for electrically disconnecting and connecting the power supply and electrode part, which is connected to the power supply via the changeover part. The method includes switching the first mode to second mode with the changeover part when a cardiac potential reaches a peak of P wave and / or Q wave obtained from an electrocardiographic waveform of the electrocardiogram, and applying a voltage to the electrode part with the power supply when the cardiac potential reaches a peak of R wave obtained from the electrocardiographic waveform of the electrocardiogram.

The invention provides an arrhythmia classification method based on dynamic cardiac beat division and feature fusion, and the method comprises the steps: carrying out the preprocessing operation of an original electrocardiosignal, and obtaining a preprocessed signal; performing dynamic threshold screening on the preprocessed signal to obtain a P-wave peak; carrying out waveform sequence interception on the P wave peak through a cardiac beat window to obtain cardiac beat data; and performing feature fusion on the heart beat data through a convolutional neural network to obtain a classification result. According to the method, the P wave is positioned through the intelligent algorithm, and dynamic cardiac beat division is performed in combination with the RR interval change rate features and the neural network, so that the detection precision, the feature fusion effect and the classification accuracy are improved.

The embodiment of the invention provides a high-precision seismic elastic wave field forward modeling method, device and equipment and a storage medium, and the method comprises the steps: carrying out the discretization of a wave equation in an elastic medium through the combined application of an improved time high-order implicit staggered grid difference method, and obtaining two types of time-space domain frequency dispersion relationships related to a P wave and an S wave based on the plane wave arrangement and simplification; solving the two types of time-space domain dispersion relations of the P wave and the S wave by combining a least square method and a Taylor series expansion method to obtain high-order difference coefficients corresponding to the P wave and the S wave; based on a wave field separation technology, a P wave field component is calculated by using a high-order difference coefficient corresponding to a P wave, an S wave field component is calculated by using a high-order difference coefficient corresponding to an S wave, and the components are superposed to obtain a high-precision full-elastic wave field. According to the technical scheme, under the same parameter condition, better simulation precision can be generated, so that an effective wave field continuation scheme is provided for subsequent elastic wave migration imaging and waveform inversion.

The present application relates to a seismic phase picking method based on time-frequency analysis and deep learning, belonging to the field of seismology, a deep neural network comprising an empirical mode decomposition module, a residual shrinkage module, an LSTM layer, a transformer layer, an up-sampling layer and a convolution layer is established, a linear rectifier function is used as an activation function, an Adam adaptive optimization function is used for back propagation, a binary cross-entropy loss function is used as a loss function, and the average absolute error function, error standard deviation, precision, recall and F1 score are used as evaluation indexes to obtain a deep neural network model. When there are differences in the frequencies of noise, P waves and S waves, empirical mode decomposition can effectively separate the three, thereby helping the network to better identify P waves and S waves. Even for wideband noise, the noise will be decomposed into different intrinsic mode functions, thereby reducing the noise content in the intrinsic mode function where the P wave and S wave are located. This time-frequency analysis method can effectively improve the noise immunity of the network.

The application discloses a kind of waveform and P wave polarity combined inversion method and system of collapse earthquake, comprising the following steps: S1, the waveform data of target area collapse earthquake event is obtained, and P wave first arrival polarity observation value is extracted from waveform data;S2, waveform data is preprocessed;S3, based on the collapse seismic source physical model of the compound of preset force couple point source and single force point source, the joint inversion operator of waveform forward modeling and P wave polarity forward modeling is established;S4, construct joint inversion objective function;The beneficial effects of the application are as follows: the application is based on the seismic source model of the compound of force couple and single force, constructs waveform and P wave polarity joint inversion operator;Create joint objective function including waveform fitting term and polarity fitting term, minimize objective function by optimization algorithm, and obtain source mechanism solution and source time function by joint inversion;By constructing the seismic source model of the compound of force couple and single force, the complex physical process of collapse earthquake is more accurately described.

The invention discloses a ship sound vibration characteristic identification and positioning method based on a DAS system. The method comprises the following steps: S1, obtaining monitoring data based on the DAS system; s2, classifying data based on two modes of frequency domain features and differential distance, and extracting underwater acoustic signals and vibration signals; s3, performing feature extraction and beam forming on the underwater acoustic signal; s4, carrying out beam forming on the vibration signal, and outputting a final beam forming result through combined judgment with beam forming of the underwater acoustic signal; obtaining the time difference between the P wave and the S wave based on the speed difference between the P wave and the S wave; and S5, finally realizing target identification based on feature extraction of the underwater acoustic signal, realizing target orientation estimation based on a combined comparison result in the step S4, and determining the distance between the target and the optical cable array based on the time difference between the P wave and the S wave. The method effectively solves the technical problems that in the prior art, accurate positioning is difficult due to single dependence on underwater acoustic signals, and target recognition cannot be achieved only according to vibration signals.

The present invention discloses a method, apparatus, terminal device and storage medium for processing ECG data for heart failure assessment, which obtains ECG data of a patient; performs peak screening processing on an ECG in the ECG data to obtain a plurality of ECG waves; determines R waves, P waves, Q waves, S waves and T waves in the plurality of ECG waves; extracts ECG acquisition parameters and ECG waveform parameters corresponding to each ECG wave in the ECG data based on generation time; establishes a coordinate system for each ECG wave with the ECG waveform parameter type of the ECG wave as the horizontal coordinate and the ECG waveform parameter value as the vertical coordinate, determines the ECG waveform type of the ECG wave based on the parameter value range of historical ECG waveform types, and summarizes the ECG waveform type and the ECG acquisition parameters to generate ECG analysis data of the patient; performs heart failure assessment based on the obtained accurate ECG analysis data in combination with the user's body impedance, which is beneficial to improving the accuracy of heart failure assessment.

The invention provides a pile foundation seismic response analysis method and system considering a transverse inertial effect, a storage medium and application of the pile foundation seismic response analysis method and system. Aiming at the problem that a pile-soil coupling vibration model in the prior art neglects the fractional order viscoelasticity characteristic of unsaturated soil and the transverse inertia effect of an uncoupled pile body, the invention provides a large-diameter end-bearing pile seismic response analysis model based on a fractional order viscoelasticity constitutive structure. A fractional derivative is introduced to describe the stress relaxation characteristic of an unsaturated soil skeleton, a Rayleigh-Love beam theory is combined to characterize the transverse deformation effect of a pile body, and a pile-soil coupling vibration control equation under the action of earthquake P waves is constructed. Laplacian transformation, potential function decomposition and operator decomposition methods are adopted to realize equation dimensionality reduction solution, a pile-soil interface continuous condition and a boundary constraint condition are synthesized, a theoretical analytical solution of vibration response in a frequency domain is obtained, and a theoretical basis is provided for aseismic design of a pile foundation in an unsaturated land area.