28 results about "Energy minimization" patented technology

The invention relates to the technical field of nanofluid enhanced heat transfer simulation, in particular to a molecular dynamics-based mixed nanofluid heat transfer characteristic analysis method, which comprises the following steps of: constructing a molecular model comprising copper, graphene nanoparticles and an ionic liquid / water mixed base solution, and calculating the heat transfer characteristic of the mixed nanofluid through a distribution potential function and a periodic boundary condition; and performing energy minimization and molecular dynamic relaxation to achieve thermodynamic equilibrium. A steady-state temperature gradient is established based on the unbalanced molecular dynamics principle, and the thermal conductivity is calculated according to the Fourier law. Microcosmic parameters such as a radial distribution function, mean square displacement and the number of hydrogen bonds are synchronously analyzed, and a coupling mechanism of interaction and heat transfer performance of the nano-particle-base fluid interface is disclosed. According to the method, the heat transfer characteristic of the mixed nanofluid can be clarified on the molecular scale, and an efficient and reliable analysis means is provided for nanofluid screening and design in application of a solar heat collector, a heat management system and the like.

The invention discloses a method for analyzing the dynamic evolution process of a methane leakage nano-channel in cold spring hydrate, and belongs to the field of natural gas hydrate science and engineering. The method comprises the following steps: constructing a hydrate molecular model and selecting a methane and water molecule force field; a hydrate block and a middle gas channel are arranged in the simulation box, and the gas channel is filled with a methane and water mixture to construct an initial simulation system; performing molecular dynamics simulation on the system under the conditions of energy minimization, pre-balance and different pressures and methane accelerations, and remodeling a gas channel dynamic evolution process in a hydrate growth process; carrying out analysis and three-dimensional rendering on the simulated track, and drawing a three-dimensional morphology graph of the gas channel changing along with time; quantitatively extracting surface area data of the gas channel, and calculating key kinetic parameters, including attenuation rate and half-life period, of the evolution process of the gas channel. According to the invention, the quantitative characterization of the dynamic evolution process of the methane leakage nano-channel in the cold spring environment is systematically realized on the molecular scale for the first time, and the contraction closing behavior of the gas channel under the condition of no methane leakage and the dynamic maintaining or adjusting process of the gas channel under the condition of continuous seepage can be accurately simulated and distinguished; the limitation of a traditional static pore model is overcome, and key basic parameters and method support are provided for revealing the interaction mechanism of methane migration and hydrate growth and constructing a cross-scale seepage model.

The invention provides a micro-scale steady-state meniscus model construction method and device, equipment and a medium, and belongs to the field of micro-scale in-pore flow simulation calculation. According to the method, a micro-scale pore meniscus molecular model containing an upper flat plate, a lower flat plate and a left fixed flat plate is constructed, and liquid is filled in the micro-scale pore meniscus molecular model; the method comprises the following steps: performing energy minimization to obtain a liquid static structure, moving a flat plate to obtain a steady-state structure, performing stress analysis on the steady-state meniscus structure, and counting horizontal tangential force borne by the flat plate under the unit length of a liquid-solid interface within preset time. By constructing the micro-scale pore meniscus molecular model, fluid dynamic behaviors of a three-phase contact line area under the micro-scale can be simulated, data which is difficult to obtain in experiments such as a three-phase interface liquid molecule migration mechanism and the like can be obtained, and the problem that a current experimental method is difficult to reveal a micro-nano liquid-solid interface flow mechanism is solved.

The application relates to a calculation method of a high-temperature alloy impurity diffusion coefficient and an atomic mobility parameter based on a molecular dynamics simulation, which comprises the following steps: parameter setting; model establishment; setting a potential function as a multicomponent EAM alloy potential; energy minimization and setting an initial atomic speed; calculating RDF values of a model system and drawing RDF curves; calculating MSD values of all atoms and respectively outputting MSD values of various atomic groups according to various different metal atom types; defining thermodynamic physical quantity output and commanding each step to output and save; giving a random initial speed to the system, and running under each system in three stages: a temperature rising stage, a temperature keeping stage and continuous temperature keeping and relaxation, and outputting MSD values; calculating the impurity diffusion coefficient of the system based on the MDS values; and taking the measured impurity diffusion coefficient as an input value, and obtaining the atomic mobility parameter through DICTRA software. Compared with the prior art, the application has the advantages of simple method, accurate numerical value and the like.

