53 results about "Akaike information criterion" patented technology

The invention discloses a coal and rock mass fracture seismic source and stress field joint inversion method and system, and the method comprises the following steps: arranging a multi-channel sonic sensor in a coal and rock mass, and collecting original acoustic emission / microseismic signal data; the head wave arrival time / amplitude of the filtered effective waveform data is automatically picked up through an improved AIC (Akaike Information Code), quantitative evaluation is carried out on the quality of a picked-up result, and high-precision positioning is carried out on waveform signals with a high signal-to-noise ratio and a low signal-to-noise ratio; constructing a tension fracture seismic source model based on a displacement discontinuous tensor, obtaining a moment tensor component of a fracture event, and calculating to obtain a spatial orientation, a tension-shear attribute, a deviation angle alpha, energy release and a moment magnitude of the fracture; establishing a stress field inversion model of the tension-shear fracture seismic source, and realizing the inversion of the stress direction and size of the tension-shear fracture seismic source; and a slip direction included angle minimization objective function is established for iteration until convergence conditions are met, and joint iteration inversion of the fracture and the stress field is realized.

The invention discloses a micro-seismic signal arrival time pickup method based on fuzzy clustering and variational mode decomposition. The method comprises the following steps: firstly, calculating characteristic functions of attribute characteristics such as a micro-seismic signal mean value, power and kurtosis and normalizing the characteristic functions to obtain a characteristic matrix, then primarily picking up a micro-seismic P-wave initial movement position and a first arrival moment through a fuzzy clustering algorithm, and extracting an effective time window by taking the first arrival moment as a reference point; a variational mode decomposition algorithm is utilized to decompose signals in a time window into K intrinsic mode function components, AIC function values of the components are calculated by means of an akaike information criterion algorithm, a minimum value point is picked up to serve as first arrival time, energy ratios of all the components are calculated, and final micro-seismic P-wave first arrival time is obtained through weighted calculation. The method effectively deals with the low signal-to-noise ratio environment of the underground coal mine, has higher pickup precision and reliability compared with a traditional pickup method, can provide accurate micro-seismic occurrence time and position information for mine safety early warning, and powerfully guarantees the safety production of the coal mine.

The invention provides a spatial variability soft rock tunnel deformation prediction method based on deep learning, and belongs to the technical field of tunnel deformation prediction, and the method comprises the steps: constructing a surrounding rock parameter database, and determining the optimal edge distribution of different surrounding rock parameters through employing an akaike information criterion; constructing a reference value initial model, and generating a related standard uniformly distributed random field; performing equal probability transformation on the correlation standard uniformly distributed random field to generate a cross-correlation non-Gaussian distributed random field; encoding the data matrix of the non-Gaussian distribution random field into a parameter field image; and performing soft rock tunnel deformation prediction on the parameter field image by adopting a trained deformation prediction model. According to the method, parameter independent hypothesis limitation is broken through, the multi-parameter coupling space variability is accurately represented, and the prediction precision is remarkably improved.

The invention relates to the technical field of GNSS deception jamming resistance, and particularly discloses a GNSS induced deception jamming resistance method and system, and the method comprises the steps: obtaining a GNSS signal, and calculating a code phase difference and a carrier phase difference; traversing a modal component number K value, performing variational modal decomposition on the code phase difference according to each K value to obtain a modal component, calculating a criterion value of a corresponding Akaike information criterion, and selecting K which enables the criterion value to be minimum as an optimal modal component number; performing variational mode decomposition on the code phase difference and the carrier phase difference based on the optimal mode component number; modal components containing the induced cheating information are screened out, and a screening result is obtained; and carrying out positioning calculation on the GNSS signal based on a screening result. The method is obviously superior to the prior art in the aspects of dynamic adaptability, feature separation precision and multi-information fusion, and a more efficient and robust solution is provided for GNSS induced deception jamming resistance.

