各地的磁倾角_试论磁倾角异常作为确定烃类渗漏带的参考系

参考系问题的提出非地震物化探探测油气田的原理,主要是根据烃类渗漏带构成“烟囱”模式或“氧化一还原灶”(Pirson,1981;Tompkins,199O),但有时也以“盖层堵塞作用”(Rosaire,1940)等模式来解释油气田上方的晕状异常、环状异常、顶端异常等。目前,非地震物探找油方法主要有自然电位、激发极化、测温、放射性、微磁、微重力、热释光等,化探方法主要有烃类测量、三角碳、汞气、示踪元素、荧光等,生物方法主要有细菌分析方法。众所周知,石油天然气运移过程可分为四种状态:其一为临界逸散状态,将向大气逸散;其二为迁移之中的有机质状态;其三为介质表面的吸附状态;其四是介质之中的固化状态。前三种多呈油气组分的相态,后者则是固化的组分,是在地史时期中经烃类一生物一化学动力学复合作用的结果,是已固化在岩石或土壤中的石油天然气作用的信息。非地震物化探探测油气田方法与这四种烃类组分的可能关系见图1。至今多种非地震找油方法,已作了很多的工作,也发现了很多异常,但在发现油气远景区以及油气储层方法,不如地震勘探的效果好。究其原因,就是对非地震物化探方法的异常认识不足,包括异常机理、异常形态、异常与油气储层的关系等。同时,在实际工作中,非地震物化探方法确定异常的规则也不尽相同,通常要由经验丰富的专家来研究异常、确定异常。由此看来,在非地震物化探找油工作之中,是否存在一种全球普遍适用的规则,作为确定烃类渗漏带异常的标准?因非地震物化探寻找油气田的方法都是利用近地表介质(一定深度的岩石、表层土壤)的特性,而有用的信息却是介质形成之后因烃类渗漏导致的各种物理的、化学的、生物的以及综合的次生作用,故作为确定烃类渗漏带的首要标准是:这种参考系应是介质形成之后形成的,与介质形成作用无关,而完全是由于烃类引起的各种次生作用所决定。至于采用哪些方法或哪些参数可以确定成为烃类渗漏带的参考系,在此不加以讨论,只对土壤磁性方法中的磁倾角作为全球均可适用的确定烃类渗漏的参考系加以初步讨论。磁倾角作为全球性确定烃类渗漏带的参考系土坯磁性异常与烃类作用的关系1979年Donovan等人在美国俄克拉何马州的Cement油气田上方发现了浅源高频微磁异常,自此,油气田与磁异常的关系为各国石油天然气的勘探者所关注。目前他们的研究结果认为,油气田上方磁异常可能起源于:烃类渗漏逸出的烃类物质与硫酸盐还原细菌的作用形成次生磁铁矿(Donovan,etal,1979;Ferguson,1979);烃类的生物降解作用形成次生磁铁矿(Elmore,etal,1987;McCabe,etal,1987);原油的微生物氧化和硫酸盐还原,形成次生磁黄铁矿(Sassen,etal,1988);烃类作用下,碳酸根浓度高时形成菱铁矿(Machel,etal,1987);油气渗漏带中的自燃变质作用形成的赤铁矿一钛铁矿一磁铁矿固溶体系列(Cisowski,etal,,1987);烃类导致陆相红层磁性下降(Kilgore,etal,1987);烃类、微生物、化学热力学的复合过程形成次生磁铁矿等矿物(Machel,etal,1991);烃类的化学吸附作用对剩磁稳定性的影响(孟小红,周海民,1997)。最近,郭友别等人(199)讨论了决定烃类渗漏带磁性异常的因素有四个:烃类物质的数量;硫酸盐还原细菌群及其生态;土壤中铁的含量及其价态;土壤接受烃类作用持续的时间等。这四个因素往往同时存在、共同决定烃类渗漏带的一切物理的、化学的、生物的反应,但它们在数量上或性质上的不同,决定了反应的发展方向和构成渗漏带磁性异常的差异。从上述数种土壤磁

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
圆柱体的磁异常可以通过计算圆柱体内部和外部的磁场差异得到。在Matlab中,可以使用数值计算工具箱中的函数来计算圆柱体的磁场分布和磁异常。以下是一个简单的示例代码: ```matlab % 圆柱体参数 R = 10; % 圆柱体半径 L = 50; % 圆柱体长度 chi = 0.01; % 圆柱体磁化率 inc = 30; % 磁倾角 dec = 0; % 磁偏角 % 地球磁场参数 B0 = 50000; % 地球磁场强度 inc0 = 60; % 地球磁场倾角 dec0 = 0; % 地球磁场偏角 % 计算磁场分布 [x, y, z] = meshgrid(-100:5:100, -100:5:100, -100:5:100); [Bx, By, Bz] = cylinder_mag(R, L, chi, inc, dec, x, y, z, B0, inc0, dec0); % 计算磁异常 [Bx0, By0, Bz0] = cylinder_mag(R, L, 0, inc, dec, x, y, z, B0, inc0, dec0); dBx = Bx - Bx0; dBy = By - By0; dBz = Bz - Bz0; dB = sqrt(dBx.^2 + dBy.^2 + dBz.^2); % 绘制磁异常图像 figure; slice(x, y, z, dB, [], [], [-50, 0, 50]); colorbar; xlabel('X'); ylabel('Y'); zlabel('Z'); title('Cylinder Magnetic Anomaly'); ``` 该代码使用了一个名为`cylinder_mag`的自定义函数,用于计算圆柱体的磁场分布和磁异常。以下是该函数的代码: ```matlab function [Bx, By, Bz] = cylinder_mag(R, L, chi, inc, dec, x, y, z, B0, inc0, dec0) % 计算圆柱体的磁场分布和磁异常 % R: 圆柱体半径 % L: 圆柱体长度 % chi: 圆柱体磁化率 % inc: 磁倾角 % dec: 磁偏角 % x, y, z: 空间坐标 % B0: 地球磁场强度 % inc0: 地球磁场倾角 % dec0: 地球磁场偏角 mu0 = 4*pi*1e-7; % 真空磁导率 % 圆柱体内部的磁场分布 r = sqrt(x.^2 + y.^2); theta = atan2(y, x); Bz1 = zeros(size(x)); Br1 = zeros(size(x)); Bx1 = Br1.*cos(theta) - Bz1.*sin(theta); By1 = Br1.*sin(theta) + Bz1.*cos(theta); Bz2 = chi*B0*L/2/mu0*(r.^2-R^2)./((r.^2+L^2/4-2*L/4*r.*cos(theta)).^(3/2)); Br2 = chi*B0*L/4/mu0*L./((r.^2+L^2/4-2*L/4*r.*cos(theta)).^(3/2)).*(L/2-r.*cos(theta)); Bx2 = Br2.*cos(theta) - Bz2.*sin(theta); By2 = Br2.*sin(theta) + Bz2.*cos(theta); Bz = Bz1 + Bz2; Bx = Bx1 + Bx2; By = By1 + By2; % 圆柱体外部的磁场分布 r = sqrt(x.^2 + y.^2); theta = atan2(y, x); phi = atan2(z, r); Bz3 = zeros(size(x)); Br3 = zeros(size(x)); Bx3 = Br3.*cos(theta) - Bz3.*sin(theta); By3 = Br3.*sin(theta) + Bz3.*cos(theta); for n = 0:20 k = n/R; Jnkr = besselj(n, k*r); Jnkp = besselj(n, k*R); Ynkp = bessely(n, k*R); Hnkr = besselh(n, 1, k*r); Hnkp = besselh(n, 1, k*R); Hnkm = besselh(n, 1, k*L/2); An = (Jnkp*Ynkp - Jnkr*Hnkp)./(Jnkp*Hnkm - Hnkp*Jnkm); Bz3 = Bz3 + k*(n+1)*An.*(chi*B0/mu0).*Jnkr.*cos(n*phi); Br3 = Br3 - k*(n+1)*An.*(chi*B0/mu0).*Jnkr.*sin(n*phi); Bx3 = Br3.*cos(theta) - Bz3.*sin(theta); By3 = Br3.*sin(theta) + Bz3.*cos(theta); end % 地球磁场 Bx0 = B0*cosd(inc0)*cosd(dec0); By0 = B0*cosd(inc0)*sind(dec0); Bz0 = B0*sind(inc0); % 总磁场 Bx = Bx + Bx3 + Bx0; By = By + By3 + By0; Bz = Bz + Bz3 + Bz0; end ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值