The invention relates to a three-dimensional modeling method for a reef limestone microstructure, which comprises the following steps of: scanning to obtain an original grayscale image of a rock core sample, and preprocessing the original grayscale image; based on the preprocessed grayscale image, extracting a sample shell mask; establishing a Gaussian mixture model of sample gray level distribution, and determining a segmentation threshold value; classifying and segmenting the volume data through the segmentation threshold to obtain a pore candidate region, a matrix candidate region and a pending region; performing morphological closed operation on the candidate region to obtain a pixel set belonging to the pore and a pixel set belonging to the matrix; and classifying and segmenting the pending region by adopting an energy minimization algorithm based on graph segmentation. According to the method, the GMM dual-threshold segmentation and shell mask construction technology is adopted to accurately extract the core area, and air at the core boundary is prevented from being wrongly segmented into pores, so that the real boundary of the core is completely reserved; according to the method, the pores can be finely judged, the automatic threshold segmentation result is stable and reliable, and manual intervention is greatly reduced.

This invention discloses a method for analyzing the wear of micro-milling tools, relating to the fields of wear mechanism research and molecular dynamics analysis in micro-milling. Specifically, it includes: constructing atomic models of the tool and workpiece, incorporating the geometric features of the cutting edge arc; setting up a simulation environment and dividing the workpiece atoms into boundary layers, isothermal layers, and Newtonian layers; assigning potential functions to the atoms and performing energy minimization and relaxation; constructing a dimensionless scale ratio loading condition based on the equivalent uncut thickness and the cutting edge arc radius; extracting the mechanothermal process quantities and quantitatively characterizing the wear state according to the number of detached atoms, the change in effective cutting edge radius, and interface mixing indices. The aim is to achieve atomic-level quantitative characterization of the wear process of micro-milling tools.

The application provides a method for predicting the ice inhibition performance of nanometer silicon dioxide based on molecular dynamics simulation. The application comprises the following steps: establishing a functionalized nanometer silicon dioxide model, and establishing an ice-containing water box model according to the size of the functionalized nanometer silicon dioxide model; assembling the nanometer silicon dioxide model and the water box model; selecting a suitable force field for the constructed model and performing energy minimization optimization. Molecular dynamics simulation calculation is performed on the silicon dioxide-ice-water model after energy minimization, and analysis is performed on the model after calculation. The root mean square deviation, root mean square fluctuation and binding energy are analyzed, and the number of hydrogen bonds is counted. Based on the analysis data, the feasibility and effect of the material for inhibiting ice are judged. The application constructs different functionalized materials and evaluates their behavior in the ice-water system, thereby greatly reducing the application cost, saving experimental resources, and guiding the development of a freezing protectant based on nanometer materials.

This invention discloses a method for obtaining interface thickness parameters based on a microscopic all-atom composite model, belonging to the field of molecular simulation. First, using MS software, a periodic model of polymer A and polymer B with local energy minimization is constructed. After combining the two, geometric optimization and kinetic equilibrium optimization are performed to obtain a globally energy-minimized all-atom composite model. Finally, using secondary programming techniques in the Perl script module of MS software, the atomic concentration distribution data points of the all-atom composite model are calculated, the atomic concentration distribution curve is plotted, and smoothing is performed. A criterion for determining the presence of interface thickness in the all-atom composite model is set to determine whether interface thickness exists. A formula for calculating the interface thickness of the all-atom composite model is defined to calculate the interface thickness of the all-atom composite model with interface thickness, thereby quickly, efficiently, accurately, and quantitatively calculating the interface layer thickness of the polymer composite material. Furthermore, this invention also has the advantages of convenient operation.

The application discloses a kind of molecular dynamics simulation methods of gas leakage in metal micro crack, the method is first in supercell iron model pre-set hollow microchannel, obtain the iron model containing micro crack;Energy minimization processing and molecular dynamics simulation are carried out based on periodic gas model, obtain the gas model after relaxation;Vacuum layer, iron model containing micro crack, relaxation gas model, isolation layer splicing are established by micro crack gas leakage interface model, then molecular dynamics simulation is carried out to interface model, according to the leakage rate of gas in micro crack is calculated according to simulation trajectory result.This application method has the advantages of high calculation efficiency, low simulation cost, can construct the complete micro crack gas leakage interface model of nanometer scale, from microcosmic angle of view analyzes the leakage rule and mechanism of gas molecule in micro crack, obtains the internal relation of micro crack size, structure and gas leakage evolution, lays the foundation for the sealing performance, life and reliability of precision instrument promotion.