The invention discloses a cascade water-wind-light multi-power system modeling method based on time-space correlation, and the method comprises the steps: building an autoregression integral moving average-generalized autoregression condition heterovariance model of wind power photovoltaic power and cascade hydropower station runoff, building an R-vine, C-vine and D-vine Copula function spatial correlation model through employing a generalized autoregression condition heterovariance value, and carrying out the modeling of a cascade water-wind-light multi-power system. And selecting an optimal spatial correlation model according to a Bayesian information criterion, an akaike information criterion and a logarithm likelihood value goodness of fit index. Based on the optimal spatial correlation model, constructing a scene generation method of Latin hypercube-Copula function sampling, performing scene generation on the wind power photovoltaic power and the runoff volume of the cascade hydropower station, and constructing a scene reduction method of contour coefficient optimization Kmeans clustering to reduce the generated scene. The invention provides a time-space correlation cascade water-wind-light multi-power system modeling method, and aims to analyze water-wind-light multi-energy complementary operation characteristics and provide reference for a scheduling mode among water-wind-light multi-energy.

The invention belongs to the technical field of ecological restoration, and provides a mountainous river-oriented fish diversified habitat natural shaping restoration method, which comprises the following steps of: determining fishes with a relative important index greater than 5% in a river as target fishes; the habitat factors of the river are screened, the habitat factors with colinearity are divided into a factor group, and only one habitat factor is reserved in each factor group to serve as a representative habitat factor; carrying out modeling on the resource distribution and representative habitat factors of the target fishes by adopting a generalized additive model; screening to obtain a generalized additive model with a minimum akaike information criterion, and determining an optimal habitat factor combination; and determining a restoration scheme according to the habitat factor combination, and adjusting the habitat of the target fish. According to the restoration method, the natural process is simulated to remodel the diversity of the mountain river habitat, the key problems of hydrological connectivity fracture, habitat simplification, engineering and ecological contradiction and the like are solved, and fish population restoration and ecological system self-maintenance are achieved.

The present invention relates to the technical field of GNSS anti-spoofing interference, and specifically discloses a GNSS anti-induced deception interference method and system, the method comprising: acquiring a GNSS signal and calculating a code phase difference and a carrier phase difference; traversing the modal component number K value, performing variational modal decomposition on the code phase difference according to each K value to obtain modal components, calculating the corresponding criterion value of the Akaike information criterion, and selecting K that minimizes the criterion value as the optimal number of modal components; performing variational modal decomposition on the code phase difference and the carrier phase difference based on the optimal number of modal components; screening out modal components containing induced deception information to obtain screening results; and performing positioning and solving on the GNSS signal based on the screening results. The present invention is significantly superior to existing technologies in terms of dynamic adaptability, feature separation accuracy, and multi-information fusion, and provides a more efficient and robust solution for GNSS anti-induced deception interference.

The application discloses a site micro-vibration vibration source identification method based on machine learning, signals measured by acceleration sensors at the edge of a measured site and acceleration sensors at the center of the measured site are processed in sequence, including denoising, transient impact signal extraction, feature extraction, clustering processing of dynamic time warping distance of the transient impact signal based on a feature matrix, GMM modeling, obtaining a model, verifying the number of independent Gaussian components of the model, selecting the number of components with the minimum value as the optimal model according to the Akaike information criterion value, and performing maximum a posteriori probability identification on the transient impact signals captured by the sensors; GMM parameter comparison is performed on sets that may belong to the same category, two sets with highly overlapped models are merged into the same set, and the set is used as vibration data sets generated around the site and having an impact on the center position of the site.

The invention discloses a distribution network line fault trip prediction method and device, a medium and a product, and relates to the field of power distribution network operation and maintenance, and the method comprises the steps: obtaining fault data of a distribution network line; constructing a fault data set; according to the fault data set, constructing a trip prediction model based on an akaike information criterion; an exogenous variable corresponding to the tripping prediction model is determined; constructing a limit gradient lifting model based on a grid search method according to the exogenous variables and the historical fault trip times; constructing a fusion model according to the trip prediction model and the limit gradient lifting model; respectively determining prediction errors of the tripping prediction model and the fusion model, and determining a final prediction model according to the prediction errors; determining the number of fault tripping times of the month to be predicted by using the final prediction model; and regularly obtaining fault data of the distribution network line to regularly update the final prediction model. The monthly fault trip times of the distribution network line can be efficiently predicted, and the operation and maintenance reliability of the distribution network is improved.

