地庙算法

‘壹’ 请教达人看看星盘，落陷、接纳、庙旺啥的

你好!
你问:请教达人看看星盘，落陷、接纳、庙旺啥的?
这是我一篇一直在用!但大多数都能接受的回答!
送给你也送给所有想算命运的人:

大多数的人往往都是在生活中!遇到不顺或对自己的未来的人生好奇或疑惑时最喜欢算命!
所以…..建议如下:
诸如道家还有西洋星象学.......都有在帮人算命或算运!

假设:今有两人帮你算星座十二生肖八字命运!一人解说你命好~运好~贵人多!命中煞中带劫!须要改运!否则………..这时那你应该要来相信那位呢?

其实真正的八字算法是很严谨及复杂!还必须从算者的周遭包括!父母!配偶!面相!手相..经度~纬度..............
因为同样的出生八字!同样的姓名!都存在不同的命运结果!
希望你明白我在说什么?

人都是喜欢听好听的!
更耽心的是回答的人也不知道真懂假懂!这不是害人吗?
所以把心结打开!别管它!个性合得来比星座八字合得来更为重要!

太过在意算命命理运势所谓的暗示! 只会害自己陷入往后生活中思想上........好坏阴影包袱之中!

大多数的宗教大都劝善止恶!主要带给人们的是心灵上的寄托!
但佛教佛法中讲求靠的主要是自修重因果!不谈怪力乱神!包括抽签~解梦~算命或算运~改运!还有撞邪!避邪之说!

又假设:某人的小孩正在学校读书!他希望他的小孩在班上读书考试能名列前矛!于是准备丰盛供品到处求神拜佛!拜托神佛能保佑他的小孩考试能考第一名!你觉得这样有效吗?佛教佛法中的种种教法教导告诉我们!种什么因得什么果!他的小孩考试要考第一名!主要靠的还是只有他自己!别人是帮不了他的!父母能做的也只能是给他好一点的学习环境!

所以不管是从科学的角度或从宗教的角度来看待命运!
只要自己能常持善心~多行善事!
自己命运!~自己决定!
随顺因缘~把心放下!
就能活得自然!
过得自在!
但愿能帮到你!
祝愿你!
吉祥!如意!

南无法阿弥陀佛

ps:最后建议你~不管别人给你算得结果如何?希望你能理性看待!!
如果你觉得能接受也认同我以上的建议说法!还望请您及时给我采纳满意答案鼓励!
让我保持更好的动力!来帮助更多的人!谢谢! 祝福你!

‘贰’ 汉诺塔的算法

算法介绍：当盘子的个数为n时，移动的次数应等于2^n–1。后来一位美国学者发现一种出人意料的简单方法，只要轮流进行两步操作就可以了。首先把三根柱子按顺序排成品字型，把所有的圆盘按从大到小的顺序放在柱子A上，根据圆盘的数量确定柱子的排放顺序：若n为偶数，按顺时针方向依次摆放A、B、C；

若n为奇数，按顺时针方向依次摆放A、C、B。

所以结果非常简单，就是按照移动规则向一个方向移动金片：如3阶汉诺塔的移动：A→C,A→B,C→B,A→C,B→A,B→C,A→C

汉诺塔问题也是程序设计中的经典递归问题。

(2)地庙算法扩展阅读

由来：

法国数学家爱德华·卢卡斯曾编写过一个印度的古老传说：在世界中心贝拿勒斯（在印度北部）的圣庙里，一块黄铜板上插着三根宝石针。印度教的主神梵天在创造世界的时候，在其中一根针上从下到上地穿好了由大到小的64片金片，这就是所谓的汉诺塔。

不论白天黑夜，总有一个僧侣在按照下面的法则移动这些金片：一次只移动一片，不管在哪根针上，小片必须在大片上面。僧侣们预言，当所有的金片都从梵天穿好的那根针上移到另外一根针上时，世界就将在一声霹雳中消灭，而梵塔、庙宇和众生也都将同归于尽。

不管这个传说的可信度有多大，如果考虑一下把64片金片，由一根针上移到另一根针上，并且始终保持上小下大的顺序。这需要多少次移动呢?这里需要递归的方法。假设有n片，移动次数是f(n).显然f(1)=1,f(2)=3,f(3)=7，且f(k+1)=2*f(k)+1。此后不难证明f(n)=2^n-1。n=64时，

