基于幂律规则的滨海矿山巷道涌水混合比分析

doi:10.11872/j.issn.1005-2518.2021.01.116

黄金科学技术 ›› 2021, Vol. 29 ›› Issue (1): 99-107.doi: 10.11872/j.issn.1005-2518.2021.01.116

基于幂律规则的滨海矿山巷道涌水混合比分析

段学良^1,^2,³(),马凤山^1,²(),郭捷^1,²,惠鑫⁴,顾鸿宇⁵,王善飞⁶

^1.中国科学院地质与地球物理研究所，中国科学院页岩气与地质工程重点实验室，北京 100029
^2.中国科学院地球科学研究院，北京 100029
^3.中国科学院大学，北京 100049
^4.北京市基础设施投资有限公司，北京 100101
^5.中国地质调查局成都地质调查中心，四川成都 610081
^6.山东黄金矿业（莱州）有限公司三山岛金矿，山东莱州 261442

收稿日期:2020-06-30 修回日期:2020-07-20 出版日期:2021-02-28 发布日期:2021-03-22
通讯作者: 马凤山 E-mail:13051876966@163.com;fsma@mail.iggcas.ac.cn
作者简介:段学良（1994-），男，河北泊头人，博士研究生，从事矿山水文地质与工程地质研究工作。13051876966@163.com
基金资助:
国家自然科学基金重点项目“海底采矿对地质环境的胁迫影响与致灾机理”(41831293);国家自然科学基金青年项目“基于‘盲源分离’的矿山突水水源识别与混合比研究”(41907174);国家重点研发计划专题“黄渤海不同类型海岸带海水入侵发生机理研究”(2016YFC0402802-01)

Analysis of Mine Water Mixing Ratio in a Coastal Deposit Based on Power Law

Xueliang DUAN^1,^2,³(),Fengshan MA^1,²(),Jie GUO^1,²,Xin HUI⁴,Hongyu GU⁵,Shanfei WANG⁶

^1.Key Laboratory of Shale Gas and Geoengineering，Institute of Geology and Geophysics，Chinese Academy of Sciences，Beijing 100029，China
^2.Innovation Academy for Earth Science，Chinese Academy of Sciences，Beijing 100029，China
^3.University of Chinese Academy of Sciences，Beijing 100049，China
^4.Beijing Infrastructure Investment Co. ，Ltd. ，Beijing 100101，China
^5.Chengdu Center，China Geological Survey，Chengdu 610081，Sichuan，China
^6.Sanshandao Gold Mine，Shandong Gold Mining（Laizhou）Company Limited，Laizhou 261442，Shandong，China

Received:2020-06-30 Revised:2020-07-20 Online:2021-02-28 Published:2021-03-22
Contact: Fengshan MA E-mail:13051876966@163.com;fsma@mail.iggcas.ac.cn

摘要/Abstract

摘要：

滨海采矿引起的突涌水灾害对矿山的安全生产构成了极大威胁，计算巷道涌水的混合比，进一步分析其演化规律对于突涌水事故的防治具有重要意义。将幂律规则应用于三山岛金矿巷道涌水混合比数据的统计分析中，根据2个已有的混合比研究结果，采用概率密度函数p（S）对相邻2个监测周期的海水比例波动事件进行了拟合，拟合的相关系数分别达到0.92和0.93，说明海水比例波动值区间间隔的均值与概率密度的分布符合幂律分布。因此，采用幂律规则对研究区的涌水混合比进行统计分析是可信的。研究结果表明：不同方法得出的混合比，在幂律规则下反映出的规律是相同的；海水比例波动值大于48%的概率小于5%，因此将48%视为预警区间的临界值；F₃断层附近的涌水点相对其他监测点海水比值波动较大，这是由于F₃为导水断层且连通了海水，受采动的影响，F₃断层周围的导水通道错综复杂且不稳定。

关键词: 滨海采矿, 突涌水灾害, 混合比, 演化规律, 统计分析, 幂律规则

Abstract:

