大森-宇津定律的一种可能机制以唐山大地震为例

为了探讨大森-宇津定律的物理机制, 本文在余震区等效黏度远低于其外部, 且构造应力场在整个余震活动时间间隔内基本保持不变的假设条件下, 提出了一个开尔文黏弹性地震震源体概念模型． 该模型可用于模拟主震后断层蠕变和震源区应力调整触发的余震序列以及蠕变停止后余震终结、 介质恢复到弹性状态、 断层重新闭锁和积累下一次地震的整个过程． 有限元方法可用来计算非均匀黏弹性地震震源体模型中主震和每次余震所引起的应力场及其随时间的演化过程． 在此基础上, 采用开尔文黏弹性地震震源体概念模型和有限元方法模拟了1976年唐山M_S7.8地震余震序列． 结果表明: 经验的大森-宇津定律可以用开尔文黏弹性震源体模型来解释, 这意味着余震衰减的频度取决于蠕变的速率; 余震序列持续时间受控于震源体的黏度, 即黏度越大, 蠕变时间越长, 余震持续的时间也就越长．      

利用GPS观测的时间序列资料反演地壳地幔黏性结构

地球介质的黏度是地球动力学研究中的重要参数. 在考虑震后震源区介质物性变化的情况下,本文提出利用GPS观测的时间序列资料反演地壳、地幔黏度的方法. 计算中,将遗传算法与Marquardt法相结合,利用遗传算法搜索结果作为Marquardt法反演的初始模型,运用微扰法形成Jacobian矩阵. 根据余震分布及地震波速度变化确定震源区的范围及介质的物性变化情况,利用黏弹性有限单元模型（Maxwell体）及台湾集集地震后GPS观测的时间序列数据,反演了台湾地区地壳、地幔的黏性系数. 反演计算及数值实验表明,该方法反演黏度的效率高,稳定性好. 反演结果初步显示台湾地区的地壳、地幔的黏度分别为12×1018Pa·s和36×1019 Pa·s.     

安徽铜陵地区燕山期侵入岩的成因及其对深部动力学过程的制约