假如每秒钟一次，共需多长时间呢？一个平年365天有31536000 秒，闰年366天有31622400秒，平均每年31556952秒，计算一下：18446744073709551615秒。

这表明移完这些金片需要5845.54亿年以上，而地球存在至今不过45亿年，太阳系的预期寿命据说也就是数百亿年。真的过了5845.54亿年，不说太阳系和银河系，至少地球上的一切生命，连同梵塔、庙宇等，都早已经灰飞烟灭。

‘叁’ 房产税的计算方法怎么算的

房产税是以房屋为征税对象，按房屋的计税余值或租金收入为计税依据，向产权所有人征收的一种财产税。征收范围限于城镇的经营性房屋。区别房屋的经营使用方式规定征税办法，对于自用的按房产计税余值征收，对于出租房屋按租金收入征税。根据税法规定，房产税的计算方法有以下两种：1.按房产原值一次减除30%后的余值计算。其计算公式为：年应纳税额=房产原值×(1-30%)×1.2% ;2.第二种按租金收入计算，其计算公式为：年应纳税额=年租金收入×适用税率(l2%)。下列房产免纳房产税：一、国家机关、人民团体、军队自用的房产;二、由国家财政部门拨付事业经费的单位自用的房产;三、宗教寺庙、公园、名胜古迹自用的房产;四、个人所有非营业使用的房产;五、经财政部批准免税的其他房产。

‘肆’ 深度优先算法的例1.迷宫问题

实现这一算法，我们要用到编程的一大利器--递归。“递归”是一个很抽象的概念， 但是在日常生活中，我们还是能够看到的。拿两面镜子来，把他们面对着面，你会看到什么？你会看到镜子中 有无数个镜子？怎么回事？A镜子中有B镜子的象，B镜子中有A镜子的象，A镜子的象就是A镜子本身的真实写照，也就是说A镜子的象包括了A镜子，还有B镜子在A镜子中的象………………好累啊，又烦又绕口，还不好理解。换成计算机语言就是A调用B，而B又调用A，这样间接的，A就调用了A本身，这实现了一个重复的功能。
再举一个例子：从前有座山，山里有座庙，庙里有个老和尚，老和尚在讲故事，讲什么呢？讲：从前有座山，山里有座庙，庙里有个老和尚，老和尚在讲故事，讲什么呢？讲：从前有座山，山里有座庙，庙里有个老和尚， 老和尚在讲故事，讲什么呢？讲：…………。好家伙，这样讲到世界末日还讲不完，老和尚讲的故事实际上就是前面的故事情节，这样不断地调用程序本身，就形成了递归。 万一这个故事中的某一个老和尚看这个故事不顺眼，就把他要讲的故事换成：“你有完没完啊！”，这样，整个故事也就嘎然而止了。我们编程就要注意这一点，在适当的时候，就必须要有一个这样的和尚挺身而出，把整个故事给停下来，或者使他不再往深一层次搜索，要不，递归就会因计算机存储容量的限制而被迫溢出，切记，切记。 这是一道很经典的题目了，我们先要明确一下思路，如何运用深度优先搜索法，完 成这道题目。我们先建立一个8*8格的棋盘，在棋盘的第一行的任意位置安放一只皇后。紧接着，我们就来放 第二行，第二行的安放就要受一些限制了，因为与第一行的皇后在同一竖行或同一对角线的位置上是不能安放 皇后的，接下来是第三行，……，或许我们会遇到这种情况，在摆到某一行的时候，无论皇后摆放在什么位 置，她都会被其他行的皇后吃掉，这说明什么呢？这说明，我们前面的摆放是失败的，也就是说，按照前面 的皇后的摆放方法，我们不可能得到正确的解。那这时怎么办？改啊，我们回到上一行，把原先我们摆好的 皇后换另外一个位置，接着再回过头摆这一行，如果这样还不行或者上一行的皇后只有一个位置可放，那怎 么办？我们回到上一行的上一行，这和老鼠碰了壁就回头是一个意思。就这样的不断的尝试，修正,我们最终会得到正确的结论的。

‘伍’ 到庙上做佛事，算法布施吗