The study area，Sanshandao gold mine，is the first coastal mine in China.It belongs to structural fissure water-filled mine，and the hydrogeological conditions are complicated.With the mining of the orebody，multiple water inrush accidents occurred in the shallow and deep parts of the mine，causing the partial roadway to be flooded，and even some sections were accompanied by sand erosion.Water inrush caused by coastal mining is an extreme threat to the safe production of the mine.Therefore，it is important to determine the mine water mixing ratios and analyze its evolutionary law for the prevention of water inrush accidents.The power-law rule is a general law shown in the occurrence of geological disasters in nature.It refers to the relationship between the frequency and the scale of disasters.The frequency of large-scale disasters is low.Conversely，disasters with a high frequency of occurrence are relatively small in scale.To determine the measure of the proportion of seawater，the power-law rule was applied to the statistical analysis of the mixing ratios of the mine water in this study area.Firstly，the results of two existing mixing ratio studies were statistically sorted out，and the probability density statistical results of seawater fluctuation events were obtained.Then，the probability density function p（S） was used to fit the fluctuation events of the seawater ratio in two adjacent monitoring periods.Finally，by integrating the fitted curve，the early warning interval of seawater fluctuation value was obtained.The research results show that the correlation coefficients of the fitting reach 0.92 and 0.93，respectively.It indicates that the distribution of the interval mean value and probability density of seawater proportional fluctuation events conformed to the power-law distribution.Thus，it is credible to use the power-law rule to analyze the mixing ratio of the mine water in the study area.For the mixing ratios obtained by different methods，the law reflected under the power law rule is the same.The fluctuation values of the seawater ratio at most monitoring sites are not large，less than 30%.It shows that the power-law rule is not affected by the calculation method of mixing ratio.Because the selected analysis index is the seawater fluctuation value，that is，for the relative value of two monitoring periods，the errors caused by different methods are eliminated.The probability that the fluctuation value of seawater ratio is greater than 48% is less than 5%，so 48% is regarded as the critical value of the warning interval.When the seawater fluctuation value is greater than this value，it should be paid attention to，and combined with the water temperature，flow rate，and other indicators of the water site for further analysis.Water samples near the F₃ fault have larger fluctuation values of seawater than that of other water samples because F₃ connects the seawater，and due to mining，the water channels around F₃ are complicated and unstable.

Key words: coastal mining, water inrush, mixing ratio, evolution law, statistical analysis, power law

中图分类号:

P642

段学良,马凤山,郭捷,惠鑫,顾鸿宇,王善飞. 基于幂律规则的滨海矿山巷道涌水混合比分析[J]. 黄金科学技术, 2021, 29(1): 99-107.

Xueliang DUAN,Fengshan MA,Jie GUO,Xin HUI,Hongyu GU,Shanfei WANG. Analysis of Mine Water Mixing Ratio in a Coastal Deposit Based on Power Law[J]. Gold Science and Technology, 2021, 29(1): 99-107.

图/表 5

图1

表1

涌水样混合比计算结果（Duan et al.，2019）"