The application discloses a kind of calculation methods for determining C-S-H pore structure in irreversible shrinkage range, comprising the following steps: establishing C-S-H pore structure model;Using Helmholtz free energy to describe the energy of C-S-H pore structure model;Establish the relationship between the energy part related to physical entropy and pore structure;Establish the relationship between the energy part related to total volume and pore structure;Establish the relationship between the energy part related to surface energy and pore structure;Determine the calculation parameter, and carry out energy calculation to C-S-H pore structure model in given calculation range;Determine the effective C-S-H pore structure energy calculation range;Based on the principle of energy minimization, determine the structure information and structure parameters corresponding to stable C-S-H pore structure model.The application solves the problem that various experimental measurement methods and structure models in the prior art are qualitatively described in the structure characterization of C-S-H pore structure under different humidity, without quantitatively describing C-S-H pore structure information.

This invention belongs to the field of seismic data processing and relates to a method and system for consistent processing and adaptive fusion of seismic attributes in multiple contiguous seismic work areas. The method includes: S1: dividing the spatial coordinates of the original seismic attribute data into a grid and establishing a gridded data index; then, based on the gridded data index, executing a localized search strategy to extract the neighboring data point set of the seismic data points to be processed; S2: eliminating systematic value range differences in the original seismic attribute data based on the neighboring data point set and stitching them together to obtain consistent fused seismic attribute data; S3: using a hierarchical interpolation strategy to fill in blank areas to obtain continuous seismic attribute interpolation data; and S4: performing global optimization based on Gaussian pyramid decomposition and an energy minimization model to obtain the final fused seismic attribute data. This method achieves better work area stitching and interpolation, efficiently avoiding the attribute non-fusion effect of traditional work area stitching methods, thereby realizing the goal of consistent seismic attribute stitching across multiple work areas.

The invention discloses a crack global optimization reconstruction method and system for structure detection. The method comprises the following steps: preprocessing a structure surface image, detecting cracks, and analyzing a connected region to obtain a discrete crack fragment set; extracting geometrical characteristics of each fragment, and calculating a fragment pair space adjacent distance to construct a candidate connection set and a crack connection graph; constructing a feature vector based on the connection probability features, and combining a likelihood model and a structure prior model to form a probability graph model; and after the problem is converted into an energy minimization problem, performing iterative optimization through a maximum posterior probability inference algorithm dynamically updated by a structural state, screening an optimal connection relation set, and combining crack segments to realize global optimization reconstruction. According to the invention, the accuracy and continuity of crack detection in a complex scene can be improved.

The invention relates to the technical field of heat insulation processes, and discloses a molecular dynamics optimization process for regulating and controlling a microporous structure of an aerogel layer, which comprises the following steps of: constructing an initial atomic model of aerogel molecules and generating a corresponding force field parameter file on the basis of a multi-software cooperative workflow; based on the three-dimensional space size of the initial atomic model, a space index data structure is initialized, and the data structure is used for efficiently storing and querying atomic space position information; and performing energy minimization calculation on the initial atomic model to obtain an initial equilibrium system. According to the method, an automatic modeling process is constructed by integrating multiple software, the calculation efficiency is greatly improved by adopting a spatial index and dynamic updating mechanism based on a hash table, super-large-scale system simulation is supported, and meanwhile, an intelligent optimization engine is introduced to automatically optimize simulation parameters; the problems of flow splitting, limited calculation scale and dependence on manual trial and error parameter adjustment in the prior art are effectively solved.

The invention discloses a method for predicting the diffusion behavior of underground hydrogen in the pore scale of water-containing clay minerals and application, and belongs to the technical field of underground hydrogen energy storage and unconventional gas migration prediction. Performing energy minimization and pre-balance on the initial pore model of the water-containing clay mineral through a molecular dynamics method; filling hydrogen molecules under given temperature and pressure conditions through a giant regular Monte Carlo method, and correcting the hydrogen pressure through an SRK-EOS state equation; according to the method, the temperature and the fugacity are used as input parameters, hydrogen diffusion simulation is carried out through a molecular dynamics method, the mean square displacement of hydrogen in the pores is calculated, and the self-diffusion coefficient is obtained, and the diffusion rate change rule of the hydrogen in the clay pores under the conditions of different water film thicknesses, salinity and temperature pressure can be quantitatively represented; diffusion inflection points, diffusion suppression stages and migration limitation mechanisms caused by water film or water bridge structures are disclosed.