The invention discloses an unbalanced power distribution network risk analysis method based on a multistage sparse polynomial, and belongs to the technical field of power distribution network risk analysis. The method comprises the steps of establishing a power distribution network input-output model; approximating the input and output model of the power distribution network as a linear regression model; calculating an ANCOVA index and a P value of the linear regression model, and sparsifying the input and output model of the power distribution network to obtain a preliminary sparsifying input and output model of the power distribution network; restraining a high-order cross term of the preliminary rarefaction power distribution network input and output model according to a preset truncation rule, shrinking a coefficient of the preliminary rarefaction power distribution network input and output model according to the optimized polynomial chaos expansion coefficient, and taking a preset residual error akaike information criterion as a model balance criterion to obtain a final rarefaction power distribution network input and output model; and obtaining a power distribution network risk analysis result according to the final sparse power distribution network input and output model. The method can provide a quantitative basis for operation and control of the power distribution network.

The invention discloses a partial discharge positioning detection method and system based on acoustoelectric joint detection and multi-dimensional feature fusion, and belongs to the technical field of intelligent power grid state detection. According to the method, an ultra-high frequency (UHF) signal is used as a discharge absolute time reference, an akaike information criterion (AIC) is adopted to adaptively pick up an ultrasonic weak wave head, and the problem of failure of a traditional threshold value method under a low signal-to-noise ratio is solved; meanwhile, a four-dimensional state vector containing space coordinates and equivalent sound velocity is constructed, a weighted Gaussian-Newton iterative algorithm is used for solving, a temperature field sound velocity correction model and a grid search initial value estimation algorithm are introduced, and the problems of iterative divergence and local optimal solution are solved. In addition, the method integrates algorithm modules of intelligent noise reduction, physical field adaptive correction, sound intensity assisted positioning and the like, and realizes high-precision and high-robustness positioning of the partial discharge source.

The invention discloses an ionosphere modeling method based on dynamic adaptive function optimization, and belongs to the technical field of earth space environment and satellite navigation, and the method comprises the steps: adaptively obtaining NeQuick model ionosphere background data of a target region in a preset time period through physical driving and a dynamic boundary; according to the NeQuick model ionized layer background data, determining an optimal polynomial order combination from a preset candidate order set by using an akaike information criterion; and on the basis of observation data of a global navigation satellite system, constructing a regional ionosphere refinement model by adopting the optimal polynomial order combination. According to the method, the simplest model is selected on the premise of ensuring the precision, and the calculation efficiency is improved.

This invention relates to a method for highlighting spike discharge data of intracranial electroencephalogram (EEG) signals. The method involves collecting background noise data from the patient's brain without neuronal discharges, preprocessing the background noise data, automatically selecting the optimal order of an autoregressive (AR) model using the Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC), estimating the AR model coefficients using the Yul-Walker equation, and constructing a background noise model. The intracranial EEG signals to be processed are then subjected to high-pass filtering. A short-time Fourier transform (STFT) and window-based frame-by-frame processing strategy are used to subtract the spectrum of the signal from the noise. By adjusting the parameters of the spectral subtraction, noise removal and preservation of neuronal signal features are achieved.

The invention relates to the technical field of global navigation satellite system data processing, and discloses a GARCH-based GNSS cycle slip detection threshold real-time determination method, which comprises the following steps: acquiring GNSS double-frequency observation data, constructing a Morse-Wiberner and geometric distance-free combination observation value, and carrying out differential processing to obtain a noise time sequence. Secondly, traversal calculation is carried out based on an akaike information criterion, and the optimal order and the sliding window length of the GARCH model are determined; then, a current epoch condition standard deviation is predicted using the model, and a multiple factor is adjusted in combination with an environmental state to construct a dynamic detection threshold. And finally, comparing an absolute value of the noise sequence with a threshold value to judge a cycle slip, and calculating an integer solution by constructing an observation equation to complete repair. According to the method, the heterovariance characteristic of the observation noise is modeled by using the GARCH model, so that the adaptive adjustment of the detection threshold along with the environmental noise level is realized, the misjudgment rate in a complex environment is effectively reduced, and the detection accuracy is improved.

This invention discloses a multi-information source mechanism calculation method based on earthquake early warning station network observations. The method includes selecting velocity and acceleration recorded waveforms to form an input database; obtaining the P-wave initial motion polarity through manual picking and template matching techniques; calculating the source mechanism using a grid scanning method, and using solutions with an inconsistency ratio below a threshold as preliminary solutions; manually screening P-wave and S-wave phase waveforms; and jointly using initial motion polarity, S / P amplitude ratio, and body wave waveform as constraints to calculate the probability value of each preliminary solution. The probability distribution of clustered solutions is obtained through the Akaike information criterion and cluster analysis to obtain the optimal source mechanism calculation result. This invention, based on joint inversion calculations using acceleration records from early warning stations and velocity records from seismic stations, lowers the lower limit of the magnitude of earthquakes with computable source mechanisms. The comprehensive use of three types of information constraints improves the reliability of the results, and the calculation of the probability distribution of clustered solutions reduces the adverse effects of multiple solutions, thus improving the stability of the results.