样品编号	海水	卤水	淡水	样品编号	海水	卤水	淡水	样品编号	海水	卤水	淡水
105-1a	0.52	0.21	0.26	375-7'f	0.81	0.07	0.12	555-5b	0.53	0.39	0.08
105-1b	0.64	0.08	0.28	375-8a	0.62	0.20	0.18	555-5c	0.61	0.28	0.11
150-1a	0.16	0.29	0.54	375-8b	0.56	0.31	0.14	555-6a	0.47	0.42	0.11
150-1b	0.30	0.19	0.51	375-8'e	0.80	0.14	0.05	555-6b	0.58	0.36	0.05
150-1c	0.23	0.22	0.55	375-8'f	0.77	0.09	0.14	555-7b	0.83	0.05	0.12
150-1d	0.35	0.29	0.35	375-9e	0.80	0.14	0.07	600-1a	0.45	0.43	0.12
195-1a	0.57	0.16	0.28	375-9f	0.83	0.07	0.10	600-1b	0.54	0.34	0.12
195-1b	0.62	0.15	0.24	375-10f	0.63	0.21	0.16	600-1c	0.61	0.30	0.09
195-1c	0.52	0.19	0.29	375-11f	0.48	0.50	0.02	600-1d	0.56	0.24	0.20
195-1d	0.52	0.20	0.29	375-12f	0.45	0.49	0.06	600-1e	0.45	0.41	0.14
240-1a	0.62	0.28	0.10	375-13f	0.31	0.63	0.06	600-2a	0.50	0.39	0.11
240-2a	0.70	0.18	0.12	420-1a	0.62	0.19	0.19	600-2b	0.63	0.31	0.06
240-3a	0.34	0.58	0.08	420-1c	0.66	0.22	0.13	600-2c	0.60	0.32	0.08
240-4a	0.72	0.16	0.12	420-2a	0.53	0.33	0.14	600-2d	0.60	0.20	0.20
240-5a	0.60	0.10	0.31	465-1a	0.54	0.34	0.12	600-3a	0.38	0.54	0.08
285-1a	0.35	0.42	0.24	465-1b	0.48	0.35	0.17	600-3b	0.50	0.44	0.07
285-1b	0.58	0.27	0.15	465-1c	0.38	0.43	0.20	600-3c	0.56	0.37	0.07
285-1c	0.57	0.29	0.15	465-2a	0.36	0.41	0.23	600-4b	0.62	0.30	0.07
285-1d	0.38	0.35	0.28	465-2c	0.41	0.37	0.22	600-4c	0.56	0.39	0.05
285-2a	0.61	0.17	0.21	510-1a	0.74	0.16	0.10	600-5b	0.45	0.46	0.09
285-2b	0.75	0.11	0.15	510-1b	0.82	0.05	0.13	600-5c	0.53	0.35	0.12
285-2d	0.45	0.33	0.22	510-1c	0.71	0.18	0.11	600-5d	0.35	0.46	0.19
285-3a	0.00	0.54	0.46	510-2a	0.79	0.10	0.11	600-6c	0.67	0.22	0.10
285-3b	0.00	0.50	0.50	510-2b	0.87	0.03	0.10	600-6d	0.12	0.62	0.26
285-3c	0.00	0.51	0.49	510-2c	0.73	0.10	0.17	600-6e	0.79	0.11	0.11
285-3d	0.00	0.58	0.42	510-3a	0.72	0.14	0.14	600-7c	0.61	0.31	0.07
320-7e	0.11	0.87	0.03	510-3b	0.74	0.12	0.14	600-8c	0.56	0.39	0.05
320-8e	0.24	0.76	0.00	510-4a	0.70	0.16	0.14	600-8d	0.57	0.27	0.16
320-9e	0.39	0.61	0.00	510-4b	0.80	0.07	0.13	600-9c	0.58	0.37	0.05
330-1a	0.55	0.31	0.14	510-5a	0.71	0.13	0.16	600-9d	0.55	0.37	0.08
330-1b	0.54	0.33	0.13	510-6a	0.72	0.10	0.18	600-10d	0.71	0.11	0.18
330-2a	0.59	0.17	0.24	510-6b	0.80	0.07	0.13	600-11e	0.00	0.81	0.19
375-1a	0.67	0.21	0.12	510-7a	0.40	0.48	0.12	600-12e	0.31	0.54	0.15
375-1b	0.68	0.17	0.15	510-8c	0.76	0.07	0.16	600-13e	0.61	0.28	0.10
375-1c	0.64	0.21	0.15	510-9d	0.64	0.13	0.23	600-14e	0.69	0.17	0.13
375-1e	0.75	0.14	0.11	510-10d	0.61	0.15	0.24	600-15e	0.82	0.07	0.11
375-2a	0.70	0.20	0.11	510-11f	0.32	0.63	0.05	600-16e	0.78	0.07	0.15
375-3a	0.78	0.10	0.13	510-12f	0.60	0.29	0.11	600-17e	0.76	0.07	0.17
375-3b	0.78	0.08	0.14	510-13f	0.64	0.28	0.09	600-17f	0.53	0.25	0.22
375-3d	0.59	0.19	0.22	510-15f	0.62	0.28	0.10	600-18e	0.81	0.04	0.14
375-4a	0.75	0.11	0.14	510-16f	0.54	0.35	0.11	600-18f	0.71	0.11	0.18
375-4b	0.71	0.12	0.17	510-16kf	0.71	0.21	0.08	600-19e	0.82	0.03	0.14
375-4c	0.69	0.18	0.13	555-1a	0.72	0.11	0.18	600-20f	0.62	0.23	0.15
375-4d	0.55	0.24	0.21	555-1c	0.71	0.15	0.15	600-21f	0.47	0.32	0.21
375-4e	0.71	0.17	0.12	555-1d	0.22	0.47	0.31	600-22f	0.77	0.08	0.15
375-5a	0.47	0.34	0.19	555-2a	0.72	0.12	0.16	600-23f	0.79	0.06	0.14
375-5b	0.52	0.30	0.18	555-3a	0.67	0.14	0.19	645-1f	0.58	0.32	0.10
375-5c	0.45	0.34	0.21	555-3b	0.67	0.17	0.16	645-2f	0.38	0.39	0.23
375-5d	0.35	0.39	0.27	555-3c	0.42	0.39	0.19	645-3f	0.54	0.34	0.12
375-5e	0.57	0.28	0.15	555-3d	0.40	0.35	0.25	645-4f	0.49	0.40	0.11
375-6a	0.51	0.35	0.14	555-4a	0.67	0.18	0.14	690-1f	0.55	0.35	0.10
375-6b	0.47	0.37	0.16	555-4b	0.80	0.09	0.11	690-10f	0.49	0.27	0.24
375-6'f	0.78	0.08	0.13	555-4c	0.76	0.11	0.13	690-2c	0.59	0.38	0.03
375-7a	0.49	0.33	0.18	555-4d	0.68	0.12	0.20	690-2f	0.63	0.26	0.10
375-7'e	0.79	0.10	0.11	555-5a	0.68	0.18	0.14	690-12f	0.67	0.24	0.09