铜陵地区燕山期侵入岩的岩石学、元素和同位素地球化学研究表明: (ⅰ) 研究区SiO₂≤55%的侵入岩主要为橄榄玄粗质系列岩石, 来自富集地幔的玄武质岩浆在上升过程中与下地壳物质发生低程度混染, 混染后的岩浆进一步发生分离结晶作用形成了这些岩石; (ⅱ) SiO₂ >55%的侵入岩主要为高钾钙碱性系列岩石, 与埃达克岩(adakite)有许多类似地球化学特征, 如富钠, 高Al₂O₃, Sr, Sr/Y与La/Yb比值, 大部分样品的Y < 18×10^-6, Yb < 1.9×10^-6, 但与埃达克岩也有不同之处, 如同位素组成((ε_Nd(t) = -9.16 ~ -16.55, (⁸⁷Sr/^86Sr)_i = 0.7068 ~ 0.7105)以及相当一部分样品的Y>18×10^-6, Yb>1.90×10^-6. 铜陵地区SiO₂ >55%的侵入岩很可能由幔源岩浆与玄武质下地壳熔融形成的埃达克质(adakite-like)岩浆混合形成. 来自地幔的橄榄玄粗质岩浆底侵可能为下地壳熔融提供了热量.     

利用ML振幅研究地壳横波Q值Ⅱ:Q横向变化特征

利用中国地震年报(1984～2004年)测定M_L所用的54000余条区域地震最大振幅数据,采用二维层析成像方法反演中国大陆地壳介质品质因子Q₀结构.反演结果表明,中国大陆地壳介质品质因子Q₀平均为370,Q₀横向变化量可达±200.在整个研究区解的分辨率能够达到3°×3°,大部分地区能够达到2°×2°.Q₀的变化显示出与地质构造分区相关的特征：低衰减 (高Q₀值) 异常区主要分布在塔里木盆地、四川盆地、鄂尔多斯高原及华南地区,鄂尔多斯高原北部Q₀值最高;高衰减异常区(低Q₀值)主要分布在青藏高原以东的川滇地区、华北平原、渤海湾,最低Q₀值在渤海湾地区,天山、祁连山和汾渭地堑Q₀也呈较低异常变化.总体看来,构造上稳定的克拉通地台的介质品质因子Q₀较高,而构造上活动的地区往往Q₀较低,显现出区域的衰减变化特征与地壳结构和区域构造之间的密切关系.     

扬子陆块西北缘碧口块体印支期花岗岩类地球化学和Pb-Sr-Nd同位素组成: 限制岩石成因及其动力学背景

对出露于扬子陆块西北缘碧口块体印支期阳坝岩体(215 Ma)、南一里岩体(224 Ma)和木皮岩体进行了岩石主量元素、微量元素和Pb-Sr-Nd同位素地球化学研究. 上述岩体花岗岩类均以高Al (Al₂O₃: 14.56~16.48%) 和Sr(352~1047 mg/g)、亏损Y(<16 mg/g)和HREE(eg. Yb<1.61 mg/g)为特征, 并具有较高的Sr/Y(36.3~150)和(La/Yb)_N(7.8~36.3)比值及强分异的稀土元素组成模式. 岩石初始Sr 同位素比值I_Sr=0.70419~0.70752, ε_Nd(t)=-3.1~-8.5,初始Pb同位素比值²⁰⁶Pb/²⁰⁴Pb=17.891~18.250, ²⁰⁷Pb/²⁰⁴Pb=15.494~15.575, ²⁰⁸Pb/²⁰⁴Pb=37.788~38.335. 地球化学特征显示阳坝、南一里和木皮岩体花岗岩类属于埃达克质(adakitic)岩石, 岩浆起源于增厚玄武质下地壳的部分熔融, 但它们具有较高的K含量(K₂O: 1.49%~3.84%)、明显演化的Nd同位素组成及较高的Nd同位素模式年龄(T_DM=1.06~1.83 Ga)清晰地不同于由俯冲洋壳或底侵玄武质岩石部分熔融形成的埃达克岩类, 而为增厚的并具有较长地壳存留年龄的玄武质下地壳部分熔融形成的埃达克质岩类. 碧口块体印支期埃达克质岩浆的产生反映了在华北板块和华南板块碰撞之后的岩石圈拆沉作用. 另一方面, 碧口块体印支期埃达克质岩石的Pb-Sr-Nd 同位素组成对岩浆源区的示踪揭示了在碧口块体的碧口群火山岩之下存在大陆型地壳基底, 这一结果不支持碧口群火山岩形成于大洋盆地或洋岛环境的认识.     

高温高压下辉长岩纵波速度和电导率实验研究

分别采用超声波透射-反射法和阻抗谱法在1～2 GPa, 室温到1100℃条件下测量了辉长岩的纵波速度(V_p)和电导率, 分析了影响辉长岩纵波速度的因素及其微观导电机制. 结果表明, 在800～850℃, 辉长岩的纵波速度开始大幅度下降, 波速下降受颗粒边界相、脱水熔融等因素的制约, 但在800～850℃, 辉长岩的电传导机制和电导率值却不会发生突变. 在680℃以下, 12～10⁵ Hz的频率范围内, 辉长岩阻抗谱上只出现代表颗粒内部传导机制的阻抗弧Ⅰ; 在680℃以上, 代表颗粒边缘传导机制的阻抗弧Ⅱ开始出现, 但岩石的总电导率由颗粒内部传导控制, 颗粒边缘传导对其影响很小. 实验测得的纵波速度与下地壳-上地幔平均V_p观测值一致, 电导率值与下地壳-上地幔辉长质岩石的组成相对应, 但与含有高导层的下地壳的大地电磁测量结果相差1～2个量级. 实验证实, 含水矿物脱水可诱发部分熔融, 通过这一机制可使地球内部的低速层和高导层具有一致性.     

华南白垩纪玄武质岩石的地球化学特征及其对太平洋板块俯冲作用的制约

通过研究白垩纪中晚期华南(浙、闽、赣、湘、粤五省)断陷红盆地和火山-沉积盆地中发育的同时代(80~110 Ma)玄武质岩石, 发现以武夷山为界华南白垩纪盆地中出露的玄武质岩石具有不同的地球化学特征. 西区白垩纪盆地(武夷山以西, 包括醴临-攸县盆地、衡阳、长沙-平江盆地、南雄、赣州、吉安等盆地)中的玄武质岩石具有低K₂O(K₂O=0.44%~3.17%, K₂O/Na₂O= 0.13~0.99)、低碱(K₂O+Na₂O=3.27%~6.80%)、低Al₂O₃(13.08%~16.75%)、高MgO(5.13%~8.78%)、高TiO₂ (1.12%~3.35%)的特征; 而东区(武夷山以东, 包括江西广丰、浙江玄坛地、文成周墩、永嘉、新昌、象山和福建永泰盆地)的玄武质岩石则以高K₂O(K₂O=0.55%~4.86%, K₂O/Na₂O= 0.19~2.22)、高碱(K₂O+Na₂O=2.95%~7.05%)、高Al₂O₃(15.80%~21.10%)、低MgO(2.63%~6.28%)、低TiO₂(1.19%~1.86%)为特征. 西区玄武质岩石富集高场强元素Nb和Ta, 而东区玄武质岩石富集大离子亲石元素K, Rb, Ba和Th, 亏损高场强元素Nb和Ta等; 稀土总量和轻重稀土比值东区(ΣREE=79.68~327.61 μg/g, (La/Yb) _N=5.35~26.37) 高于西区(ΣREE=74.66 ~287.72 μg/g, (La/Yb) _N = 4.49~22.10), 而(Gd/Yb)_N则是西区(1.30~4.62)高于东区(0.80~1.77). 西区玄武质岩石的(87Sr/86Sr) _i (平均值为0.704679)低于东区(平均为0.707658), 而ε _Nd(t)(&#8722;1.81~8.00)则高于东区(&#8722;8.50 ~ &#8722;1.22). 鉴此, 西区玄武质岩石可能形成于陆内裂解的构造环境中, 与软流圈的上涌有关, 源区有富集地幔(EMⅡ)和亏损地幔(DM)两端员混合的特征. 而东区玄武质岩石的源区则为受到了与板片俯冲有关的流体/熔体交代的岩石圈地幔, 与古太平洋板块的俯冲作用有关. 表明晚中生代存在古太平洋板块向亚欧大陆的俯冲作用, 但是俯冲作用的影响范围可能只限于武夷山以东地区.     

南秦岭早古生代地幔性质: 来自超镁铁质、镁铁质岩脉及火山岩的Sr-Nd-Pb同位素证据

南秦岭紫阳-岚皋地区早古生代晚期镁铁质岩脉及玄武岩的ε_Nd (t) = +3.28 ~ +5.02, (⁸⁷Sr/⁸⁶Sr)_i = 0.70341~0.70555, (²⁰⁶Pb/²⁰⁴Pb) _i = 17.256~18.993, (²⁰⁷Pb/²⁰⁴Pb) _i = 15.505~15.642, (²⁰⁸Pb/²⁰⁴Pb) _i = 37.125~38.968, Δ8/4 = 21.18~77.43, Δ7/4 = 8.11~18.82, 基本与南秦岭区新元古代中期以来的幔源岩石特征一致, 显示了HIMU, EMII和少量EMI富集地幔端元组分混合而成的Sr-Nd-Pb同位素组成特征, 表明与大洋地壳俯冲消减和陆缘物质再循环密切相关, 是新元古代早期扬子北缘大洋地壳俯冲消减及其携带的陆源沉积物再循环进入亏损软流圈地幔的结果.     

强震后地表变形的动力学机制研究—以1999年台湾集集地震为例

强震后地表变形的动力学机制是地球动力学研究的重要方面,现在普遍认为震后变形主要由断层的震后余滑或由介质的黏弹性松弛所至。1999年台湾集集地震GPS观测系统记录到了空前的资料,为研究震后变形的动力学机制提供了难得的机会。本研究认为集集地震后地表变形由震后断层余滑、下地壳／上地幔的黏弹性松弛、震源区介质的破裂、孔隙弹性回跳、地下流体的运移、介质孔隙度及孔隙压的变化等多种因素共同影响决定。为抓住重点,研究中将介质的破裂、地下流体的运移和孔隙弹性回跳等因素等效为震源区介质的物性变化。文中运用黏弹性有限单元模型（麦克斯威尔体）、利用GPS观测的时间序列资料对震后余滑、地壳／地幔黏度以及等效的震源区介质物性变化进行了反演。反演模型给出了震后余滑的分布及变化特征,反演结果初步显示台湾地区的下地壳／上地幔的黏度分别为2.7×10^18,4.2×10^20Pa·s。此外,反演结果还给出每种影响因素对地表变形的贡献大小,在集集地震后的450d时间里,断层的震后余滑引起的地表变形占总变形的44.6％,下地壳／上地幔的黏性松弛占34.7％,等效的震源区介质的物性变化占20.7％。     

岩石圈结构对大地震震后形变的影响——以1976年唐山大地震和2001年昆仑山大地震为例

大地震发生之后通常会诱发一系列的余震序列,对比1976年M_S7.8唐山大地震和2001年M_S8.1昆仑山大地震周边区域的地震事件可以看出,唐山大地震余震活动时间要明显长于昆仑山大地震余震活动时间.余震序列往往与震后形变密切相关,而影响震后形变的因素不仅与地震发震断层和震级有关,同时与岩石圈的结构有关.考虑到唐山大地震的发震区华北地块和昆仑山大地震的发震区青藏高原有着较大的岩石圈结构差异,本文采用PSGRN/PSCMP软件计算了岩石圈分层模型的大地震同震和震后形变,分析了地壳弹性模量、弹性厚度以及黏滞性系数对同震和震后形变的影响,进而讨论了影响唐山地震和昆仑山地震余震序列差异的原因.计算结果显示,震后形变会在黏弹性效应的作用下逐渐调整,震后形变的持续时间与地壳弹性模量、地壳弹性厚度和下地壳黏滞性系数有关.上地壳和下地壳弹性模量越大,震后形变达到稳定值的时间越短,弹性模量对震后形变稳定值影响很小.地壳弹性厚度越大,震后形变达到稳定值的时间越短,当断层面底端深度小于地壳弹性厚度时,地壳弹性厚度的增加会引起震后形变稳定值的减小;下地壳厚度对震后形变达到稳定值的时间和稳定值基本无影响.下地壳黏滞性系数越大,震后形变达到稳定值的时间越长,反之亦然.结合唐山地震区的华北地块和昆仑山地震的青藏高原深部结构发现,两者之间的上地壳弹性模型差别不大,唐山地震区地壳弹性厚度略大于昆仑山地震区,但昆仑山地震区下地壳黏滞性系数明显低于唐山地震区.这些因素均决定了昆仑山地震的震后形变持续时间短（余震时间序列短）而唐山地震的震后形变持续时间长（余震时间序列长）.由此可见,岩石圈结构差异可能是导致唐山地震和昆仑山地震余震序列差异的主要因素之一.     

华南中侏罗世玄武岩的岩石地球化学研究

华南发育一条近EW向展布的中侏罗世火山岩带, 福建永定盆地的“藩坑组玄武岩”是其中的典型代表之一. 藩坑玄武岩的TiO₂为1.92~3.21 wt%, 属“高Ti玄武岩”, 并有较高的全铁含量(11.09 wt%). 该玄武岩与东南沿海来自富集的岩石圈地幔的早白垩世玄武岩相比, 有完全不同的岩石地球化学特征: 相对富集Th, Nb, Ta, Zr和Ti等高场强元素, Zr/Ba, La/Nb, La/Ta, Zr/Y, Ti/Y, Ba/Nb, K/Ti和Rb/Zr等比值也与OIB相似; 除少数样品有较高ε_Nd(T)(1.87~3.55)外, 多数样品有类似球粒陨石的ε_Nd(T)(-0.70~0.24), 属软流圈地幔来源并记录着软流圈与岩石圈相互作用的信息. 赣南、粤北、湘东南的(早-)中侏罗世玄武岩、辉长岩与福建藩坑玄武岩有相似的OIB型岩石地球化学特征. 这种OIB型玄武岩可作为该地区早中生代时期软流圈地幔上涌的岩石学标志. 软流圈地幔的上涌导致该地区构造伸展、裂谷盆地的形成.     

希腊克里特MW6.0地震后的应力方向变化与强余震发生

地震的震源机制是地壳应力变化的指示器,而地壳应力变化与强震的发生直接相关。前人研究了地震震源机制变化在视应力较高的走滑型大震前的应力变化过程,而未见到震源机制变化对视应力较低的正断型大震发生的指示作用的研究。文章以2021年希腊克里特M_W6.0正断型地震序列为例,通过计算地震序列震源机制解与区域应力场方向之间最小空间旋转角的变化,揭示应力变化与强震发生的关系。为保证震源机制解的准确性,采用多家机构确定的震源机制得到中心震源机制作为该地震的震源机制,而后采用该地震序列精确的震源机制求解当地应力场,最后计算地震震源机制与主震震源机制及与所估计的地壳应力场方向的空间旋转角随时间的变化,探索强震发生与应力场变化的关系。结果表明:在主震发生的短期内,余震震源机制与该区域应力场方向的空间旋转角较大,与其后小震级的弱地震活动对应;随后余震震源机制与应力场方向的空间旋转角减小,对应后面发生的3次M_W＞5.0的强余震,在此之后的长时间内余震震源机制和应力场方向的空间旋转角再次增大,对应的余震震级及频度皆明显下降。文章以2021年希腊克里特M_W6.0地震序列为例,发现视应力较低的正断型地震前也存在应力方向集中现象,为探索地震应力前兆提供了范例。     

中国大陆及邻近地区Lg尾波的Q值分布

10个中国数字地震台网(CDSN)台站和5个全球地震台网(GSN)台站所记录的785个浅源地震的宽带垂直分量, 被用来研究中国大陆及邻近地区L_g尾波的衰减特性. 首先运用叠加谱比法对各记录进行处理, 得到与各路径相对应的椭圆内L_g尾波的Q₀ (1 Hz处的Q值)和频率相关因子η的平均值. 进一步应用反投影技术, 得到L_g尾波的Q₀值和η值的成像图及其误差分布. 结果表明, 在所研究的范围内, Q₀在200~500的范围内变化. 最低的Q₀发生在滇藏地区; 最高的Q₀发生在西伯利亚地台的南端. η值的变化范围在0.3~0.8之间. 对于所研究的大部分地区, η值呈现出与Q₀值的反向变化关系.     

佛冈花岗岩基及乌石闪长岩-角闪辉长岩体的形成年龄和起源

佛冈花岗岩基约6000 km², 是南岭地区最大的晚中生代岩基. 乌石闪长岩-角闪辉长岩岩体位于佛冈花岗岩基的东北部, 它和佛冈花岗岩主体一并属于高钾钙碱性系列岩石. 但乌石岩体以低Si (49%~55%), 高Fe, Mg, Ca, 稀土总量低, Eu, Ba, P, Ti亏损不明显, 而Zr和Hf亏损明显的特征不同于佛冈花岗岩主体. 锆石LA-ICP-MS定年及矿物-全岩Rb-Sr等时线年龄测定结果表明, 乌石岩体的结晶年龄为160 Ma左右, 与佛冈花岗岩主体是同时代形成的. 佛冈花岗岩主体具有较高的(⁸⁷Sr/⁸⁶Sr)_I值(0.70871~0.71570), ε_Nd(t)变化于-5.11~-8.93之间, 显示出壳源花岗岩的Sr-Nd同位素特点, 它们的两阶段Nd模式年龄介于1.37~1.68 Ga. (⁸⁷Sr/⁸⁶Sr)_I值、ε_Nd (t)值和Nd模式年龄值的不均匀性, 可能反映巨大的佛冈花岗岩主体的源区组成是不均匀的, 同时在其形成过程中有地幔物质的不均匀混合. 乌石闪长岩-角闪辉长岩是一种少见的高(⁸⁷Sr/⁸⁶Sr)_I值(0.71256~0.71318)、低ε_Nd (t值(-7.32 ~ -7.92)中基性岩浆岩, 它可能由地幔部分熔融产生的新生幔源玄武质岩浆与下地壳玄武质岩石脱水部分熔融产生的岩浆混合形成.     

西秦岭花岗岩类地球化学和Pb-Sr-Nd同位素组成对基底性质及其构造属性的限制

对西秦岭印支期5个花岗岩类岩体进行了岩石主量元素、微量元素和Pb-Sr-Nd同位素地球化学的研究, 据此限定西秦岭的地壳基底性质及其构造属性. 结果表明, 西秦岭花岗岩类总体化学成分偏基性, 岩石主要属于准铝质的高钾-钾玄质系列, A/CNK=0.90~1.05 (绝大多数样品<1.0), K₂O/Na₂O=1.04~1.86. 它们具有相似的微量元素(包括稀土元素)组成模式, (La/Yb)_N= 7.49~ 28.79, Eu*/Eu=0.39~0.76. 在Sr-Nd同位素组成上, 西秦岭花岗岩显示一定程度的不均一性, I_Sr= 0.70682~0.70845, ε_Nd(t)=-4.85~-9.17, T_DM=1.26~1.66 Ga. 西秦岭花岗岩类以高放射成因铅同位素组成为特征, 其初始铅同位素比值为: ²⁰⁶Pb/²⁰⁴Pb=17.996~18.468, ²⁰⁷Pb/²⁰⁴Pb=15.565~15.677, ²⁰⁸Pb/²⁰⁴Pb=38.082~38.587. 根据西秦岭花岗岩类的化学和Sr-Nd同位素组成, 揭示了它们的岩浆源区均来自于地壳中高K(Rb)玄武质岩石的部分熔融, 源区物质形成时代可能在900~1400 Ma之间, 由此反映在西秦岭沉积盖层之下含有大量的中、新元古代的高K (Rb)玄武质岩层. 西秦岭印支期花岗岩类与东秦岭印支期花岗岩类的Pb-Sr-Nd同位素组成对比, 指示西秦岭和东秦岭地壳具有不同的基底组成, 两者的分界线大至位于近南北向的宝成铁路线. 西秦岭花岗岩类Pb-Nd同位素组成特征表明西秦岭地壳基底具有扬子块体的构造属性.     

2019年12月26日湖北应城M4.9地震震源参数及余震活动性研究

2019年12月26日湖北应城发生M4.9有感地震,其震感波及武汉大部分地区。为了分析该地震的发震构造及余震活动性,本文利用波形拟合方法测定了不同速度模型下该地震的震源机制解和矩心深度,并用Bootstrapping抽样反演技术评价反演结果;此外,利用模板匹配技术匹配主震和目录余震波形,获取了更为完整的余震目录。结果显示,应城地震以走滑为主,矩心深度7.5km左右,矩震级M_W4.67;应城地震有1个前震和17个余震,余震序列缺少M2~4事件,表明应城地震为孤立型地震,M2以下地震的b值为0.8。     

水位潮汐周日波高值异常与中强地震的关系

地球对引潮力产生的响应反映了其内部的物性特征.潮汐观测值反映地壳在潮汐作用下发生的微形变场变化,通过对地壳介质物性参数(潮汐因子、相位滞后)的计算,分析地震前地壳岩石弹性变化,判定其由线性进入非线性状态的短期阶段.这一理论有助于追踪震源区的非线性应变过程,从而判定非线性系统失稳前兆.     

2015年新疆皮山6.5级地震序列震源机制特征分析

利用新疆区域数字地震台网波形资料,采用CAP方法反演了2015年7月3日新疆皮山6.5级地震主震及部分M_S≥3.6余震的震源机制解和距心深度。研究结果显示,皮山6.5级地震主震最佳双力偶解节面I：走向290°/倾角55°/滑动角96°;节面Ⅱ：走向101°/倾角35°/滑动角82°,最佳矩心深度16km,表明该地震是一次逆冲型事件。通过反演部分M_S≥3.6余震的震源机制解发现,早期余震的破裂方式与主震较为一致,随着时间的推移余震震源机制出现走滑型和正断型,表明早期余震的破裂受主震影响较大,随着序列的发展变化,后期震源区应力场可能出现一定程度的调整。统计皮山6.5级地震序列P轴方位发现,优势方位为NNE向,与该区域构造应力场方向较为一致。结合地震序列的震源机制及他人精定位结果和震源区地质构造情况,初步解释了导致此次地震的原因。     

根据大地震前后应力轴偏转和应力降求取应力量值的研究

为求解震源处应力量值,运用大地震造成的应力方向偏转和地震应力降给出了求取震源处应力大小的方法.地震前后的应力方向可根据大量前震和余震震源机制统计给出;地震应力降可通过地震波或地壳形变资料的震源破裂反演研究得到.该研究对地震动力学具有一定的意义.      

高演化海相碳酸盐烃源岩地球化学综合判识

针对中国海相碳酸盐岩烃源岩的成熟度高、有机质丰度低等特殊性, 以鄂尔多斯盆地下古生界为例, 通过微量元素富集系数计算、稀土元素的配分模式、元素和碳同位素等的数理统计分析及三维图解等多元综合分析, 结合必要的有机参数, 反演了盆地烃源岩发育环境, 评价碳酸盐岩烃源岩及其分布. 就鄂尔多斯盆地而言, TOC为0.2%是个重要的界限值, 研究表明: 在TOC>0.2%的层段中, Ba富集, δ¹³C_carb正偏移, δ¹³C_org<&#8722;28‰, 反映出古生产力、有机质埋藏量较高, 具较好生烃潜力, 稀土元素分异强, 指示沉降速率低, 有机质保存条件好, 利于烃源岩发育. 克里摩里组、乌拉力克组具一定规模, 可视为潜在有效烃源岩. 而TOC＜0.2%的层段, Ba含量低, δ¹³C_carb多为负值, δ¹³C_org为&#8722;24‰~&#8722;28‰, 指示有机质埋藏量低, 生烃潜力差, 稀土元素分异弱, 沉降速率较快, 多属浅水高能氧化环境, 不利于烃源岩发育. 利用无机与有机地球化学方法综合评价高演化海相有潜力的碳酸盐岩烃源岩及分布是可行的.