The application discloses a water engineering scheduling runoff intelligent simulation method, comprising: edge distribution function fitting: using the matrix method to calculate the mean, dispersion coefficient and skewness coefficient, and adopting a P-III type curve to construct a daily runoff frequency curve; wisdom knowledge set construction: using a time series model to establish wisdom knowledge sets of each parameter in a mixed Copula set, using a cooperative search algorithm, taking the minimum of the Akaike information criterion of the theoretical distribution and the empirical distribution as the target, and calibrating the wisdom knowledge set parameters; daily runoff simulation: based on the Bayesian rule, a mixed conditional Copula set is established, and Latin hypercube sampling is used to generate a long sequence of daily runoff. The method reduces the difficulty of daily runoff simulation by the mixed Copula, reduces the number of parameters, improves the fitting precision of the mixed Copula, speeds up the construction speed of the mixed Copula set, and provides a valuable technical tool for basin management.

The present application belongs to the field of blast furnace system health management and prediction, and specifically discloses a blast furnace degradation model calibration and residual life prediction method under a noise environment. The method is based on a fractional Brownian motion with historical dependence, constructs an initial degradation model of the blast furnace system, estimates the model parameters by the maximum likelihood estimation method, and determines the optimal model structure of the blast furnace system temperature change by the Akaike information criterion. In addition, the present application also designs a model calibration trigger mechanism; when the mechanism is not triggered, the Bayesian fusion particle filter is used to estimate the potential degradation state of the system and update the model parameters; when the mechanism is triggered, a prediction error model is constructed and the degradation model structure is calibrated to adapt to complex working conditions. Finally, the probability distribution function of the residual life of the blast furnace system is derived, and the residual life prediction of the blast furnace system is carried out according to the calibrated degradation model.

The invention discloses an eVTOL collaborative design optimization method based on improved genetic programming, and the method comprises the steps: firstly constructing pneumatic, propulsion and dynamics high-fidelity models of an eVTOL takeoff stage, building an original collaborative optimization problem, and solving the problem through a direct transcription method; secondly, multi-target genetic programming fusing a root mean square error, an akaike information criterion, a predicted residual sum of squares and a Bayesian information criterion is creatively adopted, and an agent model for calculating thrust is constructed to replace a high-fidelity thrust model to solve and obtain optimal design parameters. The error between the optimal solution obtained by the method and the optimal solution obtained by using the high-fidelity model is very small, the calculation cost is remarkably reduced, and multidisciplinary efficient collaborative optimization is realized.

The invention relates to the technical field of paving thickness measurement, in particular to an ultrasonic paving thickness detection method and system, and the method comprises the steps: employing an Akaike information criterion for an echo measurement sequence obtained through collection of a paving road, and obtaining an AIC value of all elements of the echo measurement sequence based on the numerical distribution of AIC values of all elements of the echo measurement sequence; an echo positioning sequence and an echo AIC statistical sequence are obtained; element distribution and change characteristics in the echo positioning sequence are analyzed, and the change trend of each element in the echo AIC statistical sequence is obtained to obtain a signal mutation trend degree; obtaining an echo development trend degree based on a straight line fitting result of each element in the echo AIC statistical sequence; and confirming a primary echo moment, and calculating the paving thickness by combining the propagation speed of the ultrasonic waves in the paving mixture. The invention aims to improve the accuracy of paving thickness measurement.

The application discloses a fault probability monitoring and regulation system and method, relates to the technical field of fault monitoring, and collects to-be-determined related parameters in each period through a collection module, determines the optimal order of an autoregressive model of each related parameter by using the Akaike information criterion, and verifies the to-be-determined related parameters in each period, decides whether to perform secondary collection and replacement based on a verification result, generates a related parameter vector of each period, and starts a window translation along a period in a reverse order from the related parameter vector of each period through a monitoring module, calculates and arranges a comprehensive feature vector in each window, extracts and generates a core feature vector set through a time series self-encoder, trains a hidden Markov model by using a forward-backward algorithm, and obtains the fault probability of each period based on Bayesian inference, compares the fault probability with a probability threshold value to determine whether a fault exists, and executes regulation when the fault exists, so that high-precision fault monitoring and regulation based on collection self-checking, multi-dimensional data fusion and probability inference are realized.