The invention discloses a method for calculating a capillary liquid bridge effect between particles based on an energy minimization principle, and relates to the technical field of wet particle system microcosmic acting force calculation. The method comprises the following steps: firstly, establishing a physical geometric model containing upper and lower spherical particles and a three-phase contact line, and calculating the total free energy of a system; then establishing volume conservation and contact angle constraint conditions, and solving a stable liquid bridge form which minimizes total free energy through iterative optimization algorithms such as gradient descent; and then capillary liquid bridge force is calculated based on the change of system energy along with the particle spacing in the stable form, and the liquid bridge fracture distance is determined by gradually increasing the spacing until a stable solution cannot be obtained. The method breaks through the idealized assumption of a traditional method on the shape and the contact angle of a liquid bridge, can accurately process the complex conditions of asymmetric and large-volume liquid bridges and the like, and provides accurate capillary force input parameters for discrete element simulation of a wet particle system.

The invention discloses a method for simulating the impact resistance behavior of a graphene-coated nanocrystalline metal composite structure based on molecular dynamics. The method is characterized by comprising the following steps: firstly, generating a nanocrystalline metal matrix model by adopting Atomsk software; then, LAMMPS software is used for constructing a target body composite model of the graphene coated nanocrystalline metal; performing energy minimization processing and relaxation operation on the target body composite model in sequence, adding a rigid projectile body into the target body composite model, constructing a complete impact simulation model, monitoring and recording motion parameters of the rigid projectile body in real time by LAMMPS software in the impact simulation process, and calculating the motion parameters of the rigid projectile body according to the motion parameters of the rigid projectile body. Visually analyzing the deformation, damage and failure evolution process of the target body composite model under the high-speed impact load through visual software; the method has the advantages that the anti-impact behavior and the damage failure mechanism of the graphene-coated nanocrystalline metal composite structure are accurately revealed from the micro-level of molecular atoms, the application range is wide, the method can adapt to different types of metal matrixes such as copper, iron and magnesium, and the anti-impact behaviors under different working conditions can be simulated by adjusting simulation parameters.

This invention discloses a method for simulating the growth of silica aerogel structures. The method includes the following steps: S1, establishing a molecular model of orthosilicic acid molecules and performing energy minimization calculations; S2, establishing an amorphous box; S3, performing energy minimization and molecular dynamics relaxation calculations on the amorphous box; S4, performing simulation condensation reaction calculations using a scripting language; S5, adjusting the cutoff radius. This method can quickly and accurately screen reactive atoms, avoid the occurrence of numerous binary ring structures, and yield more reasonable amorphous structures. The simulation process is simple and time-efficient.

The invention discloses a molecular dynamics simulation method for the impact resistance behavior of a metal-carbon honeycomb structure, and relates to the technical field of physical simulation, and the method comprises the steps: constructing a composite structure atom model and an impact projectile rigid body model; performing energy minimization and structural relaxation on the composite structure to obtain a stable initial configuration; the configuration is divided into a fixed constraint area, a temperature buffer area and a free response area along the in-plane direction of a structural plane, and fixed constraint, constant temperature control and micro-regularization system simulation are respectively applied; the composite structure is impacted by a rigid projectile at a preset speed, molecular dynamics impact simulation is carried out under a periodic boundary condition, and atomic motion data are output; by means of the data, structural deformation and damage morphology in the impact process are displayed through visual software, and by extracting a physical information curve, damage evolution and energy dissipation behaviors of a material are quantitatively analyzed, and the atomic scale impact resistance mechanism is revealed. The micro-mechanical behavior simulation of the metal-carbon honeycomb composite structure is realized.

The application discloses a polar surface liquid film thermal property prediction method based on molecular dynamics, relates to the technical field of microcosmic thermal property simulation, and comprises the following steps: constructing an initial simulation system containing a polar solid surface and polar liquid molecules, obtaining a stable system through temperature and pressure control and energy minimization, and running to thermodynamic equilibrium; extracting molecular trajectory data from the equilibrium system, calculating density distribution and molecular orientation distribution of a liquid film region, and obtaining microstructure characterization data; based on the same trajectory data, applying a Green-Kubo formula to calculate a heat flow autocorrelation function of the liquid film region, and obtaining a thermal conductivity prediction value of the region through time integral processing; and combining the microstructure data and the thermal conductivity prediction value, and analyzing an influence mechanism of a polar surface effect on liquid film thermal properties.