——感恩您对我的信任，我将悉心为您解答：
要明白这个问题，先了解一下什么是“六度波罗蜜”。佛教中有六种非常重要的修行方法称为“六度波罗蜜”（简称就是“六度”），“度”可以形象的理解为“渡河”的“渡”，佛教将我们陷身于情绪烦恼的境界称为“烦恼的彼岸”，而如何摆脱情绪烦恼的纠缠的境界称为“觉悟的彼岸”，修行就是从烦恼的彼岸达到觉悟的彼岸的过程，那么六度就是六个到彼岸的方法，犹如乘船过河一样。这六种方法分别称为“布施、持戒、忍辱、精进、禅定（止观）、般若”。
布施：布施是六度波罗蜜中的第一课，也是修行中最重要的基础。布施，即“舍弃”的意思，我们人性的本质是贪婪的，佛教称为“贪、嗔、痴”，不好的东西可以舍弃，但是好的东西、珍贵的东西就舍不得放弃，因此，学会布施，就是可以克服我们内心的贪婪和吝啬，对财物的执着，可以增长我们的福德，增加我们此生的财富，还可以带到来世去，中国的“舍得”的观念，“大舍有大得，不舍就不得”，其实就有佛教的布施的成分，但是范围窄了一些，不能替代“布施”这个概念。
布施的种类：布施，布施单纯的给予财物，它包含了“财布施”、“身布施”、“法布施”、“无畏布施”四大类型，简单的理解，“财布施”就是把自己的财产包括吃的用的住的等捐赠出来，如同512汶川地震大家踊跃捐出衣服、饮食用品等一样，“身布施”就是捐献自己的身体，小处可以利用自己的节假日来寺院做义工，或为慈善公益事业服务，大处可以捐献自己的身体，如眼角膜、献血献骨髓、移植内脏手脚乃至整个生命，这也是不容易做到的。“法布施”就是将佛法告诉别人，小则教会别人生存的技能，大则解脱别人的人生痛苦。“无畏布施”则是坚持自己的正确的行为，作为启迪别人的榜样，最重要的方式就是持“五戒”，五戒是佛教中的基本戒，近似的理解为“五种最基本的规章制度”，包括了“不杀生”、“不偷盗”、“不邪淫”、“不妄语”、“不饮酒”。当我们不杀生的时候，别人就不必担心被我们杀了。我们不偷盗，别人就不必时时刻刻提防着我们。我们不邪淫，别人就不必担心我们会强暴他的妻子。我们不妄语，别人就不必害怕被我们坑蒙拐骗。我们不饮酒，别人就不怕我们会佯装酒醉地胡搅蛮缠。
所以佛教里有这样一句话：“功德林中，布施第一”，说明了布施在整个修行的过程中占据着第一重要的位置。
——接下来，你就应该知道，你的布施是属于哪一类了吧？

‘陆’ 房产税的详细计算方法

计算方法：

1、从价计征

从价计征是按房产的原值减除一定比例后的余值计征，其公式为：

应纳税额=应税房产原值×（1-扣除比例）×年税率1.2%

2、从租计征

从租计征是按房产的租金收入计征，其公式为：应纳税额=租金收入×12%

没有从价计征的换算问题纳税义务发生时间：将原有房产用于生产经营，从生产经营之月起，缴纳房产税。其余均从次月起缴纳。

个人出租住房的租金收入计征，其公式为：应纳税额=房产租金收入×4%

(6)地庙算法扩展阅读：

根据《中华人民共和国房产税暂行条例》：

第四条房产税的税率，依照房产余值计算缴纳的，税率为1.2%；依照房产租金收入计算缴纳的，税率为12%。

第五条下列房产免纳房产税：

一、国家机关、人民团体、军队自用的房产；

二、由国家财政部门拨付事业经费的单位自用的房产；

三、宗教寺庙、公园、名胜古迹自用的房产；

四、个人所有非营业用的房产；

五、经财政部批准免税的其他房产。

第二条房产税由产权所有人缴纳。产权属于全民所有的，由经营管理的单位缴纳。产权出典的，由承典人缴纳。产权所有人、承典人不在房产所在地的，或者产权未确定及租典纠纷未解决的，由房产代管人或者使用人缴纳。

前款列举的产权所有人、经营管理单位、承典人、房产代管人或者使用人，统称为纳税义务人（以下简称纳税人）。

第三条房产税依照房产原值一次减除10%至30%后的余值计算缴纳。具体减除幅度，由省、自治区、直辖市人民政府规定。

