表征油藏随时间演变的方法和程序的制作方法

专利名称：：表征油藏随时间演变的方法和程序的制作方法
技术领域：
：本发明涉及地球物理学领域，更具体地涉及油类勘探。
背景技术：
：油类勘探的目的是基于从表面或钻井中获得的地球物理学的测量结果，测定载油层的位置。典型地，这些测量涉及向地下发送地震波，并用许多传感器测量各种反射波，其反映了地质结构-表面分离区别材料、断层等等(反射波法地震勘探)。其他的测量在钻井中进行。然后，向地下发送声波、伽马射线或电信号；再用许多传感器感测反射信号。这些传感器可置于地面或海上。典型地，这些技术涉及测量值的处理，以构建地下的图像。通常，该方法需要地震波在地下传播的速度模型。标准成像速度程序包，诸如当前引用的由Paradigm(IL)公司销售的Geodepth软件，使得可能建立这种地震速度场的模型。这些技术可被扩展至允许观测指定的储集层随时间的演变-例如，在油类生产开始之前和油类生产一段时期之后，比较测量的结果。这被称为四维地震(4Dseismic)，其涉及比较不同时间点进行的三维地震测量(3Dseismic)。目标是观测从储集层开采碳氢化合物之后发生的储集层和周围地层的状态的变化。通常，其采集需要比标准三维测量更加仔细，伴随后面的专业处理步骤。J.E.RickeR和D.E.Lumley，时移地震储集层监控的互均衡数据处理墨西哥湾的研究实例(Cross-equalizationdataprocessingfortime-lapseseismicreservoirmonitoringAcasestudyfromtheGulfofMexico)，地球物理学(Geophysics)，第66卷，第4号(2001年7月-8月)中讨论了随时间进行的地震测量中的非重复性的噪音问题。这篇文献公开了将两种实际测量进行配合。移动前数据不能获得。测量的配合包括滤波、振幅平衡和三维变形(3Dwarping)的配合。三维变形在于将视窗内的地震道互关联，以评价适于各测量之间的数据最佳配合的在x、y和t上的移动。在Hall等人，用于时移评估的采用Valhall油田三维托缆和三维海底电缆数据解析变形的相互配合(Cross-matchingwithinterpretedwarpingof3Dstreamerand3Docean-bottom-cabledateatValhallfortime-lapseassessment)，地球物理勘探(GeophysicalProspecting)，2005，53，第283-297页中，公开了传统托缆(streamer)数据和新三维海底电缆数据的相互配合，用于Valhall油田的生产造成的地质力学变化的时移分析。这篇文献主要涉及使用由不同采集方法提供的结果-在Valhall油田的例子中，三维托缆数据和三维海底电缆。该文献表明两种测量使用了类似的移动方案。该方法涉及下列步骤-体积成形，以考虑不同的采集方法；-在体积之间和内部进行振幅平衡；-光谱成形；-整体相互配合，使用局部导算子。两种测量之间的空间漂移通过以迭代法使用三维变形而解析。O.Kolbjornsen和A.R.Syversveen，时间配合-一种评估四维时间漂移的方法(Time-match-amethodofestimating4Dtimeshift)，NorksRegnesentral，noteno.SAND/03/05，2005年4月，其中讨论了一种在四维地震测量中评估时间漂移的方法。所使用的运算法则使新测量中的时间与最初测量的时间相配；按逐条地震道(tracebytrace)提供了从一种至另一种的图。更具体地，通过局部压缩和伸展如监控测量地震道的时间轴以最小化振幅之间的方差而进行配合。这篇文件讨论了合成试验案例。这些在先技术的文献教导了三维变形，比较地震测量的重新排列以补偿采集中的缺陷(或地震测量的非重复性)和地下速度的变化。基于关联的方法的一个问题是关联窗口的尺寸。如果用于关联的窗口太大，关联的精度性就要受到影响事实上，相关值不仅取决于所考察点的测量之间的不同，除所考察的点之外还取决于其他作用。如果用于关联的窗口太小，关联易受噪音和测量的非重复性的严重影响，包括由于需要观测的作用的变化。目前仍然需要能减轻这个问题的一种表征储集层的随时间的演变的方法。
发明内容因此，在一个实施例中，通过共分析沿地下传播途径的地震子波的传播时间和地震振幅的变化，本发明提供了一种表征处于生产过程中的油藏变化的方法，所述方法包括下列步骤-以与第一速度Vb关联的第一时间T时的一组地震道，提供储集层的基础测量；-在与第二速度Vm关联的第二时间T+ΔT，以与所述基础测量中相同位置相关的一组地震道，提供所述储集层的监控测量；-对于所述基础测量中的一组点，计算该一组点的一种差值的范数(norm)的总和S，所述差值的范数为在基础测量中的所述点i上的地震道的振幅bi与在监控测量中的相应时间点i’上的地震道的振幅mi’和由第一速度Vb和第二速度Vm间的差值导致的所述相应时间点i’上的局部反射率变化而造成的振幅的总和之间的差值；其中，所述相应时间点i’的随时间的漂移为沿着从表面至所述相应间点i’的传播路径的速度变化得到的时间漂移；-通过最小化所述总和S，表征该油藏的演变。在一个实施例中，由于所述相应时间点的局部反射率变化造成的振幅在测量中使用的子波长度的1至5倍的时间范围内计算。还可规定，相应时间点i’仅在时间上漂移。在这种情况中，优化步骤可按照逐个地震道的基础进行。在另一个实施例中，相应时间点i’沿传播路径在时间和空间上漂移。在计算步骤之前，所述方法可进一步包括测量的零偏移化步骤。在另一个实施例中，本发明提供了一种位于计算机可读介质上的计算机程序，包括适于在计算机上运行该方法的所有步骤的计算机程序代码装置。下面，参考附图，并通过非限制性实施例，对本发明的方法进行描述，其中图1是地震块(seimicblock)的示意图，为了清楚仅示出一条地震道；图2是本发明的一实施例中的方法的流程图；图3示出了在图2方法的二维合成测试(2Dsynthetictest)中，用于基础测量的地震块的截面；图4示出了所述合成测试中施加的速度变化；图5示出了图4的速度变化造成的振幅的变化；图6示出了图2的方法的结果；图7示出了在优化步骤期间计算的速度变化；图8示出了另一个二维合成测试的在解析的沙层中的速度变化位置；图9示出了图2的方法获得的结果；图10是在实际例子中，监控测量和基础测量之间原始差异的截面；图11示出了根据图2的方法计算的地震振幅差异；图12示出了在相同区域主要承包商(leadingcontractor)估计的阻抗差异；图13示出了在图10和11的截面的部分的成比例的速度变化在本说明书的剩余部分中，使用术语“基础测量”和“监控测量”表示储集层的地震测量。假定基础测量比监控测量在较早的时间进行。具体实施例方式本发明基于的事实是由于勘探，储集层的变化将造成速度场的变化。事实上，油类将被气体或水替代且流体压力发生变化，造成密度和模量的变化，从而速度发生变化。这些速度的变化将导致处于这些变化之下的反射体的地震响应发生时间漂移和反射率的相关变化，造成局部波形的变化。在一个具体的实施例中，本发明提出在监控测量中评估这些作用。这使之可能从对基础测量和监控测量的比较中得出速度变化场，而不用进行地震道的互相关性处理。当预期密度的变化相对较小且有效反射角度很小(和/或预期的剪切波速度的变化也相对较小)时，这个方法尤其有效。为了有利于计算，进一步有利地假设时间漂移仅来自速度的变化以及采集或处理参数上的变化可忽略。后一个假设在现代专用的四维测量(4Dsurvrys)中日益被满足。当作用于框架模量上的压力为主要时移现象时，或者类似地，当少量的气体被释放/导入至先前100％液体孔隙流中时，前一个假设可实现。如下面的实施例中所示，所述方法尤其可应用于分析邻近的偏移叠加的时间变形。图1是地震块的示意图，为了清楚，图中仅显示了一条地震道。术语地震块用于描述加工成图像后指定地理区域上的一组测量。如前所知，人们使用一组相互垂直的标准坐标，其中x轴和y轴位于水平面上。相应于时间的z轴是垂直的且向下延伸。通常对于震波测量而言，人们使用坐标(x，y，t)作为测量的时间表示，或者坐标(x，y，z)作为深度移动后的测量的空间表示。将一组传感器Ci置于地表或海上，在空间坐标(xi，yi，zi)点中，i表示传感器号的整数；尽管大部分文献似乎认同zi＝0的假设，所述传感器几乎并不刚好置于zi＝0的位置。典型地，托缆(streamer)在5-7m的深度被拖行，海底电缆的位置如其名称所示，甚至有时陆地地震检波器埋在几米深的地下。当进行测量时，各传感器Ci记录原始信号；这个原始信号表示地表下各种界面反射的地震波。随后，传感器接收的原始信号被处理以提供包括垂直地震道的集合的地表下图像，该垂直轴表示时间t或深度z。图1示出了该组坐标的轴x、y和t(或z)，以及伴有相应地震道的一传感器Ci，所述地震道在图中标记为2。为了清楚，图1仅示出了一个传感器和一条地震道，然而通常测量涉及多个传感器和超百万的多条地震道。如前所知，地震处理将地震事件尽可能精确地置于在其真实水平位置，有效地与传感器的原始位置无关。在零倾斜度和零偏移的理想情况中，由于无横向速度梯度，所述地震道将真正地位于传感器之下，但这是一个相当特殊的例子；通常，地震道还概念上地与实际传感器之间的位置相关。这些技术的详情可在如下文献中获得zdoganYilmaz，地震数据处理(SeismicDataProcessing)，地球物理学勘探协会刊物(SocietyofexplorationGeophysicists)，1987。图2是根据本发明的一个实施例的方法的流程图。在步骤12中，提供了在第一时间T时的一组地震道作为储集层的基础测量。对于一个指定的地震道，该基础测量提供了振幅b(t)，其为时间t时的振幅函数；如果该地震道被采样，简单使用一组值ti，其中i为有许多值的指数，典型地具有对于4ms采样周期的地震道采样的的约1000个值。随后，该地震道被处理为一组b(ti)或bi值。在步骤16中，以一组地震道提供了在第二时间T+ΔT进行的储集层的监控测量。在最简单的假设中，ΔT为正值，该监控测量在基础测量之后的时间进行；然而，对于本发明方法的操作而言，测量进行的顺序是不相关的，并且理论上，该时移ΔT也可为负值-其相当于将较早的测量与较后的测量进行比较。对于基础测量而言，在监控测量中采样的地震道被表示为一组m(ti)或mi值。理想地，监控测量中的地震道与基础测量中的相同位置联系。如可能，这可通过使用相同设备和相同的方法运行基础测量和监控测量而进行。具体而言，5-10m的位置之间差异仍然获得可接受的结果。一些技术，如插值法可以使用，如果监控测量和基础测量的地震道不满足这个条件(Eiken，O.等人，2003，一种获取高重复性的托缆地震数据的已验证的方法(Aprovenmethodofacquiringhighlyrepeatabletowedstreamerseismicdata)，地球物理学(Geophysics)，68，1303-1309。在这个实施例中，本发明导致估算相对慢度变化，n，其中慢度为速度的倒数，其中n＝-ΔV/V＝(Vb-Vm)/Vb在该公式中，Vb和Vm为考虑用于三维形变的概念上的局部垂直速度。通常它们不与现有地震处理方法中使用的任何速度，即移动速度或叠加速度匹配。它们的差表示基础测量和监控测量之间的地震事件的垂直时间漂移。上述相对慢度变化，n，可在地震块的各采样中，即在地震道的各采样中进行评估。为了估算所述相对慢度变化，对一组位点进行最优化，如下面的解释。在指定的地震道上，对于指定的采样i，时间漂移wi(以采样的单位)可表达如下wi=Σk=0ink]]>其中nk为采样k的相对慢度变化。这个表达式表示的事实是地震道上采样i的时间漂移wi是由该采样上的速度变化造成。严格的说，所述时间漂移是信号从源头至考察的采样以及返回的路径上慢度的积分变化。上面给出的表达式基于的假设是时间漂移来源于所考察的采样上的速度变化；这相当于从源头至反射体以及返回至传感器的垂直或准垂直的传播途径。当倾斜度为零或受限时，在零偏移下，即发送体和传感器之间的距离为零或者相比于反射体的垂直深度可以忽略不计，这个条件可被完全满足；所述倾斜度为水平面和局部反射体之间形成的角度。具体而言，因为我们通常仅关心通过储集层中的出现生产相关性变化的区域的传播，所以这些假设仅需要应用于所述储集层的厚度上；此外，“零”横向位移相当于保持在一地震箱(seismicbin)内。即，对于厚度约100m的储集层以及25m的地震箱尺寸，这些假设允许传播路径具有达12.5m的在储集层带的入口(沿自源头向下的路径)和反射点之间的横向位移，以及返回到传感器的路径上类似的位移，或者对于零偏移地震道的达14°的倾斜度。在预期地震箱之间有一些连续性的大部分情况中，以及在典型的具有100m水平尺寸的储集层模型网格单元的环境中，这些假设可进一步放松。在这些假设下，基础测量至监控测量在速度上的变化将影响指定地震道内的振幅，以使得b(ti)=m(ti+wi)+ψ*n·(ti)]]>其中ψ为地震子波，n·(ti)=n(ti)-n(ti-1).]]>第一项表示在上面给定的假设下的样点i上的速度变化导致的时间漂移。第二项表示在速度变化后的局部反射率变化对地震道的影响；在这第二项中，局部变化被认为是在与子波相同的时间范围内，即在与子波持续时间相等的时间范围内。本发明提出了评价从基础测量到监控测量的速度变化，因此通过评价多个基础测量点的以下差值的范数(norm)的总和而表征储集层的演变Δi=b(ti)-m(ti+wi)-ψ*n·(ti)]]>这个差值Δi为基础测量中位点i处的地震道的振幅bi与-在监控测量中相应时间点i’处的地震道的振幅mi’和-由于第一速度Vb和第二速度Vm之间的差值导致的所述相应时间点i’上局部反射率变化而造成的振幅扰动的总和之间的差异。相应时间点i’随时间漂移为以来源于沿表面至所述点i’及返回的传播途径的速度变化得到的时间漂移wi。不用明言，ti和t’i也表示在采样中。各点上的差值的范数的总和S，S=Σi=0N||Δi||]]>通过改变速度变化-表示为相对慢度变化ni，而使之最小化。这为各点提供了速度变化的场。该速度变化场确定了使监控测量对应基础测量的变形操作的参数，且还可用于直接表征储集层的演变。在图2的方法的步骤16中，选择一组点；所述总和S在这组点上被最小化。根据计算资源，可改变该一组点的大小，但通常会选择完全包括所考察的储集层的全部体积。在下面提供的例子中，使用来自全部基础和监控测量的地震道，所述地震道以时间作视窗横跨储集层。这提供了全部测量上的速度变化值。在所述方法的步骤18中，计算出初始的总和S值。在所述方法的步骤20中，通过改变相对速度变化的值，所述总和S值被最小化。下面提供了最优化技术的一个例子；然而，还可使用在本领域中已知的其他最优化技术，如模拟退火(simulatedannealing)。如前所述，如果在监控测量中的地震道与在基础测量中相同位置关联，位点仅发生时间漂移。那么可在逐条地震道的基础上进行计算；换句话说，最优化操作可单独对各个地震道进行。这简化了计算，且更容易作为平行任务在许多计算机上运行最优化步骤。在步骤22中，总和S被最小化，并提供了进行最优化的该组位点中不同位点的速度变化的值。最小化的总和S表征了储集层随时间的演变。步骤20中的求最小值可使用Gauss-Newton公式进行。Gauss-Newton公式为现有公知。申请人进行的测试表明收敛(convergence)通常在2-4次反复之后实现。已经确定的是如果时间漂移总计小于主频的半周期，所述方法会收敛，正如通常的情况的那样。超过这个值，会有向局部最小值收敛的风险。所述方法的上述实施例可能收敛于局部最小值的情况并不使该方法无效，因为使用，例如标准关联方法，正确选择相对慢度变化的初始值，以使得剩余的漂移小于半周期，可允许其随后的应用。可选地，使用整体优化方法可确保向真实的最小值收敛。在认为密度变化不可忽略的情况中，除统计学上与速度变化相关外，子波可相应地依比例调节。例如，如果存在正相关，以使得速度平均变化1％意味着密度上0.25％的变化，子波可以按1.25的倍数比例调节以给出最能代表由速度扰动造成的地震道的变化。图2的方法被应用于测量Girassol油田。基础测量由CGG以2ms采样的高分辨率测量在1999年进行。生产开始于2001年，在2002年使用与基础测量相同的设备进行监控测量。各测量包括2,500,000条地震道，各大约2000个采样，在如上详述的实施例中，其以500个采样作视窗，用于此处所述的方法的应用。图2的方法使用最初为基础测量造模的速度场进行；为了比较时移，处理监控测量的速度场与基础测量的相同。因为基础测量和监控测量在非常相似的条件下进行，最优化操作可在逐条地震道的基础上进行。所述方法在平行的60台处理器上进行20小时，即约1200CPU小时；在最优化过程中的3次迭代中达到收敛。在图2的实施例中，假定在监控测量中在地震道中表示的事件与在基础测量中的处于相同横向位置。从而可以认为，由于沿地震道(假定零偏移，零倾斜度)的速度的变化和反射率的相关变化，地震道上的变化仅由时间漂移造成。然而，本发明的方法还可应用于倾斜度不可忽略的情况，以使得地震道上观察到的变化不仅由时间漂移造成，还由横向漂移造成。换句话说，如果在监控测量中在地震道中表示的事件的横向位置相对于在基础测量中的那些发生扰动，本发明的方法仍然可以进行。在这种情况中，差值Δi的评价采用来自监控地震块的振幅值，其不仅被时间漂移位移，还被横向位置的变化位移。移动的图像之间的总位移量归因于沿不再垂直的传播路径的速度变化；对于一个指定的点i，那么相应的点i’沿由移动过程限定的传播路径，以相应于沿该路径的慢度变化的积分的量，在时间和空间上漂移，从而产生图像。和前面一样，由于与该点邻域中的速度变化相关的局部反射率的变化，在i’上有振幅的变化。相比于先前所述的实施例，这需要一些附加的工作以限定这些传播路径和计算相关的积分。类似的方法可应用于非零偏移的情况；甚至在倾斜度很小的情况，从而所述位移主要为垂直的，所述时间漂移与沿传播路径的慢度积分相关，所述传播路径大部分位于包含被比较的振幅采样的地震箱之外。图3-7示出了图2方法的二维合成测试；通过将速度变化应用于实际数据，基于速度变化计算该监控测量，然后将该方法应用于重构该速度变化，而进行该测试。图3示出了用于基础测量的地震块的截面。这截面是来自Girassol油田的实际数据，通过倒置基础地震数据得到的纵波(p-wave)阻抗立方的褶积造模而获得。如参考图1中所述，水平轴示出地震道数-或者沿该截面的距离；垂直轴代表时间。在图3中显示了各地震道的振幅；可以看出储集层的倾斜度几乎为零，其中倾斜度的最大值约3°且储集层中倾斜度的平均值约1.5°。图4示出了施加于该合成测试中的速度变化。该速度变化呈“蝴蝶”形，其使该倒置结果容易理解，及使该方法具有可能的分辨率和稳定性。在“蝴蝶”的左侧，速度变化为正值，具有+8％的不变值；在“蝴蝶”的右侧，速度变化为负值，具有-8％的不变值。图5示出了图4的速度变化造成的振幅的变化。具体地，图5由下列产生-由于图4的速度变化，对图3的数据引起的振幅的变化的计算；以等于基础速度场与图4的“蝴蝶”变化的总和的速度场，用传播工具计算该变化；-计算的振幅和图3的振幅之间的差值的计算。图5示出了时间漂移和反射率对振幅的影响，为基础测量和监控测量之间的原始差异。在图4的“蝴蝶”形状以上的振幅没有变化，因为速度无变化。由于反射率和时间漂移的变化，“蝴蝶”形状中的速度变化将造成振幅变化。在“蝴蝶”形状的下方，振幅上的变化由时间漂移造成-因为局部反射率相同。图6示出了图2的方法的结果。其示出了在最小化总和S之后的在基础和漂移监控测量之间振幅的差异，揭示了估计的振幅变化仅归因于反射率的变化。在该实施例中，按逐条地震道进行最优化操作，并在三次迭代中收敛。图6示出振幅的变化基本反映了图3速度的变化；这表明了该方法的有效性。图7示出了在最优化步骤期间计算的速度的变化。尽管仍有一些错误，图7示出的速度变化的值基本为图3的那些。速度扰动的光谱分析显示该倒置具有宽带性质，频率从0Hz回复至超出地震频谱的上限，因而倒置的速度变化的特性可以为定量的和易于解析的。图8和9示出了图2方法的更加真实的2D合成测试；该测试通过将速度变化应用于实际数据，以更复杂的速度变化而进行。具体而言，速度变化以解析的沙层排列，如图8所示。图9示出了本发明的方法中获得的结果；速度的变化仍基本相当于图8的那些。在图3-7的实施例中，估算的变化的频谱验证了计算的宽带性质。图8和9表明该方法还可以比图3更复杂的速度变化操作。图10-13示出了由图2的方法对Girassor油田的实际测量所获得的结果。图10为Girassol接近垂直叠加(substrack)的快速路径(fast-track)上的监控测量和基础测量之间原始差异的截面；图10是沿所述测量的交叉线4001取得；图11示出了在图2的方法中计算的地震振幅差异；相比于图10，可以看出现在的差异为图的上部的主要可见部分；这部分相应于取回石油的储集层的部分。在该储集层的较深部分，振幅的差异在图11的截面中比在图10的截面中更小；这是由于图2所述的方法考虑了在上部水平面产生的时间漂移而造成的。图13示出了在图10和11的钻井附近的部分截面上的成比例的速度变化，如在图2的最优化步骤中计算的那样；这可与图12进行比较，图12显示了在相同区域主要承包商估算的阻抗差异。源于速度变化特性的宽带性质的可解析性的改进是明显的，我们对于该特性的定量值的信心也得到提高。由于无噪音，图13还证实了由这种方法计算速度变化特性的稳定性，相比于基于其他前述方法的计算，我们认为有明显的改善。图10-13的实施例表明本发明的方法提供了立即可用且与勘探数据一致的结果。图10-12的光谱分析还表明，如先前的实施例中那样，很好地回收了较低的频率。图10-12的实施例在接近垂直叠加的Girassol油田上进行，该位置平均入射角等于12度；因而，即使变形位移可为垂直的，其并不正好等于叠加的慢度变化的垂直积分；图10-12的实施例还说明零偏移和零倾斜度的假设仍可提供代表实际勘探数据的结果。本发明的方法可以计算机程序实施。该程序适于接收基础和监控测量的数据，以及速度场的数据；这些数据采用由如前所述的现有技术计算机数据包所提供的格式。所述程序运行图2的方法的各种步骤。权利要求1.一种表征处于生产过程中的油藏演变的方法，所述方法通过共分析沿地下传播途径的地震子波的传播时间和地震振幅上的变化来实现，其包括下列步骤-以与第一速度Vb关联的第一时间T时的一组地震道，提供储集层的基础测量；-在与第二速度Vm关联的第二时间T+ΔT，以与所述基础测量中相同位置相关的一组地震道，提供所述储集层的监控测量；-对于所述基础测量中的一组点，计算该组点的差值的范数的总和S，所述差值的范数为在基础测量中的所述点i上的地震道的振幅bi与在监控测量中的相应时间点i’上的地震道的振幅mi’和由第一速度Vb和第二速度Vm间的差值导致的所述相应时间点i’上的局部反射率变化而造成的振幅的总和之间的差异；其中，所述相应时间点i’的随时间漂移为沿从表面至所述相应时间点i’的传播路径的速度变化得到的时间漂移；-通过最小化所述总和S，表征该油藏的演变。2.根据权利要求1所述的方法，其特征在于，由所述相应时间点上的局部反射率变化而造成的的振幅是在测量中使用的子波长度的1至5倍的时间范围内计算。3.根据权利要求1或2所述的方法，其特征在于，相应时间点i’仅在时间上漂移。4.根据权利要求3所述的方法，其特征在于，在逐条的地震道基础上进行优化步骤。5.根据权利要求1或2所述的方法，其特征在于，相应时间点i’沿传播路径在时间和空间上漂移。6.根据权利要求1至5中任一权利要求所述的方法，其特征在于，在计算步骤之前，还包括测量的零偏移化步骤。7.一种位于计算机可读介质上的计算机程序，其包括适于在计算机上运行权利要求1至6中任一权利要求所述的方法的所有步骤的计算机程序代码装置。全文摘要通过共分析沿地下传播途径的地震子波的传播时间和地震振幅上的变化，进行处于生产过程中的油藏的演变分析。以与第一速度V文档编号G01V1/40GK101086535SQ20071011046公开日2007年12月12日申请日期2007年6月5日优先权日2006年6月6日发明者保罗·威廉森,保罗·塞克斯顿,亚当·约翰·谢雷申请人:托塔尔公司