Abstract：Dimensional discrete fracture network(DFN) is a common approach to analyze the structure of rock masses. However，Baecher disk method which is most popular，is entirely based on statistics and the assumption that fractures are disks rather than polygons for the convenience of mathematical derivation. In this study，a polygonal DFN modeling method more really reflecting fracture feature was developed to improve the traditional Baecher disk method. A circle-controlled method was proposed to determine the shape of random polygonal fractures and its validity was proven through mathematical derivation. For determining the size distribution of fractures，an iterative inversion algorithm was presented. To ensure that the algorithm is able to rapidly converge，a set of adjusting and initialization rules of the parameters of the algorithm were designed. A case study，in which comparison between the trace maps of two DFNs generated by the proposed method and Baecher method is carried out，shows that the iterative inversion algorithm can almost achieve the same results as Baecher method when assuming that the fractures are disks. Further，a polygonal discrete fracture network is generated using the presented method. This research provides a new way to analyze the structure of rock masses.
韩 帅1，李明超1，王 刚2. 基于控制圆法和迭代反演的岩体多边形裂隙网络模拟方法[J]. 岩石力学与工程学报, 2019, 38(8): 1635-1646.
HAN Shuai1，LI Mingchao1，WANG Gang2. A polygonal DFN modeling approach based on the circle-controlled method and the iterative inversion algorithm. , 2019, 38(8): 1635-1646.
[1] PRIEST S D，HUDSON J A. Discontinuity spacings in rock[J]. International Journal of Rock Mechanics and Mining Sciences and Geomechanics Abstracts，1976，13(5)：135–148.
[2] 汪小刚，贾志欣，张发明，等. 岩体结构面网络模拟原理及其工程应用[M]. 北京：中国水利水电出版社，2010：1–405.(WANG Xiaogang，JIA Zhixin，ZHANG Faming，et al. The simulation of rock joint network and its application[M]. Beijing：China Water Power Press，2010：1–405.(in Chinese))
[3] 薛秋池，赵其华，何云松. 岩体结构面网络模拟的改进与应用[J]. 岩土工程学报，2016，38(7)：1 351–1 356.(XUE Qiuchi，ZHAO Qihua，HE Yunsong. Improvement and application of network simulation of rock mass discontinuities[J]. Chinese Journal of Geotechnical Engineering，2016，38(7)：1 351–1 356.(in Chinese))
[4] 潘别桐，井兰如. 岩体结构概率模型模拟和应用[C]// 孙 钧，林韵梅，主编. 岩石力学新进展. 辽宁：东北工学院出版社，1989：55–80.(PAN Bietong，JING Lanru. Computer simulation methods and applications of statistical models of rockmass structure[C]// SUN Jun，LI Yunmei，ed. New Progress in Rock Mechanics. Liaoning：Press of Northeast Institute，1989：55–80.(in Chinese))
[5] 邬爱清，周火明，任 放. 岩体三维网络模拟技术及其在三峡工程中的应用[J]. 长江科学院院报，1998，15(6)：15–18.(WU Aiqing，ZHOU Huoming，REN Fang. Research on 3D rock joint network simulation techniques and its application to TGP[J]. Journal of Yangtze River Scientific Research Institute，1998，15(6)：15–18.(in Chinese))
[6] LI M C，ZHOU S B，WANG G. 3D identification and stability analysis of key surface blocks of rock slope[J]. Transactions of Tianjin University，2016，22(4)：317–323.
[7] VAN EYSINGA F W B. The minor structures of deformed Rocks—A photographic atlas[J]. Earth Science Reviews，1973，(9)：75–76.
[8] BAECHER G B. Statistical analysis of rock mass fracturing[J]. Journal of the International Association for Mathematical Geology，1983，15(2)：329–348.
[9] 张 文，陈剑平，苑晓青，等. 基于三维裂隙网络的RQD尺寸效应与空间效应的研究[J]. 岩石力学与工程学报，2012，31(7)：1 437– 1 445.(ZHANG Wen，CHEN Jianping，YUAN Xiaoqing，et al. Study of size effect and spatial effect of RQD for rock masses based on three-dimensional fracture network[J]. Chinese Journal of Rock Mechanics and Engineering，2012，31(7)：1 437–1 445.(in Chinese))
[10] VENEZIANO D. Probabilistic models of joints in rock[R]. Cambridge，MA：Department of Civil and Environmental Engineering，Massachusetts Institute of Technology，1978.
[11] DERSHOWITZ W S，EINSTEIN H H. Characterizing rock joint geometry with joint system models[J]. Rock Mechanics and Rock Engineering，1988，21(1)：21–51.
[12] IVANOVA V M，SOUSA R，MURRIHY B，et al. Mathematical algorithm development and parametric studies with the GEOFRAC three-dimensional stochastic model of natural rock fracture systems[J]. Computers and Geosciences，2014，67(3)：100–109.
[13] 宋晓晨，徐卫亚. 裂隙岩体渗流模拟的三维离散裂隙网络数值模型(I)：裂隙网络的随机生成[J]. 岩石力学与工程学报，2004，23(12)：2 015–2 020.(SONG Xiaochen，XU Weiya. Numerical model of three-dimensional discrete fracture network for seepage in fractured rocks (I)：Generation of fracture network[J]. Chinese Journal of Rock Mechanics and Engineering，2004，23(12)：2 015–2 020.(in Chinese))
[14] 郭 亮，李晓昭，周扬一，等. 随机与确定耦合的裂隙岩体结构面三维网络模拟[J]. 岩土力学，2016，37(9)：2 636–2 644.(GUO Liang，LI Xiaozhao，ZHOU Yangyi，et al. Three-dimensional network simulation of fractured rock mass discontinuities based on stochastic-deterministic coupling[J]. Rock and Soil Mechanics，2016，37(9)：2 636–2 644.(in Chinese))
[15] KULATILAKE P H S W，WU T H. Relation between discontinuity size and trace length[C]// The 27th US Symposium on Rock Mechanics(USRMS). [S. l.]：[s. n.]，1986：130–133.
[16] MAULDON M. Estimating mean fracture trace length and density from observations in convex windows[J]. Rock Mechanics and Rock Engineering，1998，31(4)：201–216.
[17] 张国强，邓建辉，费文平，等. 基于面样本的完整迹长分布推断方法初探[J]. 岩石力学与工程学报，2011，30(增1)：2 720–2 729. (ZHANG Guoqiang，DENG Jianhui，FEI Wenping，et al. Preliminary research on complete trace length distribution based on areal samples[J]. Chinese Journal of Rock Mechanics and Engineering，2011，30(Supp.1)：2 720–2 729.(in Chinese))
[18] ZHANG L，EINSTEIN H H. Estimating the intensity of rock discontinuities[J]. International Journal of Rock Mechanics and Mining Sciences，2000，37(5)：819–837.
[19] ZHANG L，DERSHOWITZ W S，EINSTEIN H H. Stereological relationship between trace length and size distribution of elliptical discontinuities[J]. Geotechnique，2002，52(6)：419–434.
[20] TONON F，CHEN S. Closed-form and numerical solutions for the probability distribution function of fracture diameters[J]. International Journal of Rock Mechanics and Mining Sciences，2007，44(3)：332–350.
[21] 张 奇. 岩体随机结构面尺寸概率估算及三维网络模型应用研究[硕士学位论文][D]. 长春：吉林大学，2015.(ZHANG Qi. Probabilistic estimation of discontinuity size in rock mass and study of application for three-dimensional network model-taking dam abutment rock mass in Songta Hydroelectric station as an example[M. S. Thesis][D]. Changchun：Jilin University，2015.(in Chinese))
[22] ROBERTSON A M. The interpretation of geological factors for use in slope theory[C]// VAN RENSBURG P W J，ed. Proc. Symp. Planning Open Pit Mines. Johannesburg，South Africa：[s. n.]，1970：55–71.
[23] XU C，DOWD P. A new computer code for discrete fracture network modelling[J]. Computers and Geosciences，2010，36(3)：292–301.
[24] DERSHOWITZ W S，EINSTEIN H H. Characterizing rock joint geometry with joint system models[J]. Rock Mechanics and Rock Engineering，1988，21(1)：21–51.
[25] IVANOVA V M，YU X，VENEZIANO D，et al. Development of stochastic models for fracture systems[C]// SCHULTS D ed. The 35th US Symposium on Rock Mechanics(USRMS). [S. l.]：[s. n.]，1995：725–730.
[26] NOROOZI M，KAKAIE R，JALALI S E. 3D Geometrical- Stochastical modeling of rock mass joint networks：Case study of the right bank of rudbar lorestan dam plant[J]. Journal of Geology and Mining Research，2015，7(1)：1–10.
[27] PRATT J W. Concavity of the Log likelihood[J]. Publications of the American Statistical Association，1981，76(373)：103–106.
[28] 郭 亮，李晓昭，周扬一，等. 数学字结构面三维随机网络模拟及其有效性检验[J]. 岩石力学与工程学报，2015，34(Supp. 1)：2 854–2 861.(GUO Liang，LI Xiaozhao，ZHOU Yangyi，et al. Simulation of random 3D discontinuities network based on digitalization and its validation test[J]. Chinese Journal of Rock Mechanics and Engineering，2015，34(Supp.1)：2 854–2 861.(in Chinese)).
[29] HAN S，WANG G，LI M C. A trace map comparison algorithm for the discrete fracture network models of rock masses[J]. Computers and Geosciences，2018，115：31–41.
2019, Vol. 38 