The invention discloses a low-signal-to-noise-ratio impact signal first arrival moment extraction method, which comprises the following steps of: firstly, calculating an AIC (Akaike Information Criterion) curve of an original low-signal-to-noise-ratio impact signal according to an AIC (Akaike Information Criterion); secondly, smoothing preprocessing is conducted on the AIC curve, and then a first-order difference curve is calculated to highlight abrupt change characteristics at the first arrival moment of the signals; and finally, detecting the starting point of the first significant mutation signal through statistical analysis of the mutation characteristics of the first-order difference curve, and determining the moment corresponding to the starting point as the first arrival moment of the impact signal. According to the method, first arrival time extraction is directly carried out on the original low-signal-to-noise-ratio impact signal based on the statistical characteristics of the AIC curve of the original signal, denoising preprocessing does not need to be carried out, and the method has the technical advantages of being high in extraction precision, high in anti-interference capacity, low in calculation complexity and wide in application range; and the method shows good stability and robustness under different signal-to-noise ratio conditions.

Disclosed are methods and techniques of detecting network anomalies and responding to the anomalies once detected. The methods, for example, include receiving, by a model executed by a processor, real-time log data of an operating network; parsing, by the model executed by the processor, the log data to identify one or more metrics; determining, by the model executed by the processor, a seasonality of the one or more metrics; determining whether the model should use an autoregressive model if seasonality is detected; and on determining that an autoregressive model should be used, training a model based on determining a grid search for a parameter based on an Akaike information criterion.

The present invention relates to a method for estimating the kinetic parameters of hepatocellular carcinoma PET / CT imaging based on an improved MCMC algorithm and its application, belonging to the technical field of liver cancer auxiliary diagnosis. The present invention first obtains prior information and sample information. Secondly, based on Bayesian theory, the posterior distribution of the kinetic parameters of the hepatocellular carcinoma PET / CT imaging tracer is derived, and the Markov Chain Monte Carlo (MCMC) method is used to sample the parameter posterior distribution. The sampling algorithm adopted is an improved MCMC algorithm. The global optimization and local optimization of the flower pollination algorithm and the adaptive substitution probability are introduced into the standard MH algorithm to improve the convergence efficiency and accuracy of the MH algorithm, and estimate the kinetic parameters. Finally, the differences in the tracer kinetic parameters between the hepatocellular carcinoma tissue images and the surrounding healthy liver tissue images are analyzed and compared, and the Akaike Information Criterion (AIC) and Root Mean Square Error (RMSE) are used to compare the performance of the algorithm before and after improvement.

The present invention provides a method and device for detecting structural damage of a nonlinear system, which relates to the field of nondestructive testing technology. The method includes: applying an excitation force to a test piece to be inspected, vibrating the test piece to be inspected based on the excitation force, collecting signals from the test piece to be inspected to obtain an excitation signal and a response signal, combining the excitation signal and the response signal with multiple imported orders to construct initial NARX models of different orders, screening a target NARX model from the multiple initial NARX models based on the Akaike information criterion and optimizing it, comparing the order information in the optimized NARX model to be inspected with the order information in the reference NARX model, and determining whether the test piece to be inspected has defects based on the order comparison results. A nonlinear model is constructed and optimized based on the signal data of the test piece to be inspected, so that the model accurately represents the structural characteristics of the test piece to be inspected, thereby determining whether the test piece to be inspected has defects.

The invention relates to the technical field of intelligent power grid and distributed energy management, in particular to a distributed energy collaborative management system based on edge intelligence, which comprises an end side sensing layer, an edge intelligent node, a cloud collaborative layer, a secure communication layer and a visual interaction layer. According to the scheme, the multi-scale weighted regression model is constructed, the multi-scale space-time characteristics of the distributed energy sources are analyzed, the self-adaptive weight matrix is constructed to accurately capture the regional characteristics, the scale parameters are optimized by using the corrected akaike information criterion, and wind and light curtailment and power supply shortage are remarkably reduced; according to the scheme, the optimized self-adaptive differential evolution algorithm is adopted for regional global optimization scheduling, the scaling factor and the crossover probability are dynamically adjusted based on the historical memory bank, millisecond-level response is achieved in cooperation with a linear population reduction mechanism, the system operation cost is reduced, and the energy utilization efficiency is improved.