没有房产原值作为依据的，由房产所在地税务机关参考同类房产核定。房产出租的，以房产租金收入为房产税的计税依据。

‘柒’ 页岩气储层矿物组分及有机碳含量测井评价方法研究——以鄂西渝东建南构造东岳庙段为例

路 菁^1，2 李 军¹

（1.中国石化石油勘探开发研究院，北京 100083；

2.中国石油大学（北京）博士后流动站，北京 102249）

摘 要 鄂西渝东地区下侏罗统为四川盆地典型的陆相页岩气藏，储层矿物组分及有机碳含量是确定该类气藏工程开采难度与有效性的重要指标。为突破常规储层测井评价方法在复杂矿物储层评价中存在的多解性问题，本研究充分挖掘常规测井资料中蕴含的地质信息，以非线性反演与最优化算法为核心思想，综合评价包含有机碳在内的页岩岩石组分与含量，取得了较好的测井评价结果。研究结果完善了页岩气储层测井评价手段，为推进页岩气勘探开发相关技术发展起到了积极的作用。

关键词 常规测井响应 矿物组分 有机碳含量 非线性反演 最优化方法 测井评价

Logging Evaluation of Mineralogical Constituent and

Total Organic Contents for Gas Shale

LU Jing^1，2，LI Jun¹

（1.Exploration and Proction Research Institute，SINOPEC，Beijing 100083，China；

2.Postdoctoral Center，China University of Petroleum，Beijing 102249，China）

Abstract Formation of lower Jurassic in the western Hubei and East Chongqing is an typical continental facies shale gas reservoir in Sichuan Province.The mineral constituents and total organic contents（TOC）are important indicator of the engineering difficulty and its effectiveness for such gas reservoir.To breakthrough the problem of multi －solutions，that always occur when the conventional reservoir logging evaluation methods are used to solve the gas shale reservoir evaluation，this study fully tap the geological information concealed in conventional logging response，use the nonlinear joint inversion and optimization as the core ideas，to evaluate both the mineralogical constituent contents and TOC for gas shale，and achieve a desirable result .This research supplements the logging evaluation methods for gas shale and play a positive role in related technology developments for gas shale exploration and development.

Key words conventional logging；mineralogical constituent；TOC；nonlinear joint inversion；optimization； logging evaluation

鄂西渝东地区是四川盆地周缘页岩气藏有利目标区之一。建南构造位于四川盆地川东褶皱带石柱复向斜中北部，下侏罗统自流井组发育的深湖－半深湖页岩属于典型的陆相页岩气藏。该套页岩区域分布稳定、厚度大、埋藏浅，但相较于海相页岩具有更加频繁的相变特征，储层矿物组分复杂多变。准确把握页岩气储层矿物组分与有机碳含量是后续储层关键参数——脆性与含气性评价的重要基础，也是页岩气测井评价亟待攻克的重点及难点问题。依靠固定的解释模型，采用少部分测井曲线确定储层矿物含量的评价方法，在岩性及矿物较为单一的常规储层中评价效果较好，却无法妥善解决页岩气储层复杂矿物组分与含量的多解性问题。笔者通过深入挖掘各项常规测井资料中蕴含的丰富地质信息，分析建立储层矿物组分模型，以非线性反演与最优化算法评价包含有机碳在内的复杂岩石组分含量，突破了常规测井储层评价的思想，拓展了非常规页岩气储层矿物组分与有机碳含量的测井评价方法，通过实验室岩心全岩组分数据验证，该方法已取得了较好的评价效果。

1 东岳庙段含气页岩岩性及岩石矿物学特征

目标层下侏罗统自流井组东岳庙段泥页岩，区域横向分布稳定，厚度较大，暗色泥页岩厚约60～100m。储层岩性以含灰泥页岩为主，多见灰色粉砂质泥页岩、介壳泥页岩与介壳灰岩夹层（图1）；储层矿物成分以黏土矿物、石英及方解石为主（平均含量分别为22.49％、55.95％、17.5％），同时含有少量长石与黄铁矿。自生矿物的存在，表明东岳庙段所处的沉积环境为有利于有机质富集与保存的还原环境，实验室分析结果显示，储层有机碳以Ⅱ型干酪根为主，平均含量2％～3％；储层孔隙结构以矿物粒间孔为主，同时发育少量粒内孔及溶蚀裂隙，大量因有机质热解产生的纳米孔隙，使储层具有较好的天然气吸附与储集性能。

