电位法在油田剩余油监测中的应用研究
|
摘要
对油气田大规模及深入的开发，已经导致整装的、类型相对简单的油气藏几近殆尽，实现油气增产上主要靠寻找新的复杂油气藏和提高老油田油气采收率，九十年代以来，剩余油问题和水驱方向及前缘问题是油气田开发工作中的一个研究热点，也是有关国际会议的重要议题，预测剩余油分布的新方法和新技术亟待研发。文章对几种常用的剩余油监测方法：地震技术、示踪剂监测方法、测井储层描述方法做了简要的介绍，阐述了各自的优势和不足。与其他方法比，电位法监测剩余油具有成本低、效率高、可实施动态监测的优点。该方法是以充电法的原理和地下动态导体的基本理论为依据，通过研究地下注入工作液前后的地表电位的变化，对剩余油的分布及工作液的流动情况做出定量及定性的判断，即在注入工作液的过程中，注入的工作液会减小注入目的层的电阻率，与围岩介质的电阻率有比较大的差异。
在该方法的数值模拟中，把总电位分成了正常电位和异常电位，以点电流源电位公式为基础得到了正常电位，用有限差分法得出来异常电位的方程组，利用不完全Cholesky共轭梯度方法求解大型稀疏线性方程组，计算得出了异常电位的表达式，消除了电源附近的奇异性，而且介绍了总电位、正常电位、异常电位的边界条件问题。在正演数值计算方法的基础上，首先将反演问题线性化，然后用光滑约束最小二乘法进行了反演计算，其中雅可比矩阵的计算，从有限差分得出的异常电位公式出发，运用了Tripp提出的灵敏度定律得出雅可比矩阵，避免了多次直接求解该矩阵。运用上面的算法，模拟地表供电、井中供电以及井地供电三种方式，根据正演和反演的模型，得出供电偶极子为井地供电方式时，对模型中设置的目标体探测能力较强，而其它两种供电模式对目标体没有分辨率。
为了验证电位法的正确，我们在中原油田做了实验。中原油田的某些区块已经进入高含水开发期，储层的物性、压力场和油水关系等发生了变化，剩余油的分布情况不清楚，并且油田的开发难度越来越大，针对这些难点，我们应用了电位法监测技术，测网布置成矩形，根据测试数据绘制了视电阻率等值线图，并且与该测区的饱和度等值线图做了对比，证实了电位法监测剩余油的可行性。除此之外，我们在新疆的风城油田也进行了试验，采用环形测网，利用电位监测结果进行对比，判断蒸汽的扩散方向和运移规律，得到了满意的效果。
The oil and gas reservoir are having sharp reduction, which are exploited easilyand of relatively simple type, with large-scale and in-depth development of oil andgas fields. Increasing production of oil and gas depends mainly on the search for newcomplex reservoir and improve the oil recovery of the old oil fields. Since the1990s,the remaining oil and water flooding on the direction and the leading edge are a hotresearch topic in the oil and gas field development work, they are also about theimportant issues of the International Conference. The new methods and newtechnologies of predicting distribution of remaining oil are in urgent need of researchand development. In the article, there are also brief introductions on several othercommonly used residual oil monitoring methods: seismic techniques, tracermonitoring methods， well logging reservoir description method，and elaboratingtheir respective strengths and weaknesses. Compared with those methods, thepotential method has lower cost, higher efficiency, and dynamic monitoring could beimplemented in this method. The method is based on the law principles ofmise-a-la-masse method and the basic theory of the underground dynamicconductor. To make quantitative and qualitative judgments on the distribution ofresidual oil and the flow of working fluid through studying the change of surfacepotential after injecting the working fluid. In other words, the injected working fluidreduce the resistivity of target layer during injecting working fluid, that leads torelatively large differences to the resistivity of surrounding rock medium.
In the method of numerical simulation, the total potential is divided into thenormal potential and abnormal potential, we got the normal potential based onpotential formula of the point current source, and the equations of the abnormalpotential using finite difference method. We applied the incomplete Choleskyconjugate gradient method to solve large sparse linear equations and got theexpression of the abnormal potential to eliminate the singularity near the power source. There are discussions on the boundary conditions of total potential, normalpotential and abnormal potential. Based on forward numerical calculation method，inverse problem was linearized, and then applying the smoothness constraint leastsquares method for the calculation of inversion problem, However, by applying thereciprocity principle, the partial derivatives can be very efficiently obtained withouthaving to directly solve the matrix equations a number of times. The basic idea isdescribed in Tripp et al.（1984） and is known as the sensitivity theorem. Accordingto the above algorithm, we simulated three power supplying ways: surface-surface,borehole-borehole, surface-borehole. From the forward and inversion model,weknow that electric potential method has high detection ability on the target body inmodel when surface-borehole power supplying way was used, no detection abilitywhile the other two power supplying ways.
There is an experiment in the Zhongyuan Oilfield to verify that potentialmethod’s result is correct. Some blocks of the Zhongyuan Oilfield have entered thehigh water cut stage, the reservoir’s physical properties, pressure field and oil-waterproportion have been changing, and it is not clear about the distribution of residualoil, and exploiting the oil fields more and more difficult. Because of those difficulties，we apply the potential method in these blocks, we got the apparent resistivitycontour maps drawn according to the test data, and they are compared with thesaturation of contour map, we got this conclusion that the potential method can beused for monitoring residual oil. In addition, we also carried out test in the WindyCity Oilfield in Xinjiang using ring measured net. Potential monitoring results werecompared to determine the direction of the proliferation and migration law of thesteam, and we get satisfactory results.
在该方法的数值模拟中，把总电位分成了正常电位和异常电位，以点电流源电位公式为基础得到了正常电位，用有限差分法得出来异常电位的方程组，利用不完全Cholesky共轭梯度方法求解大型稀疏线性方程组，计算得出了异常电位的表达式，消除了电源附近的奇异性，而且介绍了总电位、正常电位、异常电位的边界条件问题。在正演数值计算方法的基础上，首先将反演问题线性化，然后用光滑约束最小二乘法进行了反演计算，其中雅可比矩阵的计算，从有限差分得出的异常电位公式出发，运用了Tripp提出的灵敏度定律得出雅可比矩阵，避免了多次直接求解该矩阵。运用上面的算法，模拟地表供电、井中供电以及井地供电三种方式，根据正演和反演的模型，得出供电偶极子为井地供电方式时，对模型中设置的目标体探测能力较强，而其它两种供电模式对目标体没有分辨率。
为了验证电位法的正确，我们在中原油田做了实验。中原油田的某些区块已经进入高含水开发期，储层的物性、压力场和油水关系等发生了变化，剩余油的分布情况不清楚，并且油田的开发难度越来越大，针对这些难点，我们应用了电位法监测技术，测网布置成矩形，根据测试数据绘制了视电阻率等值线图，并且与该测区的饱和度等值线图做了对比，证实了电位法监测剩余油的可行性。除此之外，我们在新疆的风城油田也进行了试验，采用环形测网，利用电位监测结果进行对比，判断蒸汽的扩散方向和运移规律，得到了满意的效果。
The oil and gas reservoir are having sharp reduction, which are exploited easilyand of relatively simple type, with large-scale and in-depth development of oil andgas fields. Increasing production of oil and gas depends mainly on the search for newcomplex reservoir and improve the oil recovery of the old oil fields. Since the1990s,the remaining oil and water flooding on the direction and the leading edge are a hotresearch topic in the oil and gas field development work, they are also about theimportant issues of the International Conference. The new methods and newtechnologies of predicting distribution of remaining oil are in urgent need of researchand development. In the article, there are also brief introductions on several othercommonly used residual oil monitoring methods: seismic techniques, tracermonitoring methods， well logging reservoir description method，and elaboratingtheir respective strengths and weaknesses. Compared with those methods, thepotential method has lower cost, higher efficiency, and dynamic monitoring could beimplemented in this method. The method is based on the law principles ofmise-a-la-masse method and the basic theory of the underground dynamicconductor. To make quantitative and qualitative judgments on the distribution ofresidual oil and the flow of working fluid through studying the change of surfacepotential after injecting the working fluid. In other words, the injected working fluidreduce the resistivity of target layer during injecting working fluid, that leads torelatively large differences to the resistivity of surrounding rock medium.
In the method of numerical simulation, the total potential is divided into thenormal potential and abnormal potential, we got the normal potential based onpotential formula of the point current source, and the equations of the abnormalpotential using finite difference method. We applied the incomplete Choleskyconjugate gradient method to solve large sparse linear equations and got theexpression of the abnormal potential to eliminate the singularity near the power source. There are discussions on the boundary conditions of total potential, normalpotential and abnormal potential. Based on forward numerical calculation method，inverse problem was linearized, and then applying the smoothness constraint leastsquares method for the calculation of inversion problem, However, by applying thereciprocity principle, the partial derivatives can be very efficiently obtained withouthaving to directly solve the matrix equations a number of times. The basic idea isdescribed in Tripp et al.（1984） and is known as the sensitivity theorem. Accordingto the above algorithm, we simulated three power supplying ways: surface-surface,borehole-borehole, surface-borehole. From the forward and inversion model,weknow that electric potential method has high detection ability on the target body inmodel when surface-borehole power supplying way was used, no detection abilitywhile the other two power supplying ways.
There is an experiment in the Zhongyuan Oilfield to verify that potentialmethod’s result is correct. Some blocks of the Zhongyuan Oilfield have entered thehigh water cut stage, the reservoir’s physical properties, pressure field and oil-waterproportion have been changing, and it is not clear about the distribution of residualoil, and exploiting the oil fields more and more difficult. Because of those difficulties，we apply the potential method in these blocks, we got the apparent resistivitycontour maps drawn according to the test data, and they are compared with thesaturation of contour map, we got this conclusion that the potential method can beused for monitoring residual oil. In addition, we also carried out test in the WindyCity Oilfield in Xinjiang using ring measured net. Potential monitoring results werecompared to determine the direction of the proliferation and migration law of thesteam, and we get satisfactory results.
引文
[1]徐世浙.地球物理中的有限单元法[M].北京：科学出版社，1994.
[2]傅良魁，李金铭，何继善，等.电法勘探教程[M].北京：地质出版社，1983.
[3]周熙襄，钟本善.电法勘探数值模拟技术[M].成都：四川科学技术出版社，1986：163-239.
[4]数学手册编写组.数学手册[M].北京：高等教育出版社，1979：263.
[5]姚姚.地球物理反演基本理论与应用方法[M].武汉：中国地质大学出版社，2002
[6]印兴耀，韩文功，李振春，等.地震技术新进展[M].北京：中国石油大学出版社，2006
[7]杨文采．地球物理反演的理论与方法[M].北京：地质出版社，1997．
[8]王家映．地球物理反演理论[M].北京：高等教育出版社，2002.
[9]廉洁.垂直线源井地直流电法正反演研究[D].北京：中国地质大学，2007.
[10] Loke,M.H.,Dahlin,T.A combined Gauss–Newton and quasi-Newton inversionmethod for the interpretation of apparent resistivitypseudosections.Sept8–11,1997
[C]．Aarhus， Denmark: Procs.3rd Meeting of the European Association forEnvironmental and Engineering Geophysics
[11]黄俊革.三维电阻率/极化率有限元正演模拟与反演成像[D].湖南：中南大学，2003.
[12]刘海飞.直流激电反演中的线性与非线性方法研究[D].湖南：中南大学，2007.
[13]陈德鹏，任意线源井地电法三维有限元正演研究[D].湖南：中南大学，2009.
[14]章钰，三维井-地电位测井技术的正演模拟研究[D].山东：中国海洋大学，2009.
[15]吴小平.利用共扼梯度方法的电阻率三维正、反演研究[D].合肥：中国科学技术大学，1998.
[16]王志刚.油田开发直流充电法三维模拟研究[D].北京：中国地质大学，2004.
[17]李宗田．储层三维电阻率成像技术在剩余油开发中的应用[D].北京：中国地质大学，2005．
[18]何梅兴.可控源音频大地电磁二维OCCAM反演研究[D].北京：中国地质大学.2009.
[19]牛建军，张晓培.油田注水监测中的电位测量技术[J].吉林大学学报（地球科学版），2002，32（增刊）：1-5.
[20]牛建军，张晓培，杜立志.稠油汽驱电位法动态监测技术[J].吉林大学学报（地球科学版），2008，38（5）：887-890.
[21]张世洪，张晓培，牛建军，等. SE2432型高密度电法并行采集仪及采集工作方法[J].吉林大学学报（地球科学版），2002，32（增刊）：95-97.
[22]张晓培，张琪，吴刚，等.分布式高密度电法仪的设计[J].吉林大学学报（地球科学版），2002，32（增刊）：9-10.
[23]何裕盛.地下动态导体充电法探测概论[J].物探与化探，2000，24（1）：62-67.
[24]张金成.电位法井间监测技术[J].地震地质，2001，23（2）：292-300.
[25]胡忠太，甄宝生，周怀亮，等.电位法井间监测技术[J].油气井测试，2004，（4）：62-63.
[26]梁生，韩立国，李安帮.微地震监测技术在秦家屯地区的应用[J].吉林大学学报（地球科学版），2007，37（增刊）：95-98.
[27]郑秀芬，李金铭.油气电法勘探综述[J].地球物理学进展，1997，12（2）：89-96
[28]徐凯军，李桐林，张辉，等.基于共轭梯度法的垂直有限线源三维电阻率反演[J].煤田地质与勘探，2006，34（3）：69-71
[29]徐凯军，李桐林.垂直有限线源三维地电场有限差分正演研究[J].吉林大学学报（地球科学版），2006，36（1）：137-141.
[30]安然，李桐林，徐凯军.井地三维电阻率反演研究[J].地球物理学进展，2007，22（1）：247-249.
[31]刘显，土志刚，何展翔.井地电法供电电场分布模拟研究[J].工程地球物理学报，2006，3（5）：331-336.
[32]王志刚，何展翔，刘昱.井地直流电法三维数值模拟及异常规律研究[J].工程地球物理学报，2006，3（2）：88-92
[33]贾正元，晋风明.井-地电法勘探在油藏评价中的应用[J].工程地球物理学报，2008，5（3）：326-331.
[34]屈有恒，张贵宾，赵连锋，等.井地有限线源三维电阻率反演研究[J].地球物理学进展，2007，22（5）：1393-1402..
[35]凌明友，郝新武，程同军，等.井-地电位成像技术研究剩余油分布[J].断块油气藏，2003，10（4）：55-58
[36]吴小平，徐果明.利用ICCG迭代技术加快电阻率三维正演计算[J].煤田地质与勘探，1999，（3）：62-66
[37]徐世浙，刘斌，阮百尧.电阻率法中求解异常电位的有限单元法[J].地球物理学报，1994，37（2）：511-515.
[38]阮百尧，村上裕，徐世浙.电阻率/极化率数据的二维反演程序[J].物探化探计算技术，1999，21（2）：116-125.
[39]黄俊革，阮百尧，鲍光淑．有限单元法三维电阻率最小二乘反演中存在问题的研究［J］．地质与勘探，2004，40（4）：70-75．
[40]黄俊革，阮百尧，鲍光淑．基于有限单元法的三维地电断面电阻率反演［J］．中南大学学报（自然科学版），2004，35（2）：295-299．
[41]黄俊革，阮百尧，鲍光淑．有限单元法三维电阻率最小二乘反演中存在问题的研究［J］．地质与勘探，2004，40（4）：70-75.
[42]黄俊革，阮百尧．直流电阻率测深中二维与三维反演结果的对比与分析［J］．物探与化探，2004，28（5）：447-450.
[43]阮百尧．视电阻率对模型电阻率的偏导数矩阵计算方法［J］．地质与勘探，2001，37（6）：39-41.
[44]阮百尧，熊彬，徐世浙.三维地电断面电阻率测深有限元数值模拟[J].地球科学———中国地质大学学报，2001，26（1）：73-77.
[45]刘海飞，陈德鹏，戴前伟，柳建新，阮百尧.连续电性介质线源井-地电位三维有限元数值模拟[J].桂林理工大学学报，2011，31（1）：28-37.
[46]吴小平，汪彤彤.利用共扼梯度算法的电阻率三维有限元正演[J].地球物理学报，2003，46（3）：428-432.
[47]吴小平，徐果明，李时灿.利用不完全Cho lesky共轭梯度法求解点源三维电场[J].地球物理学报，1998，41（6）：848-855.
[48]吴小平，徐果明.不完全Cholesky共扼梯度法及其在地电场计算中的应用[J].石油地球物理勘探，1998，33（1）：89-94.
[49]昊小平，徐果明，李时灿.解大型稀疏方程组的ICCG方法及其计算机实现[J].煤田地质与勘探，1999，27（6）：54-55
[50]吴小平，徐果明，李时灿.利用不完全Cholesky共扼梯度法求解点源三维地电场[J].地球物理学报，1998，41（6）：848-855.
[51]吴小平，徐果明.利用共扼梯度方法电阻率三维反演研究[J].地球物理学报，2000，43（3）：420-427
[52]吴小平，徐果明.电阻率三维反演中偏导数矩阵的求取于分析[J].石油地球物理勘探，1999，34（4）：363－372
[53]吴小平，汪彤彤．电阻率三维反演方法研究进展［J］．地球物理学进展．2002，17（1）：156-62．
[54]吴小平，徐果明．电阻率三维反演中偏导数矩阵的求取与分析［J］．石油地球物理勘探，1998，34（4）：363-372．
[55]吴小平，汪彤彤.电阻率三维复杂结构的快速反演[J].煤田地质与勘探，2001，29（2）：48-51.
[56]潭河清，沈金松，张福莱.井-地电位成像技术及其在孤东八区剩余油分布研究中的应用[J].石油大学学报（自然科学版），2004，28（2）：31-37.
[57]王贤君，李建阁，冯立.应用井-地电位测量技术判断油层高渗透带方向[J].中国石油大学学报：自然科学版，2006，30（5）：67.
[58]王志刚，何展翔，魏文博.井地电法三维物理模型实验[J]..石油地球物理勘探，2005.40（5）;594-597.
[59]邓正栋，关洪军，万乐，等.稳定点电流场三维有限差分正演模拟[J].解放军理工大学学报，2000，1（3）：45-50.
[60]刘地渊，徐凯军，赵广茂，等.任意形状线电流源三维地电场研究[J].地球物理学进展，2006，21（2）：395-399.
[61]屈有恒，张贵宾，晋风明.倾斜线源的三维电场数值模拟研究[J].物探化探计算技术，2007，29（5）：431-435.
[62]水永秀树，牛岛惠辅.流电电位法三维模拟[J].物理探查，1991，44（4）：215-226.
[63]邓正栋，吴洪军，聂永平，等.稳定地电场三维有限差分正演模拟[J].石油物探，2001，（1）：107-114
[64]甘利灯，姚逢昌，杜文辉，等.水驱油藏四维地震技术[J].石油勘探与开发，2007，37（4）：437-444
[65]宛新林，席道瑛，高尔根，等．用改进的光滑约束最小二乘正交分解法实现电阻率三维反演［J］.地球物理学报，2005，48（1）：439-444．
[66]张大海，王兴泰.二维视电阻率断面的快速最小二乘反演［J］．物探化探计算技术，1999，21（1）：2-8．
[67]刘斌，李术才，李树忱，等．基于不等式约束的最小二乘法三维电阻率反演及其算法优化［J］．地球物理学报，2012，55（1）：260-268.
[68]李天成，牛滨华，孙春岩，慕洪涛.跨孔电阻率法有限元数值模拟及反演［J］．地质与勘探，2008，44（3）：52-55．
[69]刘申芬，李天成，慕洪涛.三维井地电阻率有限元数值模拟及反演［J］．物探与化探，2009，33（1）：88-90．
[70]崔永谦，刘池洋，张以明.油藏动态监测技术：时延（四维）地震述评［J］．石油与天然气地质，2004，25（1）：81-84
[71]何琰，吴念胜.三维地震在油藏挖潜中的应用［J］．西南石油学院学报，1999，21（2）：16-19
[72]韩德辉，智卫东，张丽娜.三维地震解释技术在储层预测中的应用研究［J］．江汉石油学院学报，2001，23（增刊）：41-42
[73] Wiggins, Ralph..The general linear inverse problem： Implications of surfacewaves and free oscillations on earth structure [J].Rev. Geoplzys. SpacePhys,1972,10，in press.
[74] JACKSOND，. D.Interpretation of inaccurate， insufficient and inconsistent data[J]. Geophysical Journal of the Royal Astronomical Society,1972,28:97-110.
[75] ORISTAGLIMO， L. and WORTHINGTOMN.，H. Inversion of surface andborehole electromagnetic data for two-dimensional electrical conductivity models[J].Geophysical Prospecting,1980,28:633-657.
[76] WIGGINSR, A..The general linear inverse problem: Implication of surface wavesand free oscillations for Earth structure [J].Reviews of Geophysics and SpacePhysics,1972,10:251-285.
[77] deGroot-Hedlin,C. and Constable, S. Occam’s inversion to generatesmooth,two-dimensional models from magneto telluric data[J]. Geophysics,1990,55:1613–1624.
[78] Marquardt,D. W. An algorithm for least-squares estimation of nonlinearparameters [J]. J. Soc. Indust. Appl. Math,1963,11:431–441.
[79] McGillivray, P. R, and Oldenburg, D. W. Methods for calculating Freshetderivatives and sensitivities for the non-linear inverse problem [J]. A comparativestudy， Geophysics. Props.，1990,38:499–524.
[80] Sasaki,Y.. Two-dimensional joint inversion of magneto telluric and dipole-dipoleresistivity data [J]. Geophysics,1989,54:254-262.
[81] Loke M.H. and Barker R.D. Rapid least-squares inversion of apparent resistivitypseudo sections using a quasi-Newton method [J]. Geophysical Prospecting,1996.44:l3l-152.
[82] Sasaki Y. Two-dimensional joint inversion of magneto telluric and dipole-dipoleresistivity data [J].Geophysics,1989,54:254-262.
[83] Sasaki Y.3-D resistivity inversion using the finite-element method [J].Geophysics,1994.59:1839-1848.
[84] Constable,S. C, Parker, R. L., and Constable, C G.. Occam's inversion： Apractical algorithm for generating smooth models from electromagnetic soundingdata [J]. Geophysics，1987,52:289-300.
[85] Tripp,A. C, Hohmann,G. W.,and Swift,eM., Jr.. Two dimensional resistivityinversions [J]. Geophysics,1984,49:1708-1717.
[86] McGillivray, P.R.,Oldenburg,D.W... Methods for calculating Freshet derivativesand sensitivities for the non-linear inverse problem [J]. a comparative study.Geophysical Prospecting,1990,38:499–524.
[87] Rodi, W.， Mackie,R.L.. Nonlinear conjugate gradients algorithm for2-Dmagneto telluric inversion [J]. Geophysics,2001,66:174–187.
[88] Loke,M.H.， Dahlin, T.A comparison of the Gauss–Newton and quasi-Newtonmethods in resistivity imaging inversion [J]. Journal of AppliedGeophysics,2002,49:149–162
[89] Oldenburg， D.W.， Li， Y.，. Inversion of induced polarization data [J].Geophysics,1994,59:1327–1341.
[2]傅良魁，李金铭，何继善，等.电法勘探教程[M].北京：地质出版社，1983.
[3]周熙襄，钟本善.电法勘探数值模拟技术[M].成都：四川科学技术出版社，1986：163-239.
[4]数学手册编写组.数学手册[M].北京：高等教育出版社，1979：263.
[5]姚姚.地球物理反演基本理论与应用方法[M].武汉：中国地质大学出版社，2002
[6]印兴耀，韩文功，李振春，等.地震技术新进展[M].北京：中国石油大学出版社，2006
[7]杨文采．地球物理反演的理论与方法[M].北京：地质出版社，1997．
[8]王家映．地球物理反演理论[M].北京：高等教育出版社，2002.
[9]廉洁.垂直线源井地直流电法正反演研究[D].北京：中国地质大学，2007.
[10] Loke,M.H.,Dahlin,T.A combined Gauss–Newton and quasi-Newton inversionmethod for the interpretation of apparent resistivitypseudosections.Sept8–11,1997
[C]．Aarhus， Denmark: Procs.3rd Meeting of the European Association forEnvironmental and Engineering Geophysics
[11]黄俊革.三维电阻率/极化率有限元正演模拟与反演成像[D].湖南：中南大学，2003.
[12]刘海飞.直流激电反演中的线性与非线性方法研究[D].湖南：中南大学，2007.
[13]陈德鹏，任意线源井地电法三维有限元正演研究[D].湖南：中南大学，2009.
[14]章钰，三维井-地电位测井技术的正演模拟研究[D].山东：中国海洋大学，2009.
[15]吴小平.利用共扼梯度方法的电阻率三维正、反演研究[D].合肥：中国科学技术大学，1998.
[16]王志刚.油田开发直流充电法三维模拟研究[D].北京：中国地质大学，2004.
[17]李宗田．储层三维电阻率成像技术在剩余油开发中的应用[D].北京：中国地质大学，2005．
[18]何梅兴.可控源音频大地电磁二维OCCAM反演研究[D].北京：中国地质大学.2009.
[19]牛建军，张晓培.油田注水监测中的电位测量技术[J].吉林大学学报（地球科学版），2002，32（增刊）：1-5.
[20]牛建军，张晓培，杜立志.稠油汽驱电位法动态监测技术[J].吉林大学学报（地球科学版），2008，38（5）：887-890.
[21]张世洪，张晓培，牛建军，等. SE2432型高密度电法并行采集仪及采集工作方法[J].吉林大学学报（地球科学版），2002，32（增刊）：95-97.
[22]张晓培，张琪，吴刚，等.分布式高密度电法仪的设计[J].吉林大学学报（地球科学版），2002，32（增刊）：9-10.
[23]何裕盛.地下动态导体充电法探测概论[J].物探与化探，2000，24（1）：62-67.
[24]张金成.电位法井间监测技术[J].地震地质，2001，23（2）：292-300.
[25]胡忠太，甄宝生，周怀亮，等.电位法井间监测技术[J].油气井测试，2004，（4）：62-63.
[26]梁生，韩立国，李安帮.微地震监测技术在秦家屯地区的应用[J].吉林大学学报（地球科学版），2007，37（增刊）：95-98.
[27]郑秀芬，李金铭.油气电法勘探综述[J].地球物理学进展，1997，12（2）：89-96
[28]徐凯军，李桐林，张辉，等.基于共轭梯度法的垂直有限线源三维电阻率反演[J].煤田地质与勘探，2006，34（3）：69-71
[29]徐凯军，李桐林.垂直有限线源三维地电场有限差分正演研究[J].吉林大学学报（地球科学版），2006，36（1）：137-141.
[30]安然，李桐林，徐凯军.井地三维电阻率反演研究[J].地球物理学进展，2007，22（1）：247-249.
[31]刘显，土志刚，何展翔.井地电法供电电场分布模拟研究[J].工程地球物理学报，2006，3（5）：331-336.
[32]王志刚，何展翔，刘昱.井地直流电法三维数值模拟及异常规律研究[J].工程地球物理学报，2006，3（2）：88-92
[33]贾正元，晋风明.井-地电法勘探在油藏评价中的应用[J].工程地球物理学报，2008，5（3）：326-331.
[34]屈有恒，张贵宾，赵连锋，等.井地有限线源三维电阻率反演研究[J].地球物理学进展，2007，22（5）：1393-1402..
[35]凌明友，郝新武，程同军，等.井-地电位成像技术研究剩余油分布[J].断块油气藏，2003，10（4）：55-58
[36]吴小平，徐果明.利用ICCG迭代技术加快电阻率三维正演计算[J].煤田地质与勘探，1999，（3）：62-66
[37]徐世浙，刘斌，阮百尧.电阻率法中求解异常电位的有限单元法[J].地球物理学报，1994，37（2）：511-515.
[38]阮百尧，村上裕，徐世浙.电阻率/极化率数据的二维反演程序[J].物探化探计算技术，1999，21（2）：116-125.
[39]黄俊革，阮百尧，鲍光淑．有限单元法三维电阻率最小二乘反演中存在问题的研究［J］．地质与勘探，2004，40（4）：70-75．
[40]黄俊革，阮百尧，鲍光淑．基于有限单元法的三维地电断面电阻率反演［J］．中南大学学报（自然科学版），2004，35（2）：295-299．
[41]黄俊革，阮百尧，鲍光淑．有限单元法三维电阻率最小二乘反演中存在问题的研究［J］．地质与勘探，2004，40（4）：70-75.
[42]黄俊革，阮百尧．直流电阻率测深中二维与三维反演结果的对比与分析［J］．物探与化探，2004，28（5）：447-450.
[43]阮百尧．视电阻率对模型电阻率的偏导数矩阵计算方法［J］．地质与勘探，2001，37（6）：39-41.
[44]阮百尧，熊彬，徐世浙.三维地电断面电阻率测深有限元数值模拟[J].地球科学———中国地质大学学报，2001，26（1）：73-77.
[45]刘海飞，陈德鹏，戴前伟，柳建新，阮百尧.连续电性介质线源井-地电位三维有限元数值模拟[J].桂林理工大学学报，2011，31（1）：28-37.
[46]吴小平，汪彤彤.利用共扼梯度算法的电阻率三维有限元正演[J].地球物理学报，2003，46（3）：428-432.
[47]吴小平，徐果明，李时灿.利用不完全Cho lesky共轭梯度法求解点源三维电场[J].地球物理学报，1998，41（6）：848-855.
[48]吴小平，徐果明.不完全Cholesky共扼梯度法及其在地电场计算中的应用[J].石油地球物理勘探，1998，33（1）：89-94.
[49]昊小平，徐果明，李时灿.解大型稀疏方程组的ICCG方法及其计算机实现[J].煤田地质与勘探，1999，27（6）：54-55
[50]吴小平，徐果明，李时灿.利用不完全Cholesky共扼梯度法求解点源三维地电场[J].地球物理学报，1998，41（6）：848-855.
[51]吴小平，徐果明.利用共扼梯度方法电阻率三维反演研究[J].地球物理学报，2000，43（3）：420-427
[52]吴小平，徐果明.电阻率三维反演中偏导数矩阵的求取于分析[J].石油地球物理勘探，1999，34（4）：363－372
[53]吴小平，汪彤彤．电阻率三维反演方法研究进展［J］．地球物理学进展．2002，17（1）：156-62．
[54]吴小平，徐果明．电阻率三维反演中偏导数矩阵的求取与分析［J］．石油地球物理勘探，1998，34（4）：363-372．
[55]吴小平，汪彤彤.电阻率三维复杂结构的快速反演[J].煤田地质与勘探，2001，29（2）：48-51.
[56]潭河清，沈金松，张福莱.井-地电位成像技术及其在孤东八区剩余油分布研究中的应用[J].石油大学学报（自然科学版），2004，28（2）：31-37.
[57]王贤君，李建阁，冯立.应用井-地电位测量技术判断油层高渗透带方向[J].中国石油大学学报：自然科学版，2006，30（5）：67.
[58]王志刚，何展翔，魏文博.井地电法三维物理模型实验[J]..石油地球物理勘探，2005.40（5）;594-597.
[59]邓正栋，关洪军，万乐，等.稳定点电流场三维有限差分正演模拟[J].解放军理工大学学报，2000，1（3）：45-50.
[60]刘地渊，徐凯军，赵广茂，等.任意形状线电流源三维地电场研究[J].地球物理学进展，2006，21（2）：395-399.
[61]屈有恒，张贵宾，晋风明.倾斜线源的三维电场数值模拟研究[J].物探化探计算技术，2007，29（5）：431-435.
[62]水永秀树，牛岛惠辅.流电电位法三维模拟[J].物理探查，1991，44（4）：215-226.
[63]邓正栋，吴洪军，聂永平，等.稳定地电场三维有限差分正演模拟[J].石油物探，2001，（1）：107-114
[64]甘利灯，姚逢昌，杜文辉，等.水驱油藏四维地震技术[J].石油勘探与开发，2007，37（4）：437-444
[65]宛新林，席道瑛，高尔根，等．用改进的光滑约束最小二乘正交分解法实现电阻率三维反演［J］.地球物理学报，2005，48（1）：439-444．
[66]张大海，王兴泰.二维视电阻率断面的快速最小二乘反演［J］．物探化探计算技术，1999，21（1）：2-8．
[67]刘斌，李术才，李树忱，等．基于不等式约束的最小二乘法三维电阻率反演及其算法优化［J］．地球物理学报，2012，55（1）：260-268.
[68]李天成，牛滨华，孙春岩，慕洪涛.跨孔电阻率法有限元数值模拟及反演［J］．地质与勘探，2008，44（3）：52-55．
[69]刘申芬，李天成，慕洪涛.三维井地电阻率有限元数值模拟及反演［J］．物探与化探，2009，33（1）：88-90．
[70]崔永谦，刘池洋，张以明.油藏动态监测技术：时延（四维）地震述评［J］．石油与天然气地质，2004，25（1）：81-84
[71]何琰，吴念胜.三维地震在油藏挖潜中的应用［J］．西南石油学院学报，1999，21（2）：16-19
[72]韩德辉，智卫东，张丽娜.三维地震解释技术在储层预测中的应用研究［J］．江汉石油学院学报，2001，23（增刊）：41-42
[73] Wiggins, Ralph..The general linear inverse problem： Implications of surfacewaves and free oscillations on earth structure [J].Rev. Geoplzys. SpacePhys,1972,10，in press.
[74] JACKSOND，. D.Interpretation of inaccurate， insufficient and inconsistent data[J]. Geophysical Journal of the Royal Astronomical Society,1972,28:97-110.
[75] ORISTAGLIMO， L. and WORTHINGTOMN.，H. Inversion of surface andborehole electromagnetic data for two-dimensional electrical conductivity models[J].Geophysical Prospecting,1980,28:633-657.
[76] WIGGINSR, A..The general linear inverse problem: Implication of surface wavesand free oscillations for Earth structure [J].Reviews of Geophysics and SpacePhysics,1972,10:251-285.
[77] deGroot-Hedlin,C. and Constable, S. Occam’s inversion to generatesmooth,two-dimensional models from magneto telluric data[J]. Geophysics,1990,55:1613–1624.
[78] Marquardt,D. W. An algorithm for least-squares estimation of nonlinearparameters [J]. J. Soc. Indust. Appl. Math,1963,11:431–441.
[79] McGillivray, P. R, and Oldenburg, D. W. Methods for calculating Freshetderivatives and sensitivities for the non-linear inverse problem [J]. A comparativestudy， Geophysics. Props.，1990,38:499–524.
[80] Sasaki,Y.. Two-dimensional joint inversion of magneto telluric and dipole-dipoleresistivity data [J]. Geophysics,1989,54:254-262.
[81] Loke M.H. and Barker R.D. Rapid least-squares inversion of apparent resistivitypseudo sections using a quasi-Newton method [J]. Geophysical Prospecting,1996.44:l3l-152.
[82] Sasaki Y. Two-dimensional joint inversion of magneto telluric and dipole-dipoleresistivity data [J].Geophysics,1989,54:254-262.
[83] Sasaki Y.3-D resistivity inversion using the finite-element method [J].Geophysics,1994.59:1839-1848.
[84] Constable,S. C, Parker, R. L., and Constable, C G.. Occam's inversion： Apractical algorithm for generating smooth models from electromagnetic soundingdata [J]. Geophysics，1987,52:289-300.
[85] Tripp,A. C, Hohmann,G. W.,and Swift,eM., Jr.. Two dimensional resistivityinversions [J]. Geophysics,1984,49:1708-1717.
[86] McGillivray, P.R.,Oldenburg,D.W... Methods for calculating Freshet derivativesand sensitivities for the non-linear inverse problem [J]. a comparative study.Geophysical Prospecting,1990,38:499–524.
[87] Rodi, W.， Mackie,R.L.. Nonlinear conjugate gradients algorithm for2-Dmagneto telluric inversion [J]. Geophysics,2001,66:174–187.
[88] Loke,M.H.， Dahlin, T.A comparison of the Gauss–Newton and quasi-Newtonmethods in resistivity imaging inversion [J]. Journal of AppliedGeophysics,2002,49:149–162
[89] Oldenburg， D.W.， Li， Y.，. Inversion of induced polarization data [J].Geophysics,1994,59:1327–1341.
© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号 地址：北京市海淀区学院路29号 邮编：100083 电话：办公室：(+86 10)66554848；文献借阅、咨询服务、科技查新：66554700 |