The application discloses a signal first arrival time picking method and system, electronic equipment and a storage medium. The method comprises the following steps: obtaining a target signal by collecting and preprocessing to-be-processed data; performing amplitude-aware permutation entropy calculation on the target signal to determine a target interval containing a first arrival time; picking an improved Akaike information criterion minimum value point in the target interval as the first arrival time; and completing a geophysical data processing process according to the first arrival time. The embodiment of the application can improve the picking efficiency and the signal picking accuracy under the condition of low signal-to-noise ratio, and can be widely applied to the near-surface engineering geophysical technology field.

The invention belongs to the technical field of deflection calculation, and relates to an elastic winding drum deflection inversion calculation method and device based on strain analysis. The method comprises the following steps: performing strain analysis on an elastic drum to obtain elastic drum fixed end strain data and elastic drum free end deflection data; generating a strain and deflection data sample set of the elastic winding drum according to the strain data of the fixed end of the elastic winding drum and the deflection data of the free end of the elastic winding drum; according to the strain and deflection data sample set, obtaining a polynomial relationship between the strain and the deflection, and writing the polynomial relationship into a matrix form to obtain a matrix relationship between the strain and the deflection; calculating coefficient vectors of the matrix relationship under different orders, and determining the order of a polynomial relationship in combination with an akaike information criterion to obtain a polynomial equation of strain and deflection; and according to the polynomial equation and the obtained strain data, carrying out inversion calculation to obtain the free end deflection of the elastic winding drum. According to the invention, deflection inversion calculation can be carried out.

The present invention discloses a method and system for joint inversion of coal rock fracture source and stress field. The inversion method comprises the following steps: deploying a multi-channel acoustic wave sensor inside the coal rock mass to collect original acoustic emission / microseismic signal data; automatically picking up the first wave arrival time / amplitude of the filtered effective waveform data by improving the Akaike Information Criterion (AIC), quantitatively evaluating the quality of the picking results, and performing high-precision positioning of waveform signals with high signal-to-noise ratio and low signal-to-noise ratio; constructing a tensor-based fracture source model based on a displacement discontinuity tensor to obtain the moment tensor component of the fracture event, and calculating and obtaining the spatial orientation, tension-shear properties, and offset angle of the fracture. α , energy release and moment magnitude; establish a stress field inversion model of the tensile shear rupture source to realize the inversion of the stress direction and magnitude of the tensile shear rupture source; establish an objective function for minimizing the slip direction angle and iterate until the convergence condition is met to realize the joint iterative inversion of the crack and stress field.

The invention relates to the technical field of signal processing, in particular to a pseudo-period detection method and device for a coordinate time series, electronic equipment and a medium, and the method comprises the steps: carrying out the data preprocessing of the coordinate time series of a target navigation satellite system, carrying out the periodic signal detection, and obtaining a plurality of target periodic signals meeting a preset detection condition, calculating a reference Akaike information criterion of the preprocessed coordinate time sequence under zero hypothesis, and performing fitting based on the single-cycle model of each target cycle signal and the preprocessed coordinate time sequence to obtain a first new coordinate time sequence, and obtaining a real periodic signal of the coordinate time sequence based on a comparison result between the target Akaike information criterion corresponding to each target periodic signal and the reference Akaike information criterion. Therefore, the problems of pseudo-period interference and the like caused by gross error residue, incomplete elimination of step signals and poor precision of a period detection method in the prior art are solved.

Embodiments of the present application relate to the field of microseismic monitoring, and disclose a microseismic first arrival picking method, device, equipment and storage medium. The center time point of the time window of the microseismic signal is calculated according to the propagation path of the microseismic signal and the preset layer velocity, and the length of the time window is determined according to the center time point and the preset number of time sampling points, so as to constrain the time window of the first arrival point of the microseismic signal picked by the Akaike information criterion method, and the layer velocity is constantly updated in an iterative manner, and the time window is further updated, so as to determine the best first arrival point of the microseismic signal. Compared with the traditional Akaike information criterion method, the application has better noise immunity, can avoid the "multi-peak" phenomenon of the Akaike information criterion method under the condition of low signal-to-noise ratio, and thus avoids the misjudgment of the first arrival point.