图1 研究区东岳庙段泥页岩典型岩性

2 常规测井响应评价储层岩石组分

测井响应是被测地层物理特性的宏观表现^［1］，在排除井眼与泥浆侵入等影响的情况下，测井响应本质是测井仪器探测范围内所有岩石微观组分物理特性的综合表现，故各类测井响应实际上涵盖了被测地层所有组分的岩石物理信息。充分挖掘、利用常规测井响应中蕴藏的储层信息评价页岩岩石组分，提供了一条除实验室分析和元素俘获能谱（ECS）测井之外的储层评价思路，同时，弥补了岩心实验室分析无法全井段连续、ECS测井数据采集与解释评价成本高昂等问题^［2，3］。

2.1 常规曲线非线性联合最优化反算法

2.1.1 目标函数

区别于利用单一或少数测井曲线与储层某一矿物含量建立函数关系、用以评价其含量的方法，利用常规测井信息开展非线性联合最优化反演评价储层矿物组分的方法与步骤，可简要概括如下：首先，需要对实测响应进行预处理，以期得到接近原始储层真实物理特性的校正测井响应；其次，依据岩心观察与常规评价结果得到的初步认识，圈定解释评价井段内存在的岩石组分类型，并确定其初始含量，形成完整的基于原始假设的储层岩石物理体积模型；再次，依据地区经验或理论参数合理选取各组分的测井响应骨架值，以非线性测井响应方程正演各个常规测井响应，并计算关于校正曲线与正演模拟曲线如式（1）所示的目标函数T（X^j）；最后，通过反复迭代调整各矿物组分含量，使目标函数T（X）达到最小值，并将此时的岩石组分与含量模型作为反演的最终结果，即通过解决图2所示的最优化问题，达到求解复杂矿物储层岩石组分与含量问题的目的^［4］。

油气成藏理论与勘探开发技术（五）

图2 非线性联合最优化反算法简图

式中：logging_s为第j次迭代后产生的正演曲线组；logging_c为实测曲线经校正产生的校正曲线组；X^j为第j次迭代确定的各个岩石组分含量；W为各测井曲线在目标函数中的权重；α为迭代稳定性控制参数；T（X^j）为反映正演曲线与校正曲线相似程度的目标函数，当该函数达到最小值时，表明正演曲线已逼近校正曲线，此时，即可认为模型求解得到的岩石组分与含量与地层真实情况最为接近。需要说明的是，采用更丰富的测井响应信息，以及岩心分析、常规储层评价取得的地层初步认识等，能够在更大的程度上降低反算法的多解性。

2.1.2 共轭梯度最优化算法

从上述分析可知，求解页岩复杂岩石组分的测井评价问题，已被转化为求解目标函数T（X^j）最小值的最优化问题。本研究综合考虑目标函数属于多元函数，且测井响应的非线性关系决定了目标函数的非线性特性，故采用共轭梯度法解决目标函数的最优化问题^［5］。

对目标函数T（X^j），在极值点X^*处作Taylor展开，忽略高效项时，有

油气成藏理论与勘探开发技术（五）

式中：H＝▽²T（X^*）为T（X）在X^*处的二阶偏导数矩阵。因为X^*为极值点，故▽T（X^*）＝0，因而

油气成藏理论与勘探开发技术（五）

可见，任何次的函数T（X）在其极值点附近具有二次函数的特征。设T（X）可以表示为如下所示二次函数

油气成藏理论与勘探开发技术（五）

可以证明具有N阶正定矩阵A的n元二次函数，最多可在n维空间中找到n个彼此关于A的共轭方向（向量），且从任意的初始点出发，依次沿这n个共轭方向作不超过n次的一维搜索，就可以求得目标函数T（X）在n维空间内的极小点。采用上述共轭梯度算法回避了因牛顿法及其改进算法需要计算二阶偏导数矩阵的逆矩阵而带来的巨大运算量，且克服了最速下降法在接近极小点时收敛速度很慢的缺陷，妥善地解决了研究建立的非线性反算法的求解问题。

2.2 东岳庙段页岩岩石组分反演

2.2.1 初始模型假设

图3为研究区某井东岳庙段泥页岩常规测井响应，该井含气页岩岩石组分评价的目的在于，明确包括有机碳在内的岩石重要组分的具体含量。初始模型假设的建立，需要分别确定待求解的储层岩石组分及其初始含量，以及参与岩石组分评价的测井曲线。