The invention belongs to the technical field of metal matrix composite modeling, and particularly relates to a particle reinforced metal matrix composite modeling method for molecular dynamics simulation, which comprises the following steps: establishing a metal matrix atomic model based on crystal structure information of a target matrix material; defining and removing a preset spherical region in the model to form a spherical cavity; constructing a spherical reinforced particle atomic model matched with the size of the cavity, and embedding the spherical reinforced particle atomic model into the corresponding cavity to generate a composite material initial model; and sequentially performing structural relaxation treatment on the initial model under the conditions of energy minimization and constant temperature and constant pressure to obtain a stable atomic model reaching a thermodynamic equilibrium state. By accurately controlling geometric parameters of the cavities and the particles, controllable design of the number, the position, the size and the volume ratio of the reinforced particles is realized, and the problem of atom overlapping is fundamentally avoided; in addition, the multi-particle uniform distribution modeling is supported, and the multi-scale combination design of reinforced particles with different sizes can be realized.

This application provides a molecular dynamics-based method and apparatus for modeling the self-crosslinking of epoxy resin, relating to the field of molecular dynamics. The method includes: constructing epoxy resin monomer molecules and labeling target atoms in the three-membered ring as reaction atom pairs; constructing an amorphous structure model based on the reaction atom pairs to form a molecular dynamics simulation system; entering a crosslinking cycle after energy minimization and kinetic equilibrium; searching and filtering reaction atom pairs with a first cutoff radius to obtain target reaction atom pairs; applying harmonic oscillator constraints to flexibly bring the reaction atoms closer and construct crosslinking chemical bonds, while simultaneously performing ring-opening operations on epoxy groups; calculating the current crosslinking conversion rate and comparing it with a preset crosslinking conversion rate threshold; if it is less than the threshold, expanding the cutoff radius to continue the crosslinking cycle; if it is greater than or equal to the threshold, performing abnormal coordination atom detection and repair. This application solves the problem that existing epoxy resin crosslinking simulation methods are difficult to realistically reproduce the actual crosslinking process of epoxy resin.

The application discloses a CO2-ion water flooding molecular dynamics simulation method and belongs to the technical field of CO2 oil displacement. The method comprises the following steps: establishing oil molecule models, water molecule models, CO2 molecule models, CaCl2 molecule models and graphene plate models corresponding to the characteristics of a target oil reservoir; optimizing the geometric structures of the models; assembling the optimized molecular models into a CO2-CaCl2-H2O ternary displacement system and performing molecular dynamics simulation; selecting an NVT ensemble during simulation, applying a periodic boundary condition, adopting a COMPASS II force field and configuring atomic potential parameters, running energy minimization and temperature initialization, and then performing displacement process simulation; quantitatively evaluating the oil displacement effect of the displacement medium of the CO2-CaCl2-H2O ternary system based on a trajectory file output by simulation; and the application can better explain the behavior of molecules under the synergistic action of multiple displacement media, reflect the oil displacement microscopic mechanism under the synergistic action of multiple displacement media, and has important significance for optimizing the development effect of tight oil reservoirs and improving the recovery rate.

The invention belongs to the field of cable defect data enhancement, and relates to a cable defect data enhancement method and system based on physical consistency constraint, and the system outputs a balanced training set and a test set through a data preprocessing module; the generative model construction module is used for embedding a double-field potential function in a potential space based on a conditional variation self-encoding network encoding observable vector and a conditional label, constructing a Hamiltonian energy function fusing the constraint of the double-field potential function, introducing an adaptive kinetic energy item, and unifying a PCDH-VAE generative model training target as a Hamiltonian energy minimization problem; a dynamic situation trap mechanism is introduced into the sample generation module, adaptive quality Hamiltonian Monte Carlo sampling is adopted, and a defect enhancement sample meeting multi-physics field constraint is generated through Hamiltonian dynamic evolution and Metropolis correction; and the sample evaluation module evaluates the generated defect enhancement sample set. The problem of scarcity of cable defect data samples can be relieved.

