Distribution and genesis of global Na-carbonate deposits and its prospecting potential
-
摘要:研究目的
天然碱矿主要用于制纯碱,和合成碱相比,具有绿色环保、成本低的优势。其下游产业遍布各个领域,如玻璃、医药等。光伏等新能源产业的蓬勃发展给纯碱带来了新需求。中国是纯碱消费大国,天然碱或成为稀缺资源。天然碱矿是一种蒸发岩矿床,成矿过程具有重要的研究意义,国内蒸发岩矿床研究领域对天然碱矿的关注较少。理清世界天然碱矿产资源类型、分布特征,总结碱矿成因对深化中国天然碱矿成矿规律研究、预测找矿远景、寻找天然碱矿床具有积极意义。
研究方法搜集已公开发表或出版的天然碱矿床资料,对世界典型矿床区域概况、物源、成因等内容进行系统总结。
研究结果绝大部分天然碱矿床分布于北美洲、亚洲和非洲,主要形成于新生代。类型以现代天然碱矿床居多,古代天然碱矿床较少,但资源规模大。古代天然碱矿以美国绿河组、河南泌阳凹陷安棚天然碱矿最为典型,现代天然碱矿以肯尼亚马加迪湖和内蒙古查干诺尔泡碱矿床为典型。
结论天然碱矿需要在封闭的构造环境、干旱的气候条件、充足的碳酸钠型卤水补给等成矿因素耦合条件成矿。与其他蒸发岩不同的是,维持碳酸钠型卤水所需的CO2来源具有复杂性。中国天然碱矿找矿远景应聚焦内蒙古第四纪盐湖、白垩纪封闭凹陷,以及河南始新世白云岩、油页岩发育的封闭盆地等地。积极开展天然碱矿成矿理论研究和找矿工作具有重要意义。
创新点:(1)系统综述了天然碱矿资源特征、成因要素;(2)根据天然碱矿资源分布特征和成因要素,预测了找矿远景。
Abstract:This paper is the result of mineral exploration engineering.
ObjectiveNatural soda-ash deposits, also called sodium carbonate evaporites, are mainly used to produce soda ash. In the soda industry, compared with synthetic soda, this process has advantages for environmental protection and lower cost. The downstream industries contain various fields, including glass, medicine, and so on. The booming of new energy industries such as photovoltaic has brought new demand for soda ash. China consumes loads of soda ash annually; thus, soda may become a scarce resource in the future. Mineralization process of natural soda-ash evaporites is of great significance in evaporite community. In China however, soda-ash deposits have received little attention and have great research potential. Investigations on types and distribution characteristics of global soda-ash deposits and their genesis are helpful and useful for deepening research and prospecting of Nacarbonate deposits in China.
MethodsCompiling published data and systematically summarizing the regional setting, provenance, genesis of typical soda-ash deposits.
ResultsMost deposits are Cenozoic in age and located in North America, Asia and Africa. Modern soda-ash deposits predominate in number, but resources size of ancient deposits is tremendous. Specific ancient deposits include the Green River Formation, United States and the Anpeng soda deposit in Biyang Depression in Henan, the counterparts, typical modern alkaline lake deposits formed in Lake Magadi, Kenya and Chaganor, Inner Mongolia.
ConclusionsThe requirement for Na-carbonate deposits precipitation are (1) hydrologically-closed basins, (2) arid climate, and (3) sufficient sodium carbonate supply. Unlike other evaporites, the sources of CO2 required to maintain Na-carbonate brines are complicated. Efforts for prospecting in China should focus on Quaternary saline lakes, Cretaceous restrict basins in Inner Mongolia, and tectonically-closed depression in which dolomite and oil shale develops in Henan. It is of great significance to carry out the metallogenic theory research and mine prospecting.
-
1. 引 言
世界卫生组织及中国饮用水标准规定砷浓度不可超过10 μg/L(WHO, 2017)。长期饮用高砷地下水可导致慢性砷中毒及皮肤癌等疾病,全球有70多个国家,超过1.5亿人的饮用水安全受到高砷地下水的威胁(韩双宝等,2010;郭华明等,2013;Wang et al., 2020;曹文庚等,2022; 张卓等,2023a)。沉积物中的固相砷是地下水中砷的主要来源。多数岩石中砷含量范围为0.5~2.5 μg/g(Mandal and Suzuki, 2002),松散沉积物中砷的含量范围通常为3~10 μg/g(Smedley and Kinniburgh, 2002; 何锦等,2020;马雪梅等,2020),富含砷矿物的沉积物中砷含量可达170 μg/g(Cook et al., 1995)。研究含水层中砷的迁移转化,除了查明沉积物总固态砷的含量,还需分析砷在沉积物中的赋存形态(van Herreweghe et al., 2003;朱丹尼等,2021;Drahota et al., 2021)。沉积物中固相砷赋存形态的微小差别可能引起地下水砷浓度的显著差异(Meharg et al., 2006; 张卓等,2023b)。分步提取实验是获取沉积物中砷赋存状态信息的主要手段。在之前的研究中,已经在分步提取过程中研究了萃取剂溶液的最优选择性(Paul et al., 2009;Eiche et al., 2010)。国外学者就河流三角洲沉积物中砷的赋存形态开展了大量研究。Eiche et al.(2008)研究表明,磷酸盐提取释放的强吸附砷是越南红河三角洲沉积物中砷的主要赋存形态。印度孟加拉三角洲平原的含水层中也发现了类似的结果(Neidhardt et al., 2014)。然而在内陆盆地,有关沉积物砷赋存形态的系统性研究相对缺乏。
河套盆地是中国西北地区典型的内陆盆地,地下水As浓度高达857 μg/L,远超中国饮用水标准(Guo et al., 2008)。因此,本研究选取河套盆地,通过刻画岩性与地球化学特征和开展砷的分步提取与解吸附实验,对比分析低砷和高砷含水层中沉积物砷的赋存形态与吸附特征。研究结果将有助于查明内陆盆地高砷地下水的形成机理,为合理开发可饮用地下水提供科学依据。
2. 研究区概况
河套盆地地处阴山隆起与鄂尔多斯台地之间,西界和北界均为狼山山前断裂,南界为鄂尔多斯北缘断裂,东界为乌梁素海断裂。研究区位于河套盆地西北侧,地处狼山山脉与主排干渠之间,包括山前冲洪积扇区和南部平原区,地理坐标为40°55′31″N~41°08′15″N,106°46′30″E~107°03′28″E(图1)。受沉积条件制约,研究区含水层具有明显的分带性。山前冲洪积扇区含水层沉积物主要由中砂、细砂组成,黏土在其中所占比重小于5%;平原区含水层沉积物主要由细砂、粉砂、粉质黏土和偶有泥炭夹层的淤泥质黏土组成,粉土和不同种类的黏土是其中的主要组成部分。
研究区浅层地下水受到大气降雨入渗补给、灌溉水补给和渠水的侧渗补给,深层地下水受到山前裂隙水的侧向补给和浅层地下水的垂向入渗。浅层地下水的排泄途径是蒸发作用、人工抽取、流入排干沟和垂向入渗到深层地下水,深层地下水的排泄路径是农业开采。原来研究区地下水流向大体是由西北向东南,但过度开采导致地下水流向逐渐转变为山前冲洪积扇由北向南、平原区由南向北的流动方向。地下水水化学类型受地势地貌、气候条件影响明显,具有显著的差异性。浅层地下水受强烈蒸发运移影响,水化学类型有HCO3−(Cl)−Na、Cl−HCO3−Na·Mg和Cl−SO4(HCO3)−Na·Mg型。深层地下水由山前冲洪积扇的Cl−HCO3−Ca·Mg型转变为平原区的Cl−Na型。高砷地下水主要分布在平原区(Zhang et al., 2020)。
3. 材料与方法
3.1 沉积物样品采集与测试
本研究从钻孔K02和K01中分别取出25和26个沉积物样品(图1)。其中,K2钻孔位于山前冲洪积扇区,坐标为41°01′07.37″N、106°57′41.41″E,钻孔深度约为80 m;K1钻孔位于平原区,位置坐标为41°00′13.73″N、106°58′16.85″E,钻孔深度约为81 m。获取的沉积物去掉外层沉积物后,马上用锡箔纸包裹,密封在装有纯N2(> 99.999%)的无菌塑料袋中,尽可能减少与O2的接触,并在−20℃的条件下保存。带回到实验室后,样品分装为两份,一份储存于−20℃的冰箱中,另一份进行冷冻干燥。
在色度分析和含水率测试之前,−20 ℃条件下保存的样品放入厌氧箱解冻。色度分析采用光谱色度计(CM-700d,Konica Minolta),测试之前对光谱色度计进行白板校正和零点校正。测试过程中保证切面平整,并在切口表面铺上一层高净度聚乙烯薄膜,每个样品测试3次。测试结束后计算出530 nm和520 nm的光谱反射差(R530-520),该差值能够指示沉积物的氧化还原环境(Horneman et al., 2004)。含水率测试采用通用的烘干法,用铝盒准确称取烘干前的原状土样质量,放入105℃恒温干燥箱中烘干后放入干燥器冷却,准确称量烘干后的土样质量,通过计算得出含水率。
沉积物电导率和pH的测量采用Bélanger and VanRees(2007)的方法。冷冻干燥后的沉积物与去离子水以1∶5的比例置于PE离心管中,25℃状态下以150 rpm转速震荡1 h。震荡完毕后,将离心管置于离心机中以5000 rpm转速离心20 min并取上清液用0.22 μm纤维滤膜过滤。所得部分滤液通过电导率仪(DDS-307A, SHKY)进行电导率的检测,所得电导率值可以反映出沉积物的可溶性组分含量。沉积物样品与超纯水以1∶2.5比例充分混合后,摇匀,静置1 h使用pH检测仪(HI 8424,HANNA)对其进行pH测定。
沉积物样品中的主量和微量元素的测定采用手持便携式XRF仪(XL3t800, Thermo Niton)进行测定,测试元素主要包括Ca、Sr、As、Fe和Mn。测试之前将样品冷冻干燥,并研磨至200目,取适量于专用测量杯中,压实后放置在手持XRF仪光源处,每个样品测试3次。2个标准物质(GBW07303,GBW07305)用于确保数据的准确性,测试偏差均小于20%,其中As元素的测试偏差均小于5%。
3.2 室内实验
3.2.1 分步提取实验
为查明沉积物中砷的赋存状态,本研究开展了分步提取实验(Sequential extraction procedure,SEP)。分步提取方法参照Eiche et al.(2008, 2010)的研究,该提取方法也是基于Keon et al.(2001)和Wenzel et al.(2001)等研究的改进(表1)。每个新鲜沉积物样称取0.5 g,放入离心管中,加入适量的提取剂。由于分步提取后提取液盐度较高,需稀释测试,这就要求测试仪器需要较低的检出限和较高的分析精度。ICP−MS的分析精度为±3.0%,检出限为0.01 μg/L,能够满足测试要求。其中分步提取第六步(F6)的提取液中含有高浓度的HF,会损坏仪器影响测试精度。因此,F6的提取液在测试之前,需要在电热板加热进行赶酸处理。
表 1 分步提取实验具体步骤Table 1. Sequential extraction procedure步骤 目标物 提取剂 条件 F1 弱吸附态砷 0.05 mol/L (NH4)2SO4 25 mL,25℃,4 h,重复一次,水洗一次 F2 强吸附态砷 0.5 mol/L NaH2PO4 40 mL,25℃,16 h及24 h各一次,每个时间段重复一次,水洗一次 F3 与可挥发硫化物、碳酸盐、锰氧化物和完全无定形态的铁氧化物或氢氧化物共存的砷 1 mol/L HCl 40 mL,25℃,1 h,重复一次,水洗一次 F4 与无定形态铁氧化物或氢氧化物共存的砷 0.2 mol/L NH4H2C2O3 40 mL,25℃,2 h,pH=3,黑暗条件下,重复一次,水洗一次 F5 与结晶态铁氧化物或氢氧化物
共存的砷0.5 mol/L NaC6H8O7
1 mol/L NaHCO3,Na2S2O4XH2O35 mL NaC6H8O7+2.5 mL NaHCO3(加热至85℃),加0.5 g Na2S2O4XH2O,15 min在85℃,重复一次,水洗一次 F6 与硅酸盐有关的砷 10 mol/L HF,H3BO3 40 mL,25℃,1 h、24 h、16 h后各加5 g硼酸,每个时间段重复一次,热水洗一次 F7 含砷硫化物,与硫化物和有机质
共沉淀的砷16 mol/L HNO3,30% H2O2 先加入10 mL HNO3,反应过后加入多次30%过氧化氢,加热,冷却后稀释到100 mL,离心、过滤、测试 3.2.2 解吸附实验
本研究从钻孔K02和K01各选取一个典型沉积物进行pH和反离子效应对砷的解吸附影响的批实验。该实验主要包括三部分内容:解吸附动力学实验、pH对解吸附影响的实验、反离子效应(Na/Ca0.5(M/M))对砷解吸附影响的实验。
(1)解吸附动力学实验
为查明砷解吸附达到平衡的时间,本研究开展了解吸附动力学实验。分别称取0.6 g新鲜沉积物放入厌氧瓶中,然后加入24 mL、125 mmol/L NaCl和1.5 mmol/L CaCl2的混合溶液,用橡胶塞封闭,整个过程在厌氧箱中操作,设置3个平行样。混合溶液离子强度约为130 mmol/L,Na/Ca0.5比值约为102,pH值为7.6。为保证沉积物颗粒与溶液均匀混合,超声15 min后放入150 r/min的恒温振荡箱中。取样间隔为1 h、3 h、5 h、7 h、10 h、14 h、20 h、28 h、36 h、48 h和60 h。取样之前保证溶液混合均匀,每次取样量为2 mL,用0.22 μm过滤器过滤到2 mL离心管中,放入4℃冰箱中保存,一周之内完成测试工作。
(2)pH对解吸附影响的实验
控制Na/Ca0.5(M/M)比值约为102和离子强度约为130 mmol/L,探究不同pH值对沉积物中砷解吸附的影响。将Na/Ca0.5比值为102的NaCl和CaCl2的混合溶液分装为5份,并将溶液pH值分别调到5.4、6.7、7.6、8.6和9.6。在5个厌氧瓶中,分别称取0.6g新鲜沉积物,并加入24 mL不同pH值梯度的NaCl和CaCl2的混合溶液,用橡胶塞封闭,整个过程在厌氧箱中操作,设置3个平行样。所有加入沉积物和混合溶液的厌氧瓶,超声15 min后放入150 r/min的恒温振荡箱中。60 h后取样,用0.22 μm过滤器过滤到离心管中,放入4℃冰箱中保存,一周之内完成测试工作。
(3)反离子效应对砷解吸附影响的实验
控制离子强度为(130±5)mmol/L,通过改变NaCl和CaCl2的浓度来改变Na/Ca0.5比值(表2)。在7个厌氧瓶中,分别称取0.6 g新鲜沉积物,并分别加入24 mL不同Na/Ca0.5比值梯度的NaCl和CaCl2的混合溶液,用橡胶塞封闭,整个过程在厌氧箱中操作,设置3个平行样。所有加入沉积物和混合溶液的厌氧瓶,超声15 min后放入150 r/min的恒温振荡箱中。60 h后取样,用0.22 μm过滤器过滤到离心管中,放入4℃冰箱中保存,一周之内完成测试工作。
表 2 离子强度为(130±5)mmol/L条件下,不同浓度NaCl和CaCl2混合液的Na/Ca0.5(M/M)比值Table 2. Na/Ca0.5(M/M) ratio of the mixed solution of different concentrations of NaCl and CaCl2 under the condition of ionic strength of about (130±5) mmol/LNaCl/(mmol/L) CaCl2/(mmol/L) Na/Ca0.5 2 43 0.3 5 42 0.7 10 40 1.6 30 35 5.0 60 23 13 110 7 42 125 1.5 102 4. 结果与讨论
4.1 沉积物的岩性特征
研究区的山前冲洪积扇区钻孔K02和平原区钻孔K01沉积物的岩性特征如图2所示。钻孔K02沉积物的组成是从粗砂到黏土,而钻孔K01主要从中砂到黏土。对于钻孔K02,14 m以上的沉积物主要由砂质黏土和粉质黏土组成,14~42 m主要以砂质含水层为主。在42~44 m存在约2 m厚的黏土层,42 m以下主要以砂质含水层为主同时伴有砂质黏土互层(图2a)。与钻孔K02不同,位于平原区的钻孔K01沉积物颗粒整体较细且含有大量的黏土互层。其中,8 m以上主要以黏土为主,8~40 m则主要以砂质含水层为主并且常常伴有砂质黏土互层,40~42 m出现黏土层,42 m以下为颗粒较细的细砂含水层,这个研究结果与Shen et al.(2018)一致。总体来看,研究区近表层沉积物主要以粉质黏土为主,地表以下10~40 m是砂质含水层,地表以下40 m处存在1~2 m厚的相对连续的黏土层将40 m以上和约42 m以下的含水层隔开。
沉积物的色度特征能够指示沉积物的氧化还原环境和铁氧化物的还原程度(Horneman et al., 2004)。钻孔K02和K01沉积物色度随深度的变化均是由浅黄色变为深灰色,说明深部含水层处于一个相对还原的环境当中,铁氧化物的还原程度也较强。而从整体来看,两个钻孔的色度特征有较大差异,相对于钻孔K02,钻孔K01的沉积物色度更深,这可能是因为平原沉积物颗粒较细,含水层处于更封闭的还原环境,铁锰氧化物的还原程度更强(van Geen et al., 2013)。
沉积物含水率主要受其岩性控制。两个钻孔表层5 m以上沉积物尽管颗粒较细,含水率仍然较低,主要由于其处于非饱和带。而在饱和带,沉积物含水率随深度的变化主要受岩性影响,沉积物岩性颗粒越细,含水率越高。两个钻孔沉积物电导率在近地表较高(图2),主要是因为研究区为干旱半干旱气候,蒸发蒸腾作用较强,使得近地表沉积物含有大量的可溶盐(Yuan et al., 2017)。沿深度随沉积物岩性的变化而波动,沉积物岩性越细,电导率越大,这是由于颗粒较细的黏土颗粒表面有大量可交换的离子。此外,由于钻孔K01位于平原区,沉积物颗粒整体较细且地下水水位埋深较浅蒸发作用强,导致其沉积物电导率(均值为395 μS/cm)大于钻孔K02(均值为308 μS/cm)。
4.2 沉积物的地球化学特征
研究区沉积物中0~10 m、40~45 m和75~80 m含水层位的Ca和Sr的含量明显高于其他含水层(图3)。微量元素As、Fe和Mn也有相似的分布特征。沉积物的岩性特征表明,10 m以上的沉积物主要以黏土和粉质黏土为主,40~45 m是不连续的黏土层,而75~80 m也是颗粒较细的黏土层。对比钻孔的黏土层和砂层沉积物的地球化学特征发现,K02钻孔黏土层沉积物Ca含量中值为53.6 mg/g,而砂层沉积物Ca含量中值为33.0 mg/g;K01钻孔中两者中值分别为48.3 mg/g和31.6 mg/g。黏土层和砂层沉积物中微量元素的含量差异更为明显,K02钻孔黏土层沉积物As含量中值为17.6 μg/g,而砂层沉积物As含量中值为8.6 μg/g;K01钻孔中两者中值分别为20.1 μg/g和7.9 μg/g。这主要是因为砂层沉积物中富含石英,含Ca和Sr矿物的含量低于黏土层(李晓峰,2018)。其次是因为黏土层表面吸附能力强,能够吸附As、Fe和Mn等微量元素(崔邢涛等,2015)。
两个钻孔沉积物的地球化学特征也有一定的差异。普遍表现为钻孔K02的Ca、Sr、As、Fe和Mn含量大于钻孔K01,且在深层沉积物中表现更为明显(图3)。钻孔K02沉积物中Ca的含量范围为12.2~86.9 mg/g,平均值为37.9 mg/g,钻孔K01沉积物中Ca的含量范围为9.6~68.7 mg/g,平均值为35.7 mg/g。K02钻孔沉积物中As的浓度范围为4.6~33.1 μg/g,平均值13.1 μg/g;K01钻孔沉积物中As的浓度范围为5.3~34.0 μg/g,平均值12.9 μg/g,表明冲洪积扇边缘地区沉积物总As的含量略大于平原区。两个钻孔沉积物中Fe和Mn含量的差异更为明显,钻孔K02沉积物中Fe的含量比K01高13.7%,其Mn的含量比K01高14.1%。这主要是由于钻孔K01位于平原区,沉积物经历了更强的风化作用,且积物颗粒整体较细,地下水流速慢,水岩作用强烈,有利于沉积物中化学组分向地下水中释放(张文凯等,2020)。此外,平原区含水层较为封闭,沉积物的色度特征也表明含水层长期处于较为还原的环境中,变价微量元素被还原为较低价态,易于向地下水中迁移。因此,钻孔K02和K01沉积物地球化学的微小差异主要受沉积环境和水动力条件控制。
4.3 沉积物中砷的赋存形态
山前冲洪积扇的含水层的沉积物岩性主要以中砂、细砂和黏土为主,平原区含水层的沉积物则以细砂、粉砂和黏土为主。因此,本研究从钻孔K02和K01各选取3个不同岩性的代表性沉积物用于分步提取实验(SEP)(表3)。实验过程选用GBW07303和GBW07305作为标准样品检验回收率,结果表明:对于GBW07303不同状态As的提取实验的回收率分别为81%,GBW07305不同状态As的提取实验的回收率分别为88%。分步提取实验获取的7种形态砷的总和与XRF测得的总固相砷的相对偏差均小于10%。
表 3 用于分步提取的沉积物信息Table 3. Sediment information for SEP编号 岩性 采样深度/m K02−M 中砂 38.35 K02−F 细砂 62.25 K02−C 黏土 41.95 K01−F 细砂 55.15 K01−S 粉砂 30.95 K01−C 黏土 37.85 分步提取结果表明,K02钻孔中砂、细砂和黏土沉积物固相砷主要以与可挥发硫化物、碳酸盐、锰氧化物和完全无定形态的铁氧化物或氢氧化物共存的砷(F4)为主,占比分别为33%、40%和43%(图4a、b、c)。其次是结晶态铁氧化物或氢氧化物结合态(F5)和强吸附态砷(F2)。砂层沉积物中与无定形态铁氧化物或氢氧化物结合的固相砷(F3)占比大于与硅酸盐结合的砷(F6),前者占比均大于10%,后者均小于5%,而黏土沉积物中两者的占比分别为7%和12%。最容易释放到地下水中的弱吸附态砷(F1)和最顽固的与硫化物和有机质共沉淀的固相砷(F7)占比较小,均低于5%。钻孔K01细砂沉积物的固相砷以F4为主(35%),其次分别是F2(32%)和F6(16%)(图4d)。粉砂和黏土沉积物则以F2为主(分别为43%和40%),其次以F4为主(分别为12%和18%);两个沉积物中F3所占的比例均超过10%(图4e、f)。细砂、粉砂和黏土沉积物中F1和F7均小于5%。
对比山前冲洪积扇的钻孔K02和平原区的钻孔K01发现,前者沉积物中固相砷主要以F4为主,后者则主要以F2为主。钻孔K02黏土沉积物中F4达到11.3 μg/g,明显高于K01的4.6 μg/g。而钻孔K02黏土沉积物中F2仅有5.8 μg/g,低于钻孔K01的10.3 μg/g(图4c、f)。钻孔K01砂层沉积物中的F2也明显大于K02。此外,平原区沉积物的F3含量也大于山前冲洪积扇沉积物。这说明平原区沉积物经历更强的风化侵蚀作用后,固相砷活性增强,向更具迁移性的吸附态和完全无定形铁氧化物或氢氧化物结合态转化。大量研究表明吸附态的砷迁移性较强,通过竞争解吸附或者弱碱条件下的解吸附,更容易释放到地下水中,而无定形态铁氧化物或氢氧化物结合态砷相对稳定,需要通过还原性溶解才能释放到地下水中(Smedley and Kinniburgh, 2002)。这也解释了为何平原区地下水砷浓度普遍高于山前冲洪扇的地下水(李晓峰,2018; Zhang et al., 2020)。除了含水层沉积物本身物源的影响,含水层所处的环境和地下水的化学特征也会影响砷的解吸附。
4.4 砷的解吸附
以往的研究表明,研究区地下水pH和Na/Ca0.5(M/M)与砷浓度均有较好的正相关关系(Zhang et al., 2020),因此,本研究选取钻孔K02和K01的沉积物(表3),分别探讨了pH和Na/Ca0.5(M/M)对砷解吸附的影响。动力学实验结果表明,在pH为7.6、离子强度为130 mmol/L和Na/Ca0.5比值为102的条件下,砷解吸附能够48 h时基本达到平衡(图5a)。为确定砷解吸附达到平衡,实验设定反应时间为60 h。
4.4.1 pH的影响
实验设定离子强度为130 mmol/L,Na/Ca0.5比值为102。pH条件分别设定为5.4、6.7、7.6、8.6和9.6。当pH为5.4时,K02−F和K01−F沉积物释放的砷占总吸附砷的比值分别为0.54和0.44;当pH升高至6.7时,砷释放量所占总吸附砷比值分别降为0.32和0.30(图5b),这可能是因为较低的pH可能使铁氧化物发生少量溶解导致砷的释放。pH从6.7上升至8.6的过程中,沉积物砷的释放量并没有明显增加,仅上升0.03左右。而pH由8.6上升至9.6,沉积物砷的释放量显著增加,释放量上升0.15。这是由于随着pH升高沉积物颗粒表面带负电荷,与含砷阴离子形成静电斥力导致吸附态的砷发生解吸附,进入水溶液中(Masue et al., 2007)。
4.4.2 反离子效应的影响
许多学者认为,沉积物颗粒表面存在扩散双电子层(Dzombak and Morel, 1990; 刘新敏,2014),相比于以Na+为主的地下水系统,以Ca2+为主的地下水系统能够导致带负电的沉积物颗粒表面与带负电的含砷弱阴离子之间的斥力减小,有利于砷的吸附,这种现象被称为反离子效应(Masue et al., 2007; Fakhreddine et al., 2015)。当水中离子强度一定时,带有两个正电荷Ca2+被单电荷Na+替换时,即Na/Ca0.5比值增加时,这种反离子效应就会减弱,促进吸附态的砷释放到地下水中。
实验过程中保持pH和离子强度不变,通过调节溶液中Na+和Ca2+浓度改变Na/Ca0.5比值。结果表明,砷解吸附的量随Na/Ca0.5比值的增加而增加(图5c)。当Na/Ca0.5比值为0.3时,K02−F和K01−F沉积物砷的解吸附量占总吸附态砷的比值分别为0.12和0.11。而当Na/Ca0.5比值增加到102时,K02−F沉积物砷的解吸附量占总吸附态砷的比值能够达到0.37,在K01−F沉积物中这一比值为0.47。
4.5 地下水开发利用建议
河套盆地是中国的塞上粮仓,对水资源的需求较大。研究区地势较高,引黄河入河套盆地并难以满足居民的农业和生活需求,因此,居民普遍开采地下水用于农业灌溉和日常生活,这虽然解决水量的问题,却忽视了原生劣质地下水的危害。根据国家《生活饮用水卫生标准》(GB 5749—2022)和《地下水质量标准》(GB/T14848—2017),砷浓度大于10 μg/L的地下水为高砷地下水,摄入后对人体有害。以往的研究发现高砷地下水主要集中在平原区,浓度高达857 μg/L(Guo et al., 2008)。本研究发现,山前冲洪积扇区的含水层沉积物固相砷相对稳定,而平原区的含水层沉积物固相砷迁移性相对较强,且平原区沉积物吸附态砷在弱碱性和高Na/Ca0.5摩尔比值条件下,容易向地下水迁移,导致砷的富集。因此,当地居民种植农作物时避免使用碱性复合肥,从而减少碱性水的向下补给。此外,生活污水中Na+较高,建议适当处理后排放。用于日常生活的地下水,建议采用混凝沉淀或吸附法降砷。
5. 结 论
山前冲洪积扇区含水层处于相对氧化的环境中,其沉积物以细砂和中粗砂为主,而平原区含水层处于封闭的还原环境中,沉积物以粉细砂为主。两者沉积物总固相砷含量相差不大,但固相砷的赋存形态差别较大。山前冲洪积扇区含水层沉积物固相砷以与可挥发硫化物、碳酸盐、锰氧化物和完全无定形态的铁氧化物或氢氧化物共存的砷为主(33%~43%),平原区含水层沉积物固相砷则以强吸附态砷为主(32%~43%),后者沉积物的中固相砷迁移性更强,容易通过解吸附释放到地下水中。此外,当pH值由6.1上升到9.6时,山前和平原区沉积物解吸附砷占总吸附砷的比值分别上升0.16和0.22。同时,Na/Ca0.5摩尔比值的增加,会导致反离子效应减弱,比值由0.3增加到102时,山前沉积物和平原区解吸附砷占总吸附砷的比值分别上升0.26和0.36。可见含水层中pH的升高和Na/Ca0.5摩尔比值的增加,都会促使沉积物中的砷发生解吸附,导致地下水中砷的富集。因此,当地居民应减少碱性以及富含Na+的生产生活用水的排放,同时平原区用于日常生活的地下水,建议当地居民采用混凝沉淀或吸附法降砷。
-
图 1 世界主要天然碱矿分布图(据Warren, 2010)
Figure 1. Distribution of major soda-ash deposit in the world (after Warren, 2010)
图 2 形成天然碱矿的不饱和卤水演化路径(改自Hardie and Eugster, 1970)
Figure 2. Brine evolution pathways with different brine types for Na-carbonate evaporites formation (modified from Hardie and Eugster, 1970)
图 3 苏打石-天然碱-泡碱随温度、CO2的变化的稳定图(改自Jagniecki et al., 2015)
Figure 3. Stability fields of nahcolite, natron, and trona as function of temperature and pCO2 (modified from Jagniecki et al., 2015)
图 5 贝帕扎里盆地地质简图及天然碱矿床地层柱状简图(据Helvaci, 2010)
Figure 5. Geological profile of the Beypazari Basin and schematic stratigraphic column of soda deposits (after Helvaci, 2010)
图 6 泌阳凹陷地质简图及安棚天然碱矿床岩性柱状简图(据杨江海等, 2014; 易承龙, 2016)
Figure 6. Geological map of the Biyang depression and schematic stratigraphic column of soda deposits (after Yang Jianghai et al., 2014; Yi Chenglong, 2016)
图 7 东非大裂谷东支区域图(a)与马加迪湖区地质图(b)(改自Eugster and Blair, 1968; Schubel and Simonson, 1990; Damnati et al., 1992; Owen et al., 2019)
Figure 7. The eastern branch of the East African Rift Valley (a) and geological map of Lake Magadi (b) (modified from Eugster and Blair, 1968; Schubel and Simonson, 1990; Damnati et al., 1992; Owen et al., 2019)
图 8 查干诺尔区域地质简图和泡碱矿床岩性柱状简图(改自张晨鼎, 2013)
Figure 8. Schematic map of Chaganor lake area and stratigraphic column of soda deposits (modified from Zhang Chending, 2013)
表 1 纯碱矿物产量及碳酸钠储量(据USGS, 2023)
Table 1 Soda-ash mine production and reserves (after USGS, 2023)
表 2 世界典型天然碱矿床(据张晨鼎, 2013; 叶铁林, 2013; Warren, 2016)
Table 2 Typical Na-carbonate deposits in the world (after Zhang Chending, 2013; Ye Tielin, 2013; Warren, 2016)
-
Baker B H, Mitchell J G. 1976. Volcanic stratigraphy and geochronology of the Kedong-Olorgesailie area and the evolution of the South Kenya rift valley[J]. Journal of the Geological Society, 132: 467-484. doi: 10.1144/gsjgs.132.5.0467
Baker B H. 1986. Tectonics and volcanism of the southern Kenya Rift Valley and its influence on rift sedimentation[J]. Geological Society London Special Publications, 25(1): 45-57. doi: 10.1144/GSL.SP.1986.025.01.05
Chen Jianli. 2013. Geological characteristics, genesis and ore prediction of natural soda deposit in Biyang depression: Taking Anpeng trona deposit as an example[J]. Contributions to Geology and Mineral Resources Research, 28(3): 393-400 (in Chinese with English abstract).
Crossley R. 1979. The Cenozoic stratigraphy and structure of the western part of the Rift Valley in southern Kenya[J]. Journal of the Geological Society, 136(4): 393-405. doi: 10.1144/gsjgs.136.4.0393
Culbertson W C. 1972. Trona and Halite resources in Wilkins Peak Member of Green river Formation, Green River Basin, Wyoming: Abstract[J]. AAPG Bulletin, 56: 612-612.
Cumming V M, Selby D, Lillis P G. 2012. Re-Os geochronology of the lacustrine Green River Formation: Insights into direct depositional dating of lacustrine successions, Re-Os systematics and paleocontinental weathering[J]. Earth and Planetary Science Letters, 359/360: 194-205. doi: 10.1016/j.epsl.2012.10.012
Cumming V M, Selby D, Lillis P G, Lewan M D. 2014. Re-Os geochronology and Os isotope fingerprinting of petroleum sourced from a Type Ⅰ lacustrine kerogen: Insights from the natural Green River petroleum system in the Uinta Basin and hydrous pyrolysis experiments[J]. Geochimica et Cosmochimica Acta, 138: 32-56. doi: 10.1016/j.gca.2014.04.016
Damnati B, Taieb M, Williamson D. 1992. Laminated deposits from Lake Magadi (Kenya); climatic contrast effect during the maximum wet period between 12, 000-10, 000 yrs BP[J]. Bulletin De La Société Géologique De France, 163(4): 407-414.
Demicco R V, Lowenstein T K, Hardie L A. 2003. Atmospheric pCO2 since 60 Ma from records of seawater pH, calcium, and primary carbonate mineralogy[J]. Geology, 31(9): 793-796. doi: 10.1130/G19727.1
Dickinson W R, Klut M A, Hayes M J, Janecke S U, Lundin E R, McKittrick M A, Olivares M D. 1988. Paleogeographic and paleotectonic setting of Laramide sedimentary basins in the central Rocky Mountain region[J]. Geological Society of America Bulletin, 100: 1023-1039. doi: 10.1130/0016-7606(1988)100<1023:PAPSOL>2.3.CO;2
Earman S, Phillips F M, McPherson B J O L. 2005. The role of "excess" CO2 in the formation of trona deposits[J]. Applied Geochemistry, 20(12): 2217-2232. doi: 10.1016/j.apgeochem.2005.08.007
Eugster H P. 1967. Hydrous sodium silicates from lake Magadi, Kenya: Precursors of bedded chert[J]. Science, 157(3793): 1177-1180. doi: 10.1126/science.157.3793.1177
Eugster H P, Blair F J. 1968. Gels composed of sodium-aluminum silicate, Lake Magadi, Kenya[J]. Science, 161(3837): 160-163. doi: 10.1126/science.161.3837.160
Eugster H P. 1969. Inorganic bedded cherts from the Magadi area, Kenya[J]. Contributions to Mineralogy and Petrology, 22(1): 1-31. doi: 10.1007/BF00388011
Eugster H P, Hardie L A. 1975. Sedimentation in an ancient PlayaLake Complex: The Wilkins peak member of the Green River Formation of Wyoming[J]. Geological Society of America Bulletin, 86(3): 319-334. doi: 10.1130/0016-7606(1975)86<319:SIAAPC>2.0.CO;2
Eugster H P, Jones B F, Shirley L R. 1977. Hydrochemistry of the Lake Magadi basin, Kenya[J]. Geochimica et Cosmochimica Acta, 41(1): 53-72. doi: 10.1016/0016-7037(77)90186-7
Eugster H P. 1980. Hypersaline Brines and Evaporitic Environments[M]. Netherland: Elsevier Science.
Garcia-Veigas J, Gundogan I, Helvacı C, Prats E. 2013. A genetic model for Na-carbonate mineral precipitation in the Miocene Beypazari trona deposit, Ankara province, Turkey[J]. Sedimentary Geology, 294: 315-327. doi: 10.1016/j.sedgeo.2013.06.011
Hardie L A, Eugster H P. 1970. The evolution of closed-basin brines[J]. Mineralogical Society of America Special Publication, 3: 253-273.
Helvaci C, Inci U, Yilmaz H, Yaǧmurlu F. 1989. Geology and Neogene trona deposit of the Beypazari region, Turkey[J]. Turkish Journal of Engineering and Environmental Sciences, 13: 245-256.
Helvaci C. 1998. The Beypazari Trona Deposit, Ankara Province, Turkey[C]//Wyoming State Geological Survey Public Information Circular 40.
Helvaci C. 2010. Geology of the Beypazari trona field, Ankara, Turkey. Mid-congress field excursion guide book, tectonic crossroads: Evolving orogens of Eurasia-Africa-Arabia[C]//Tectonic Crossroads: Evolving Orogens of Eurasia-AfricaArabia, Ankara, Turkey.
Hyland E G, Sheldon N. 2013. Coupled CO2-climate response during the Early Eocene climatic optimum[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 369: 125-135. doi: 10.1016/j.palaeo.2012.10.011
Inci U. 1991. Miocene alluvial fan-alkaline playa lignite-trona bearing deposits from an inverted basin in Anatolia: Sedimentology and tectonic controls on deposition[J]. Sedimentary Geology, 71(1): 73-97.
Jagniecki E A, Lowenstein T K, Jenkins D M, Demicco R V. 2015. Eocene atmospheric CO2 from the nahcolite proxy[J]. Geology, 43(12): 1075-1078.
Jagniecki E A, Lowenstein T K, Demicco R V, Baddouh M, Carroll A R, Beard B L, Johnson C M. 2021. Spring origin of Eocene carbonate mounds in the Green River Formation, Northern Bridger Basin, Wyoming, USA[J]. Sedimentology, 68(6): 2334-2364. doi: 10.1111/sed.12852
Jin Qiang, Xiong Shousheng, Lu Peide. 1998. Volcanic activity in major source rocks in faulted basins of China and its significance in main source rocks of fault basins in China [J]. Geological Review, 44: 136-142 (in Chinese with English abstract). doi: 10.3321/j.issn:0371-5736.1998.02.004
Lee H, Muirhead J. D, Fischer T P, Ebinger C J, Sharp Z. D, Kianji G. 2016. Massive and prolonged deep carbon emissions associated with continental rifting[J]. Nature Geoscience, 9: 145-149. doi: 10.1038/ngeo2622
LeGall B, Nonnotte, P, Rolet J, Benoit M, Guillou, H, MousseauNonnotte M, Albaric J, Deverchère J. 2008. Rift propagation at craton margin. Distribution of faulting and volcanism in the North Tanzanian Divergence (East Africa) during Neogene times[J]. Tectonophysics, 448, 1-19. doi: 10.1016/j.tecto.2007.11.005
Lowenstein T K, Demicco R V. 2006. Elevated Eocene atmospheric CO2 and its subsequent decline[J]. Science, 313(5795): 1928-1928. doi: 10.1126/science.1129555
Lowenstein T K, Lauren A D, García-Veigas, J. 2016. Influence of magmatic-hydrothermal activity on brine evolution in closed basins: Searles Lake, California[J]. Geological Society of America Bulletin, 128: 1555-1568. doi: 10.1130/B31398.1
Lowenstein T K, Jagniecki E A, Carroll A R, Smith M E, Renaut R W, Owen R B. 2017. The Green River salt mystery: What was the source of the hyperalkaline lake waters?[J]. Earth-Science Reviews, 173: 295-306. doi: 10.1016/j.earscirev.2017.07.014
Lowenstein T K, Demicco R V. 2019. When evaporites are not formed by evaporation-The role of temperature and pCO2 on saline deposits of the Eocene Green River Formation, Colorado, USA[J]. Geological Society of America Bulletin, 132: 1365-1380.
McNulty E. 2017. Lake Magadi and the Soda Lake Cycle: A Study of the Modern Sodium Carbonates and of Late Pleistocene and Holocene Lacustrine Core Sediments[D]. Binghamton: Binghamton University.
Muirhead J D, Simon A K, Hyunwoo L, Sara M, Brent D T, Tobias P F, Kianji G W, Sarah S D. 2016. Evolution of upper crustal faulting assisted by magmatic volatile release during early-stage continental rift development in the East African Rift[J]. Geosphere, 12: 1670-1700. doi: 10.1130/GES01375.1
Ogola J S, Behr H J. 2000. Mineralogy and Trona Formation in lake Magadi, Kenya[C]//Applied Mineralogy in Research, Economy, Technology, Ecology and Culture: Proceedings of the Sixth International Congress on Applied Mineralogy. Göttingen, GER: ICAM, 383-386.
Olson K J, Lowenstein T K. 2021. Searles Lake evaporite sequences Indicators of late Pleistocene-Holocene lake temperatures, brine evolution, and pCO2[J]. Geological Society of America Bulletin, 133(11): 2319-2334.
Owen R B, Muiruri V M, Lowenstein T K, Renaut R W, Rabideaux N, Luo S, Deino A L, Sier M J, Dupont-Nivet G, McNulty E P, Leet K, Cohen A, Campisano C, Deocampo D, Shen C, Billingsley A, Mbuthia A. 2018. Progressive aridification in East Africa over the last half million years and implications for human evolution[J]. Proceedings of the National Academy of Sciences, 115(44): 11174-11179. doi: 10.1073/pnas.1801357115
Owen R B, Renaut R W, Muiruri V M, Rabideaux N M, Lowenstein T K, McNulty E P, Leet K, Deocampo, D, Luo S, Deino A L, Cohen, A, Sier M J, Campisano, C, Shen C, Billingsley A, Mbuthia A, Stockhecke M. 2019. Quaternary history of the Lake Magadi Basin, southern Kenya Rift: Tectonic and climatic controls[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 518: 97-118. doi: 10.1016/j.palaeo.2019.01.017
Qi Bingde, Xie Xing, Deng Liangwu, Zeng Qingliang, Ma Lanlan. 2011. Analysis on metallogenic conditions of Chaganur natural alkali deposit in Sunite Right Banner, Inner Mongolia[J]. Journal of Sichuan Geology, 31(S1): 8-10 (in Chinese).
Renaut R W, Ashley G M. 2002. Sedimentation in Continental Rifts[M]. US: SEPM Society for Sedimentary Geology.
Renaut R W, Owen R B, Lowenstein T K, De C G, Mcnulty E, Scott J J, Mbuthia A. 2021. The role of hydrothermal fluids in sedimentation in saline alkaline lakes: Evidence from Nasikie Engida, Kenya Rift Valley[J]. Sedimentology, 68(1): 108-134. doi: 10.1111/sed.12778
Schubel K A, Simonson B M. 1990. Petrography and diagenesis of cherts from Lake Magadi, Kenya[J]. Journal of Sedimentary Petrology, 60(5): 761-776.
Smith M E, Carroll A R, Scott J J, Singer B S. 2014. Early Eocene carbon isotope excursions and landscape destabilization at eccentricity minima: Green River Formation of Wyoming[J]. Earth and Planetary Science Letters, 403: 393-406. doi: 10.1016/j.epsl.2014.06.024
USGS. 2023. Mineral Commodity Summaries[R]. United States Geological Survey.
Wang Aiyun, Chen Wenxi. 2022. Trona: From ancient washing powder to the mother of chemical industry[J]. Earth, (1): 6-11 (in Chinese).
Wang Jiuyi, Liu Chenglin, Wang Chunlian, Yu Xiaocan, Yan Kai, Gao Chao. 2021. Tectono-paleoclimatic coupling process for mineralization of Late Cretaceous-Paleogene evaporites in South China[J]. Acta Geologica Sinica, 95(7): 2041-2051 (in Chinese with English abstract).
Warren J K. 2010. Evaporites through time: Tectonic, climatic and eustatic controls in marine and nonmarine deposits[J]. EarthScience Reviews, 98(3/4): 217-268.
Warren J K. 2016. Evaporites: Sediments, Resources and Hydrocarbons[M]. Switzerland: Springer International Publishing AG.
Xu Hong, Peng Qiming, Martin R Palmer. 2004. Origin of tourmalinerich rocks in a Paleoproterozoic terrene (N.E. China): Evidence for evaporite-derived boron[J]. Geology in China, (3): 240-253 (in Chinese with English abstract).
Xu Yang, Cao Yangtong, Liu Chenglin. 2021. Whether the middle Eocene salt-forming brine in the Kuqa Basin reached the potashforming stage: Quantitative evidence from halite fluid inclusions[J]. Geofluids, 2: 1-12.
Xu Yang, Liu Chenglin, Cao Yangtong. 2021. Salt-forming evolution characteristics of Middle Eocene in the Kuqa basin, Xinjiang: A case study of borehole KL4[J]. Acta Geologica Sinica, 95(7): 2183-2192 (in Chinese with English abstract).
Yan Kai, Liu Chenglin, Wang Chunlian, Fan Meiling, Xu Haiming, Wang Jiuyi. 2021. Mineral deposition and paleoenvironment of Cretaceous evaporite in Southwestern Congo[J]. Acta Petrologica et Mineralogica, 40(3): 525-534 (in Chinese with English abstract).
Yang Jianghai, Yi Chenglong, Du Yuansheng, Zhang Zongheng, Yan Jiaxin. 2014. Geochemical characteristics of Paleogene alkalibearing rock series in Biyang Depression and their indicative significance for alkali-formation[J]. Scientia Sinica (Terrae), 44(10): 2172-2184 (in Chinese). doi: 10.1360/zd-2014-44-10-2172
Ye Tielin. 1978. Trona and its genesis briefly[J]. Industrial Minerals and Processing, 6: 18-23 (in Chinese).
Ye Tielin. 2013. Trona Resources, Geology, Mining and Processing, 3rd Edition[M]. Beijing: Chemical Industry Press (in Chinese).
Yi Chenglong. 2016. Sequence stratigraphy characteristics and its significance of alkaliferous strata of the Paleogene Hetaoyuan Formation in Anpeng area, Biyang sag in Henan Province[J]. Journal of Palaeogeography (Chinese Edition), 18(1): 93-100 (in Chinese with English abstract).
Zachos J C, Dickens G R, Zeebe R E. 2008. An Early Cenozoic perspective on greenhouse warming and carbon-cycle dynamics[J]. Nature, 451: 279-283. doi: 10.1038/nature06588
Zhang Chending. 1979. Soda-ash deposits[J]. Soda Industry, 3: 56-64(in Chinese).
Zhang Chending. 2004. Development of the Trona deposit in Beypazari, Turkey [J]. Soda ash Industry, (2): 15-18 (in Chinese).
Zhang Chending. 2013. Development of the Trona Deposit[M]. Beijing: China Petrochemical Press (in Chinese).
Zhang Tianfu, Zhang Yun, Cheng Xianyu, Sun Lixin, Cheng Yinhang, Zhou Xiaoxi, Wang Shaooyi, Ma Hailin, Lu Chao. 2020. Borehole databases and 3D geological model of Jurassic-Cretaceous strata in Dongsheng area, North Odors Basin[J]. China Geology, 47(S1): 220-245 (in Chinese with English abstract).
Zhong Yisi, Wang Licheng, Dong Haowei. 2022. Evaporite sedimentary characteristics and environment: A review[J]. Acta Sedimentologica Sinica, 40(5): 1188-1214 (in Chinese with English abstract).
陈建立. 2013. 泌阳凹陷碱矿床地质特征、成因及成矿预测——以安棚碱矿床为例[J]. 地质找矿论丛, 28(3): 393-400. https://www.cnki.com.cn/Article/CJFDTOTAL-DZZK201303010.htm 金强, 熊寿生, 卢陪德. 1998. 中国断陷盆地主要生油岩中的火山活动及其意义[J]. 地质论评, 44: 136-142. doi: 10.16509/j.georeview.1998.02.005 齐兵德, 谢星, 邓良武, 曾清亮, 马兰兰. 2011. 内蒙古苏尼特右旗查干诺尔天然碱矿成矿条件浅析[J]. 四川地质学报, 31(S1): 8-10. https://www.cnki.com.cn/Article/CJFDTOTAL-SCDB2011S1005.htm 王爱云, 陈文西. 2022. 天然碱: 从古代洗衣粉到化工之母[J]. 地球, (1): 6-11. https://www.cnki.com.cn/Article/CJFDTOTAL-DIQU202201003.htm 王九一, 刘成林, 王春连, 余小灿, 颜开, 高超. 2021. 晚白垩世-古近纪华南蒸发岩矿床形成的构造和气候耦合控制[J]. 地质学报, 95(7): 2041-2051. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXE202107004.htm 许虹, 彭齐鸣, Martin R Palmer. 2004. 辽宁古元古代地体中富电气石岩石的成因: 蒸发岩硼源的证据(英文)[J]. 中国地质, (3): 240-253. http://geochina.cgs.gov.cn/geochina/article/abstract/20040302?st=search 徐洋, 刘成林, 曹养同. 2021. 新疆库车盆地中始新世成盐演化特征——以KL4钻孔为例[J]. 地质学报, 95(7): 2183-2192. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXE202107017.htm 颜开, 刘成林, 王春连, 范美玲, 徐海明, 王九一. 2021. 刚果盆地西南部白垩纪蒸发岩矿物与古环境特征[J]. 岩石矿物学杂志, 40(3): 525-534. https://www.cnki.com.cn/Article/CJFDTOTAL-YSKW202103006.htm 杨江海, 易承龙, 杜远生, 张宗恒, 颜佳新. 2014. 泌阳凹陷古近纪含碱岩系地球化学特征对成碱作用的指示意义[J]. 中国科学: 地球科学, 44(10): 2172-2184. https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201410006.htm 叶铁林. 1978. 天然碱及其成因简述[J]. 化工矿山技术, 6: 18-23. https://www.cnki.com.cn/Article/CJFDTOTAL-HGKJ197806003.htm 叶铁林. 2013. 天然碱资源·地质·开采·加工第3版[M]. 北京: 化学工业出版社. 易承龙. 2016. 河南省泌阳凹陷安棚地区古近系核桃园组含碱地层层序特征及其意义[J]. 古地理学报, 18(1): 93-100. https://www.cnki.com.cn/Article/CJFDTOTAL-GDLX201601008.htm 张晨鼎. 1979. 天然碱矿床[J]. 纯碱工业, 3: 56-64. 张晨鼎. 2004. 土耳其贝帕扎里天然碱矿床的开发[J]. 纯碱工业, (2): 15-18. 张晨鼎. 2013. 天然碱矿床开发[M]. 北京: 中国石化出版社. 张天福, 张云, 程先钰, 孙立新, 程银行, 周小希, 王少轶, 马海林, 鲁超. 2020. 鄂尔多斯盆地北部东胜地区侏罗系-白垩系钻孔数据库与三维地质模型[J]. 中国地质, 47(S1): 220-245. http://geochina.cgs.gov.cn/geochina/article/abstract/2020S120?st=search 钟逸斯, 王立成, 董浩伟. 2022. 蒸发岩沉积特征及环境综述[J]. 沉积学报, 40(5): 1188-1214. https://www.cnki.com.cn/Article/CJFDTOTAL-CJXB202205003.htm