依据上节所述实验室全岩分析结果（图3），初始模型假设页岩中不存在除干酪根之外的其他固体有机碳；脆性矿物包括石英、方解石、长石，塑性矿物即为黏土；另外，由于相关研究表明，页岩成岩过程中自生的黄铁矿常结晶于储层层理界面之间，在一定程度上有利于水力压裂形成网状缝，且黄铁矿物具有极好的导电特性、极高的光电俘获截面指数以及较高的密度，即使含量较小，对电阻率、光电截面指数与体积密度等测井响应的影响也十分明显，因此，作为影响页岩力学性质与岩石物理特性的重要矿物，黄铁矿在岩石组分模型中不可忽略；最后，由于该段泥页岩黏土矿物含量较高、有效孔隙度较低，且地层水矿化度不高，自由水对测井响应影响不大，故模型仅考虑黏土束缚水存在且假设页岩储层有效孔隙全部被游离气占据的情况。综合上述考虑，最终确定该井东岳庙段泥页岩需要反演计算的岩石组分如图4所示，依次包含黏土（含黏土束缚水）、石英、方解石、长石、黄铁矿、孔隙（游离气）与有机碳（干酪根）。

图3 建南地区某井东岳庙段泥页岩常规测井响应特征与岩心分析结果

图4 页岩岩石体积模型

综合考查本井可参考的测井曲线条数，以及上述页岩岩石体积模型需要涵盖的组分种类，确定利用光电截面指数（PEF）、自然伽马（GR）、中子孔隙度（NPHI ）、体积密度（DEN）、声波时差（DT）、浅侧向电阻率（LLS）、深侧向电阻率（LLD）、铀（URAN）、钍（TH）共9条曲线（图3），反演8种岩石（图4中）组分的含量。可以注意到，如不考虑欠定求解，参加非线性反演的测井曲线条数理论上最多可处理10种岩石组分含量的求解问题，此数大于本模型求解的岩石组分数量，故模型求解结果属于非线性超定解，能够有效降低评价结果的多解性，确保评价结果更加接近页岩气储层的真实情况。

依靠常规储层评价方法，如自然伽马泥质含量V_sb评价方法^［6］、密度中子孔隙度P_hi评价方法^［7］、Pessay有机碳TOC含量评价方法等^［8］，可以取得黏土、孔隙度、干酪根含量的初步评价结果，对本井岩石组分初始含量 进行赋值， 剩余组分的初始含量——石英含量 根据实验室岩心分析确定的平均含量（石英55.9％、方解石17.5％、长石7.2％、黄铁矿4.5％）按比例分配，结果如图5第2～8道内实线所示。可以注意到，各岩石组分初始含量（棕色实线）与岩心分析结果（黑色圆点）相比，均存在不同程度的偏差。其中方解石、长石两种矿物含量的偏差最为明显；利用中子－密度孔隙度评价的孔隙度结果也明显偏高；此外，利用自然伽马泥质含量评价方法计算的黏土矿物含量，以及电阻率－声波重叠Passey法计算的有机碳含量，在局部深度上还存在一定误差。本研究将通过随后的反演计算逐步降低这些误差，以得到最接近真实地层岩石组分的评价结果。

图5 建南地区某井东岳庙段泥页岩非线性反演初始模型

2.2.2 模型反演结果

经过非线性反演计算，最终确定该井岩石组分的含量如图6所示，图中第2～8道依次为黏土矿物（含黏土束缚水）、石英、方解石、长石、黄铁矿、孔隙及有机碳含量的评价结果（实线）与对应组分实验室分析结果（黑色圆点），图中第9与第10道分别为页岩岩石组分非线性反演结果与岩心实验室分析结果。

图6 建南地区某井东岳庙段泥页岩岩石组分非线性反演成果图

通过图7各组分初始评价结果（黑色方块）与非线性反演计算结果（三角）的对比分析可以发现，非线性反演结果与实验室分析结果具有更好的线性相关性，与初始评价结果相比更集中于45°对角线附近。图6与图7均显示，非线性反算法显着提高了石英与方解石含量的评价精度；使孔隙度评价结果更加接近实验室分析结果；此外，黏土矿物与有机碳含量各自在局部位置上的误差也得到了较好的修正；在初始模型中，以平均含量为依据粗略估算的长石与黄铁矿含量，这里也得到了进一步细化，评价结果与实验室分析结果在整体趋势上更为吻合。至此，本研究利用建立的非线性反演方法，同时完成了研究区东岳庙段页岩气储层复杂矿物组分与有机碳含量测井评价两个问题，且取得了较高的评价精度，本研究将进一步定量分析测井评价结果，以验证该方法的可靠性与有效性。