表1

表2

表3

图2

参考文献 0

	Bak P，Tang C，1989.Earthquakes as a self-organized critical phenomenon［J］.Journal of Geophysical Research，94（B11）：15635-15637.
	Bak P，Tang C，Wiesenfeld K，1987.Self-organized criticality：An explanation of 1/f noise［J］.Physical Review Letters，59（4）：381-384.
	Carlson J M，Langer J S，1989.Mechanical model of an earthquake fault［J］.Physical Review A，40：6470-6483.
	Chung J S，1996.Deep-ocean mining：Technologies for manganese nodules and crusts［J］.International Journal of Offshore Polar Engineering，6（4）：244-254.
	Commeau R F，Clark A，Johnson C，al et，1984.Ferromanganese crust resources in the Pacific and Atlantic Oceans［C］//Proceedings of the Oceans Conference （IEEE）.New York：IEEE.
	Duan X L，Ma F S，Guo J，al et，2019.Source identification and quantification of seepage water in a coastal mine，in China［J］.Water，11（9）：1862.
	Dussauge C，Grasso J R，Helmstetter A，2003.Statistical analysis of rockfall volume distributions： Implications for rockfall dynamics［J］.Journal of Geophysical Research Solid Earth，108（B6）：1-11.
	Gao Song，Zhang Junjin，Sun Shanshan，al et，2016.Hydrogeological characteristics of gold deposit in north sea area of Sanshandao［J］.Gold Science and Technology，24（1）：11-16.
	Gu H Y，Ma F S，Guo J，al et，2017.Hydrochemistry，multidimensional statistics，and rock mechanics investigations for Sanshandao gold mine，China［J］.Arabian Journal of Geosciences，10（3）：62.
	Gu H Y，Ma F S，Guo J，al et，2018.Assessment of water sources and mixing of groundwater in a coastal mine：The Sanshandao gold mine，China［J］.Mine Water and the Environment，37：351-365.
	Guo J，Zhao H J，Ma F S，al et，2015.Investigating the permeability of fractured rock masses and the origin of water in a mine tunnel in Shandong Province，China［J］.Water Science and Technology，72（11）：2006-2017.
	Guzzetti F，Malamud B D，Turcotte D L，al et，2002.Power-law correlations of landslide areas in central Italy［J］.Earth and Planetary Science Letters，195（3/4）：169-183.
	Hui X，Ma F S，Guo J，al et，2018.Power-law correlations of mine subsidence at a metal mine in China［J］.Environmental Geotechnics：1-12.doi：10.1680/jenge.18.00039. doi: 10.1680/jenge.18.00039
	Hui X，Ma F S，Zhao H J，al et，2019.Monitoring and statistical analysis of mine subsidence at three metal mines in China［J］.Bulletin of Engineering Geology & the Environment，78（6）：3983-4001.
	Hurst M D，Ellis M A，Royse K R，al et，2013.Controls on the magnitude-frequency scaling of an inventory of secular landslides［J］.Earth Surface Dynamics，1（1）：67-78.
	Ito K，Tsuzaki M，1990.Earthquakes as self-organized critical phenomena［J］.Journal of Geophysical Research，95（B5）：6853.
	Liu Yushan，Wu Bihao，2005.Exploitation of marine mineral resources：Review and prospects［J］.Mineral Deposits，24（1）：81-84.
	Ma F S，Zhao H J，Guo J，2015.Investigating the characteristics of mine water in a subsea mine using groundwater geochemistry and stable isotopes［J］.Environmental Earth Sciences，74（9）：6703-6715.
	Nakanishi H，1990.Cellular-automaton model of earthquakes with deterministic dynamics［J］.Physical Review A，41（12）：7086-7089.
	Peng K，Li X B，Wang Z W，2015.Hydrochemical characteristics of groundwater movement and evolution in the Xinli deposit of the Sanshandao gold mine using FCM and PCA methods［J］.Environmental Earth Sciences，73（12）：7873-7888.
	Qiu Haijun，Cao Mingming，Liu Wen，2013.Power-law correlations of landslides：A case of Ningqiang County［J］.Geolo-gical Science and Technology Information，32（3）：183-187.
	Rona P A，2003.Resources of the sea floor［J］.Science，299（5607）：673-674.
	Teixeira S B，2006.Slope mass movements on rocky sea-cliffs：A powerlaw distributed natural hazard on the Barlavento Coast，Algarve，Portugal［J］.Continental Shelf Research，26（9）：1077-1091.
	Wang Shanfei，2001.Analysis of hydrogeology for deep mining in Sanshandao gold mine［J］.Nonferrous Mines，30（3）：9-12.
	Xu Qiang，Huang Runqiu，1997.Power law between volume and frequency of geological hazards［J］.Journal of Chengdu University of Technology，24（Supp.1）：93-98.
	Ye Bailong，Peng Ensheng，1994.Study on the conducting-water structure model in the Sanshandao gold deposit［J］.Journal of Central South Universtiy（Science and Technology），25（2）：146-150.
	Zhang Shouquan，Huang Wei，1994.Hydrogeological and engineering features of F₃ fracture zone in Sanshan Island gold mine district and prevention of hazards ［J］.Journal of Engineering Geology，2（1）：62-72.
	高松，张军进，孙珊珊，等，2016.三山岛北部海域金矿区水文地质特征分析［J］.黄金科学技术，24（1）：11-16.
	刘玉山，吴必豪，2005.海底金属矿产资源的开发——回顾与未来展望［J］.矿床地质，24（1）：81-84.
	邱海军，曹明明，刘闻，2013.地质灾害的幂律相依性：以宁强县为例［J］.地质科技情报，32（3）：183-187.
	王善飞，2001.三山岛金矿深部开采水文地质浅析［J］.有色矿山，30（3）：9-12.
	许强，黄润秋，1997.地质灾害发生频率的幂律规则［J］.成都理工学院学报，24（增1）：93-98.
	叶柏龙，彭恩生，1994.三山岛金矿导水构造模式研究［J］.中南大学学报（自然科学版），25（2）：146-150.
	张寿全，黄巍，1994.三山岛金矿F₃断裂带的水文地质工程地质特征及灾害防治［J］.工程地质学报，2（1）：62-72.

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed

统计区间	区间均值	事件个数	累计个数	发生频率	概率密度
[0，10）	5	45	45	0.625	0.063
[10，20）	15	19	64	0.264	0.026
[20，30）	25	5	69	0.069	0.007
[30，40）	35	2	71	0.028	0.003
[40，50）	45	1	72	0.014	0.001

统计区间	区间均值	事件个数	累计个数	发生频率	概率密度
[0，10）	5	46	46	0.613	0.061
[10，20）	15	22	68	0.293	0.029
[20，30）	25	4	72	0.053	0.005
[40，50）	45	1	73	0.013	0.001
[50，60）	55	1	74	0.013	0.001
[60，70）	65	1	75	0.013	0.001

[1]	鞠培姣. 太行山北段梨园金矿构造地球化学特征[J]. 黄金科学技术, 2021, 29(1): 64-73.
[2]	崔宇,李夕兵,董陇军,白吕. 玲珑金矿微震监测台网布设优化[J]. 黄金科学技术, 2019, 27(3): 417-424.
[3]	段学良,马凤山,赵海军,郭捷,顾鸿宇,刘帅奇. 滨海矿山矿坑涌水源识别与混合比研究[J]. 黄金科学技术, 2019, 27(3): 406-416.
[4]	殷大发,孟祥瑞,殷大澍. 矿井采掘计划及智能决策支持系统[J]. J4, 2004, 12(6): 38-42.

基于幂律规则的滨海矿山巷道涌水混合比分析

Analysis of Mine Water Mixing Ratio in a Coastal Deposit Based on Power Law

RichHTML

PDF (PC)

摘要/Abstract

引用本文

使用本文

图/表 5

参考文献 0

相关文章 4

Metrics

本文评价

推荐阅读 0