Abstract：Equivalent construction method in statistics of coal?s pore structure was discussed；and the spatial- temporal evolution process of fluid flow was simulated at pore scale. First of all，the porous media model of coal medium was constructed using modified quartet structure generation set(QSGS) algorithm which integrates the nonlinear and anisotropic characteristics of coal?s pore structure. And then，single-phase gas flow through pores was simulated using lattice Boltzmann method. Based on the numerical simulation results，the main influential factors and the evolution pattern of coal?s permeability before fluid flow achieve equilibrium state were analyzed. The analysis results show that：(1) Permeability in porous media is dominated by the channels connected by pores with large size，but the migration in the paths connected by small pore and micropore was concentration diffusion behavior；(2) In the same spatial distribution variability pattern，there is a power law relation between porosity and permeability，but the power rate coefficient increases with the decrease of spatial variability；(3) Permeability is reduced gradually over time in porous media and the time needed before the flow reaches its stable state meets negative correlation power rate relation with porosity.
[1] 李仁民，刘松玉，方 磊，等. 采用随机生长四参数生成法构造黏土微观结构[J]. 浙江大学学报：工学版，2010，44(10)：1 897–1 901. (LI Renmin，LIU Songyu，FANG Lei，et al. Micro-structure of clay generated by quartet structure generation set[J]. Journal of Zhejiang University：Engineering Science，2010，44(10)：1 897–1 901.(in Chinese))
[2] 王唯威. 分形多孔介质内导热与流动数值模拟研究[硕士学位论文][D]. 北京：中国科学院工程热物理研究所，2006：19–32. (WANG Weiwei. Numerical study of heat conduction and fluid flow in fractal porous media[M. S. Thesis][D]. Beijing：Institute of Engineering Thermophysics，Chinese Academy of Sciences，2006：19–32.(in Chinese))
[3] TACHER L，PERROCHET P，PARRIAUX A. Generation of granular media[J]. Transport in Porous Media，1997，26(1)：99–107.
[4] PILOTTI M. Generation of realistic porous media by grains sedimentation[J]. Transport in Porous Media，1998，33(3)：257–278.
[5] 刘勇军，朱岳明，丰土根，等. 随机土体结构模型理论及应用[J]. 岩土工程学报，2001，23(3)：354–357.(LIU Yongjun，ZHU Yueming，FENG Tugen，et al. Study of random soil frame model and its application[J]. Chinese Jounal of Geotechnical Engineering，2001，23(3)：354–357.(in Chinese))
[6] WANG M，PAN N. Numerical analyses of effective dielectric constant of multiphase microporous media[J]. Journal of Applied Physics，2007，101(11)：114102-1–114120-8.
[7] NARSILIO G A，BUZZI O，FITYUS S，et al. Upscaling of Navier-Stokes equations in porous media：theoretical，numerical and experimental approach[J]. Computers and Geotechnics，2009，36(7)：1 200–1 206.
[8] YOUNGSEUK K. Computational rock physics：transport properties in porous media and applications[Ph. D. Thesis][D]. Stanford：Stanford University，2003.
[9] HARDY J，POMEAU Y，，PAZZIS O D. Time evolution of a two-dimensional classical lattice system[J]. Physics Review Letters，1973，31(5)：276–279.
[10] HARDY J，POMEAU Y，PAZZIS O D. Time evolution of a two-dimensional model system. I. Invariant states and time correlation functions[J]. Journal of Mathematical Physics，1973，14(12)：1 746–1 759.
[11] FRISCH U，HASSLACHER B，POMEAU Y. Lattice-gas automata for the Navier-Stokes equation[J]. Physics Review Letters，1986，56(14)：1 505–1 508.
[12] WOLFRAM S. Cellular automaton fluids 1：Basic theory[J]. Journal of Statistical Physics，1986，45(3/4)：471–526.
[13] FRISCH U，D'HUMIÈRES D，HASSLACHER B，et al. Lattice gas hydrodynamics in two and three dimensions[J]. Complex Systems，1987，1：649–707.
[14] D'HUMIÈRES D，LALLEMAND P，FRISCH U. Lattice gas models for 3D hydrodynamics[J]. Europhysics Letters，1986，2(4)：291–297.
[15] ROTHMAN D H，KELLER J M. Immiscible cellular-automaton fluids[J]. Journal of Statistical Physics，1988，52(3/4)：1 119–1 127.
[16] BENZI R，SUCCI S，VERGASSOLA M. The lattice Boltzmann equation：theory and applications[J]. Physics Reports，1992，222(3)：145–197.
[17] MCNAMARA G R，ZANETTI G. Use of the Boltzmann equation to simulate lattice-gas automata[J]. Physics Review Letters，1988，61(20)：2 332–2 335.
[18] HIGUERA F J，SUCCI S，BENZI R. Lattice gas-dynamics with enhanced collisions[J]. Europhysics Letters，1989，9(4)：345–349.
[19] CHEN S Y，CHEN H D，MARTINEZ D，et al. Lattice Boltzmann Model for simulation of magneto-hydrodynamics[J]. Physical Review Letters，1991，67(27)：3 776–3 779.
[20] KOELMAN J M V A. A simple lattice Boltzmann scheme for Navier-Stokes fluid-flow[J]. Europhysics Letters，1991，15(6)：603–607.
[21] QIAN Y H，D'HUMIÈRES D，LALLEMAND P. Lattice BGK models for Navier-Stokes equation[J]. Europhysics Letters，1992，17(6)：479–484.
[22] ABE T. Derivation of the lattice Boltzmann method by means of the discrete ordinate method for the Boltzmann equation[J]. Journal of Computational Physics，1997，131(1)：241–246.
[23] HE X Y，LUO L S. Theory of the lattice Boltzmann method：from the Boltzmann equation to the lattice Boltzmann equation[J]. Physical Review E，1997，56(6)：6 811–6 817.
[24] SHAN X W，HE X Y. Discrimination of the velocity space in the solution of the Boltzmann equation[J]. Physical Review Letters，1998，80(1)：65–68.
[25] 滕桂荣，谭云亮，高 明，等. 基于LBM方法的裂隙煤体内瓦斯抽放的模拟分析[J]. 煤炭学报，2008，33(8)：914–919.(TENG Guirong，TAN Yunliang，GAO Ming，et al. Simulation of gas drainage in fissured coal based on lattice Boltzmann method[J]. Journal of China Coal Society，2008，33(8)：914–919.(in Chinese))
[26] 滕桂荣，谭云亮，高 明. 基于Lattice Boltzmann方法对裂隙煤体中瓦斯运移规律的模拟研究[J]. 岩石力学与工程学报，2007，26(增1)：3 503–3 508.(TENG Guirong，TAN Yunliang，GAO Ming，et al. Simulation of gas seepage in fissured coal based on Lattice Boltzmann method[J]. Chinese Journal of Rock Mechanics and Engineering，2007，26(Supp.1)：3 503–3 508.(in Chinese))
[27] 朱益华，陶 果，方 伟，等. 低渗气藏中气体渗流Klinkenberg效应研究进展[J]. 地球物理学进展，2007，22(5)：1 591–1 596. (ZHU Yihua，TAO Guo，FANG Wei，et al. Research progress of the Klinkenberg effect in tight gas reservoir[J]. Progress in Geophysics，2007，22(5)：1 591–1 596.(in Chinese))
[28] 邢玉飞. 基于LBM地下煤层瓦斯运移仿真实现[硕士学位论文][D]. 西安：西安建筑科技大学，2009.(XING Yufei. Simulation of gas flow in underground coal seam based on LBM[M. S. Thesis][D]. Xi?an：Xi?an University of Architecture and Technology，2009.(in Chinese))
[29] 陆秋琴. 地下煤矿瓦斯运移数值模拟及积聚危险性评价研究[博士学位论文][D]. 西安：西安建筑科技大学，2010.(LU Qiuqin. Research on numerical simulation of gas migration and risk evaluation of gas aggregation in underground coal mine[Ph. D. Thesis][D]. Xi?an：Xi?an University of Architecture and Technology，2010.(in Chinese))
[30] 傅雪海，秦 勇，薛秀谦，等. 煤储层孔、裂隙系统分形研究[J]. 中国矿业大学学报，2001，30(3)：225–228.(FU Xuehai，QIN Yong，XUE Xiuqian，et al. Research on fractals of pore and fracture- structure of coal reservoirs[J]. Journal of China University of Mining and Technology，2001，30(3)：225–228.(in Chinese))
[31] FU X H，QIN Y，ZHANG W H，et al. Fractal classification and natural classification of coal pore structure based on migration of coal bed methane[J]. Chinese Science Bulletin，2005，50(Supp.1)：66–71.
[32] FRIESEN W I，MIKULA R J. Fractal dimensions of coal particles[J]. Journal of Colloid and Interface Science，1987，120(1)：263–271.
[33] 郭照立，郑楚光. 格子Boltzmann方法的原理及应用[M]. 北京：科学出版社，2008：17–21.(GUO Zhaoli，ZHENG Chuguang. Theory and applications of lattice Boltzmann method[M]. Beijing：Science Press，2008：17–21.(in Chinese))
[34] SUCCI S. The lattice Boltzmann equation for fluid dynamics and beyond[M]. New York：Oxford University Press，2001：304.
2013, Vol. 32 