2.2.3 非线性反演结果分析

考虑到各项实测测井响应其本质是被测储层岩石组分反映在各类物理场中的宏观物理特性，因此，为验证非线性反算法及其反演结果的可靠性与有效性，本研究同时分析了非线性反演结果并在反演结果下模拟了测井响应的误差。

图8展示了非线性反演结果下的模拟测井响应（虚线）与环境校正后的测井响应（黑色实线），涉及的测井项目依次为自然伽马GR、铀Uran、钍Th、中子孔隙度Nphi、体积密度DEN、宏观截面指数U、声波时差DT、冲洗带电导率CXO与原状地层电导率CT。从两组测井响应的对比看，非线性反演结果下的模拟测井响应与实测测井响应具有良好的一致性。表1中定量评分析了两组测井响应间的相关系数，各项测井响应的相关系数在0.867～0.996之间，相关系数均值达到0.921，充分反映了反演结果下的岩石组分宏观物理特性与真实储层物理特性的相似性，即说明通过非线性反演得到的岩石组分及其含量已十分接近页岩气储层的实际情况。此外，以实验室分析结果为标准，表2分别统计分析了图8中初始评价结果与非线性反演结果对实验室结果的相关系数，两组相关系数的对比可以说明，本研究建立的非线性反算法明显提高了页岩各岩石组分评价的精确度。因此，上述两方面分析充分证明，本研究建立的非线性反算法在解决页岩储层复杂岩石组分与含量评价问题方面的可靠性与有效性。

该方法能够同时解决页岩气储层岩石矿物组分与有机碳含量评价的两大问题，这两项问题的顺利解决对于后续储层脆性、吸附气含量等重要储层参数评价提供了科学的依据与技术保障。

图7 页岩岩石组分初始评价结果与非线性反演计算结果对比

表1 模拟测井响应与实测响应相关系数

图8 建南地区某井东岳庙段泥页岩复杂岩石组分反演质量控制

表2 初始评价及非线性反演评价较岩心分析结果的相关性对比

3 结论

本研究以非线性反演与最优化算法为核心思想建立的页岩气储层岩石组分测井评价方法，在鄂西渝东建南构造东岳庙段的页岩气储层评价中取得了较好的评价效果。该方法充分挖掘了常规测井资料中蕴含的丰富地质信息，同时解决了页岩储层重要矿物与有机碳含量评价两大问题，弥补了岩心分析深度不连续、ECS测井代价高昂的弊端，且极大地提高了测井评价结果的精度，为后续储层脆性与含气性的综合评价提供了科学的依据与重要的技术保障。

参考文献

［1］洪有密.测井原理与综合解释.东营：石油大学出版社，1993：1 ～365.

［2］Pemper R R，Sommer A，Guo P，et al.A new pulsed neutron sonde for derivation of formation lithology and mineralogy.SPE Annual Technical Conference and Exhibition，San Antonio，Texas，USA，2006：1～13.

［3］Sondergeld C H，Newsham K E，Comisky J T，et al.Petrophysical considerations in evaluating and procing shale gas resources.SPE Unconventional Gas Conference，Pittsburgh，Pennsylvania，USA，2010：1～34.

［4］Heidari Z，Torres－Verd C，Preeg W E.Quantitative method for estimating total organic carbon and porosity，and for diagnosing mineral constituents from well logs in shale－gas formations.2011：1 ～15.

［5］雍世和.最优化测井解释.东营：石油大学出版社，1995：87～119.

［6］潘仁芳，伍媛，宋争.页岩气勘探的地球化学指标及测井分析方法初探.中国石油勘探.2009，3：6～28.

［7］朱华，姜文利，边瑞康，等.页岩气资源评价方法体系及其应用——以川西坳陷为例.天然气工业.2009，29（12）：130～134.

［8］Passey Q R，Creaney S，Kulla J B，et al.A practical model for organic richness from porosity and resistivity logs.AAPG Bulletin.1990，74：1777～1794.