The application discloses a polar surface liquid film thermal property prediction method based on molecular dynamics, relates to the technical field of microcosmic thermal property simulation, and comprises the following steps: constructing an initial simulation system containing a polar solid surface and polar liquid molecules, obtaining a stable system through temperature and pressure control and energy minimization, and running to thermodynamic equilibrium; extracting molecular trajectory data from the equilibrium system, calculating density distribution and molecular orientation distribution of a liquid film region, and obtaining microstructure characterization data; based on the same trajectory data, applying a Green-Kubo formula to calculate a heat flow autocorrelation function of the liquid film region, and obtaining a thermal conductivity prediction value of the region through time integral processing; and combining the microstructure data and the thermal conductivity prediction value, and analyzing an influence mechanism of a polar surface effect on liquid film thermal properties.

The application provides a simulation method and system for diffusion behavior of nanoparticles in tissue gaps, and belongs to the field of molecular dynamics simulation, and comprises the following steps: constructing a multi-component coarse-grained molecular model; performing energy minimization on the multi-component coarse-grained model by using a conjugate gradient method, and obtaining an initial configuration with stable energy by setting a convergence tolerance; maintaining a constant temperature condition by using a Langevin heat bath under a canonical ensemble, performing molecular dynamics equilibrium simulation on the initial configuration to obtain the diffusion behavior of nanoparticles in tissue gaps, and recording motion trajectory data of the nanoparticles; obtaining mean square displacement (MSD) of the nanoparticles according to the motion trajectory data, and plotting MSD-t curves by corresponding MSD and simulation time t; performing linear fitting on the MSD-t curves and solving, to obtain a diffusion coefficient of the nanoparticles in the tissue gaps. The application analyzes the interaction between the nanoparticles and the solid-liquid two-phase medium from the molecular level, and provides a powerful microscopic research tool for elucidating the diffusion mechanism of the nanoparticles.

The application relates to a stereo matching method and device, computer equipment and a storage medium. The method comprises the following steps: performing grid processing on a to-be-matched image to obtain image blocks of the to-be-matched image, and for each image block, performing inverse search optimization operation on the gray values of each pixel point on the image block and the gray values of first candidate matching points corresponding to the pixel points on a reference image by using a pre-constructed inverse search optimization model to determine the first inverse depth value of the image block, and then determining the second inverse depth value corresponding to each pixel point according to the first inverse depth value of each image block, obtaining the depth map of the to-be-matched image by using a pre-constructed energy minimization model according to the second inverse depth value corresponding to each pixel point, and improving the calculation efficiency.

The present application relates to a kind of reef limestone mesostructure three-dimensional modeling method, comprising the following steps: scanning and obtaining the original gray image of core sample, and pre-processing;Based on the gray image after pre-processing, extract sample shell mask;Establish the Gaussian mixture model of sample gray distribution, determine the segmentation threshold;Through the classification segmentation of body data by segmentation threshold, obtain pore candidate area, matrix candidate area and undecided area;Morphological closing operation is carried out to candidate area, and the pixel set belonging to pore and the pixel set belonging to matrix are obtained;For undecided area, energy minimization algorithm based on graph cut is used for classification segmentation.The present application adopts GMM double-threshold segmentation and shell mask construction technology to accurately extract core area, avoid air at core boundary being misclassified as pore, so that the real boundary of core is completely retained;The present application can finely determine pore, and the automatic threshold segmentation result is stable and reliable, greatly reduces manual intervention.

The application relates to a prediction method and application of underground hydrogen adsorption behavior in a pore scale of a clay mineral containing impurities and water, and belongs to the technical field of underground hydrogen energy storage and unconventional energy development. The key points are as follows: an initial pore model of the clay mineral containing impurities and water is subjected to energy minimization and pre-equilibrium through a molecular dynamics method to obtain a stable pore model of the clay mineral containing impurities and water; based on the stable pore model of the clay mineral containing impurities and water, hydrogen fugacity is calculated through an SRK-EOS state equation, the corresponding relationship between the hydrogen fugacity and pressure is calibrated, and an initial hydrogen adsorption simulation model is obtained; based on the initial hydrogen adsorption simulation model, hydrogen molecules are introduced into the pore through a grand canonical Monte Carlo method with temperature and fugacity as input parameters, hydrogen insertion and deletion simulation is carried out, and a hydrogen adsorption model is obtained; and the hydrogen adsorption characteristics are calculated according to the hydrogen adsorption model, so that the prediction of the hydrogen adsorption behavior in the pore scale of an underground reservoir is realized.