‘捌’ 深度优先算法解决八数码问题

首先我们来想象一只老鼠，在一座不见天日的迷宫内，老鼠在入口处进去，要从出 口出来。那老鼠会怎么走？当然是这样的：老鼠如果遇到直路，就一直往前走，如果遇到分叉路口，就任意选 择其中的一个继续往下走，如果遇到死胡同，就退回到最近的一个分叉路口，选择另一条道路再走下去，如果 遇到了出口，老鼠的旅途就算结束了。深度优先搜索法的基本原则就是这样：按照某种条件往前试探搜索，如 果前进中遭到失败（正如老鼠遇到死胡同）则退回头另选通路继续搜索，直到找到条件的目标为止。
实现这一算法，我们要用到编程的另一大利器--递归。递归是一个很抽象的概念， 但是在日常生活中，我们还是能够看到的。拿两面镜子来，把他们面对着面，你会看到什么？你会看到镜子中 有无数个镜子？怎么回事？A镜子中有B镜子的象，B镜子中有A镜子的象，A镜子的象就是A镜子本身的真实写 照，也就是说A镜子的象包括了A镜子，还有B镜子在A镜子中的象………………好累啊，又烦又绕口，还不好理 解。换成计算机语言就是A调用B，而B又调用A，这样间接的，A就调用了A本身，这实现了一个重复的功能。再 举一个例子；从前有座山，山里有座庙，庙里有个老和尚，老和尚在讲故事，讲什么呢？讲：从前有座山，山 里有座庙，庙里有个老和尚，老和尚在讲故事，讲什么呢？讲：从前有座山，山里有座庙，庙里有个老和尚， 老和尚在讲故事，讲什么呢？讲：…………。好家伙，这样讲到世界末日还讲不玩，老和尚讲的故事实际上就 是前面的故事情节，这样不断地调用程序本身，就形成了递归。 万一这个故事中的某一个老和尚看这个故事不 顺眼，就把他要讲的故事换成：“你有完没完啊！”，这样，整个故事也就嘎然而止了。我们编程就要注意这 一点，在适当的时候，就必须要有一个这样的和尚挺身而出，把整个故事给停下来，或者使他不再往深一层次 搜索，要不，我们的递归就会因计算机存储容量的限制而被迫溢出，切记，切记。
我们把递归思想运用到上面的迷宫中，记老鼠现在所在的位置是(x,y),那它现在有 前后左右4个方向可以走，分别是(x+1,y),(x-1,y),(x,y+1),(x,y-1),其中一个方向是它来时的路，我们先不 考虑，我们就分别尝试其他三个方向，如果某个方向是路而不是墙的话，老鼠就向那个方向迈出一步。在新的 位置上，我们又可以重复前面的步骤。老鼠走到了死胡同又是怎么回事？就是除了来时的路，其他3个方向都是 墙，这时这条路就走到了尽头，无法再向深一层发展，我们就应该沿来时的路回去，尝试另外的方向。
例：八皇后问题：在标准国际象棋的棋盘上（8*8格）准备放置8只皇后，我们知 道，国际象棋中皇后的威力是最大的，她既可以横走竖走，还可以斜着走，遇到挡在她前进路线上的敌人，她 就可以吃掉对手。要求在棋盘上安放8只皇后，使她们彼此互相都不能吃到对方，求皇后的放法。
这是一道很经典的题目了，我们先要明确一下思路，如何运用深度优先搜索法，完 成这道题目。我们先建立一个8*8格的棋盘，在棋盘的第一行的任意位置安放一只皇后。紧接着，我们就来放 第二行，第二行的安放就要受一些限制了，因为与第一行的皇后在同一竖行或同一对角线的位置上是不能安放 皇后的，接下来是第三行，……，或许我们会遇到这种情况，在摆到某一行的时候，无论皇后摆放在什么位 置，她都会被其他行的皇后吃掉，这说明什么呢？这说明，我们前面的摆放是失败的，也就是说，按照前面 的皇后的摆放方法，我们不可能得到正确的解。那这时怎么办？改啊，我们回到上一行，把原先我们摆好的 皇后换另外一个位置，接着再回过头摆这一行，如果这样还不行或者上一行的皇后只有一个位置可放，那怎 么办？我们回到上一行的上一行，这和老鼠碰了壁就回头是一个意思。就这样的不断的尝试