1.一种用于在介质中从非线性相互作用产生微地震现象以便描述所述介质的特征的方法，所述方法包括：
由第一声源产生包含排列为时间序列的第一多个脉冲的第一编码声信号，所述第一多个脉冲在时间上分开，每个脉冲都包含在中心频率的被调制信号，其中两个相继脉冲的中心频率不同；
由第二声源产生包含排列为时间序列的第二多个脉冲的第二编码声信号，所述第二多个脉冲在时间上分开，其中两个相继脉冲的中心之间的时间分离与所述第一多个脉冲中两个对应脉冲的中心之间的时间分离相同，其中在传播所述第二多个脉冲的起始时间与传播所述第一多个脉冲的起始时间之间提供了起始时间差，其中每个脉冲都包含被调制信号并且所述第二多个脉冲中的每个脉冲内所述被调制信号的中心频率是所述第一多个脉冲中的对应脉冲的被调制信号的中心频率的选定分数d；
其中所述第一声源和所述第二声源是可控的，使得所述第一声信号和所述第二声信号的轨迹在所述介质内的混合区中相交;
由接收器接收包括在所述混合区中由非线性过程从所述第一声信号和所述第二声信号产生的第三信号的检测信号；
由处理器对所接收的信号进行数据处理，或与编码信号模板相关，或者兼而有之，以提取噪音之上或由线性相互作用过程产生的信号之上，或者兼而有之，由所述非线性混合过程产生的第三信号，以获得在所述第一声信号和第二声信号的混合区发生的仿真的微地震现象信号；以及
根据所述仿真的微地震现象信号，描述所述介质性质的特征或产生所述介质性质的3D图像，或者兼而有之。
2.根据权利要求1的方法，其中，产生所述第一编码声信号包括产生所述第一声信号，使得所述第一多个脉冲中两个相继脉冲的中心之间的时间分离大于每个脉冲的持续时间。
3.根据权利要求1的方法，其中，由所述非线性混合过程产生并在接收器接收的所述第三信号包括以时间序列到达并在时间上分开的第三多个脉冲，其中两个相继脉冲的中心之间的时间分离与所述第一多个脉冲中两个相继脉冲的中心之间的时间分离相同；
其中所述第三多个脉冲中的每个脉冲都包含被调制信号，该被调制信号的每个脉冲的中心频率都等于所述第一多个脉冲的对应脉冲的中心频率与所述第二多个脉冲的对应脉冲的中心频率之间的差；
其中所述第三多个脉冲的每个脉冲到达接收器的时间是相对于产生所述第一多个脉冲的对应脉冲延迟从所述第一声源到所述混合区中心的传播时间和从所述混合区中心到接收器的传播时间的总和的时间。
4.根据权利要求1的方法，进一步包括对起始时间差的范围和频率分数d的范围重复产生所述第一信号、产生所述第二信号、接收所述第三信号，和执行数据处理或与编码信号模板相关，或者兼而有之，以便从多个混合区产生的仿真微地震现象获得信号。
5.根据权利要求1的方法，其中，产生所述第一信号包括产生作为多个脉冲之和的第一信号，其中每个脉冲具有的信号振幅都等于包络函数与被调制信号函数的乘积；以及其中产生所述第二信号包括产生作为多个脉冲之和的第二信号，其中每个脉冲具有的信号振幅都等于包络函数与被调制信号函数的乘积。
6.根据权利要求1的方法，其中，产生所述第一信号和所述第二信号包括产生第一信号u1(t)和第二信号u2(t)，形式如下：
以及
其中m是与脉冲相关联的索引号；
∑表示对索引m=1至M求和，M是等于或大于1的整数；
E1m(t-tm)是所述第一声信号u1的脉冲m的振幅包络；
E2m(t-tm-δ)是所述第二声信号u2的脉冲m的振幅包络；
ωm是所述第一声信号u1的脉冲m的调制信号的中心频率；
(d＊ωm)是所述第二声信号u2的脉冲m的调制信号的中心频率；
d是所述第一多个脉冲中所述频率ωm与所述第二多个脉冲中所述频率d＊ωm之间的频率比，其中d是正实数；
δ是产生第一声信号u1与所述第二声信号u2之间的所述起始时间差；
exp(iωm*(t-tm)是所述第一声信号u1的所述脉冲m内的调制信号；
exp(id*ωm(t-tm-δ)是第二声信号u2的所述脉冲m内的调制信号；
tm是在所述第一声信号u1中产生脉冲m的时间；
tm+δ是在所述第二声信号u2中产生脉冲m的时间；以及
exp(iζm)是在所述第一信号u1或第二信号u2内每个脉冲m的相位项。
7.根据权利要求1的方法，其中，由非线性混合过程产生的所述信号u3(t)形式如下：
其中：
E3m（t-tm-T1-T3）是通过非线性混合所述第一信号和所述第二信号而产生的所述第三多个脉冲信号u3的脉冲的包络；
exp(i*(1-d)*ωm*(t-tm–T1-Τ3))对应于所述第三多个脉冲信号中第m个脉冲内的第三调制信号；
(1-d)*ωm是所述第三多个脉冲中第m个脉冲内的所述调制信号的频率，其是所述第一多个脉冲的第m个脉冲内的所述第一调制信号的中心频率ωm与所述第二个多个脉冲的第m个脉冲内的所述第二调制信号的中心频率d＊ωm之间的差；
d是选定的频率比；
t是所述信号时间；
R是接收器的标志或接收器的位置；
（T1+T3）是所述第三多个脉冲中第m个脉冲到达接收器的时间；以及
ζm是每个脉冲m的相位而exp(iζm)是每个脉冲m的相位项。
8.根据权利要求1的方法，进一步包括由所述处理器提取由所述非线性混合过程产生的信号，以便通过将所述检测信号与模板信号相关，从所述第一声信号和第二声信号的混合区的所述仿真的微地震现象获得所述信号，其中所述模板信号包括多个脉冲，其中所述模板信号中所述多个脉冲的中心之间的时间分离与所述第一声信号的多个脉冲中的中心之间的时间分离相同，并且其中在所述第三信号中所述多个脉冲的中心频率是第一声信号中所述多个脉冲的中心频率和所述频率分数d的函数。
9.根据权利要求8的方法，其中，进一步包括增大所述第一多个脉冲的数量和增大所述第二多个脉冲的数量，对噪音之上或由线性相互作用过程产生的信号，或者兼而有之，由所述非线性混合过程产生的第三信号增强分辨力。
10.根据权利要求1的方法，进一步包括由所述处理器提取由所述非线性过程产生的第三信号，以便通过将所述检测信号与模板信号相关，从所述第一声信号和第二声信号的混合区的所述仿真的微地震现象获得所述信号，其中所述模板us形式如下：
其中：
Wm(t-tm)是所述模板信号中第m个脉冲的振幅包络；
exp(i*g(ωm)*(t-tm))对应于所述模板信号的多个脉冲中第m个脉冲内的调制信号；
g(ωm)是在所述模板的多个脉冲中第m个脉冲内，取决于频率分数d的调制信号的选定函数；
t是传播时间；
tm是所述多个脉冲中第m个脉冲被模拟到达接收器的时间；以及
ζm是每个脉冲m的相位而exp(iζm)是每个脉冲m的相位项。
11.根据权利要求10的方法，进一步包括：
通过使在接收器Ri的所述检测信号与编码模板信号us(t)相关，由所述处理器从所述第一声信号和第二声信号的混合区的所述仿真的微地震现象提取所述第三信号；以及对于所述第二多个脉冲与所述第一多个脉冲之间的多个起始时间差δ以及所述第二多个脉冲中的每个脉冲内的所述调制信号中心频率与所述第一多个脉冲中的每个对应脉冲内的所述调制信号中心频率之间的多个频率比d，重复所述模板信号与所述检测信号之间的相关，对于每个起始时间差δ和每个频率比d产生对每个接收器Ri的相关信号，表示为M(Ri,t)，它包含在根据非线性相互作用的选择规则的特定混合位置由非线性相互作用产生的仿真的微地震现象。
12.根据权利要求11的方法，进一步包括使用基于压缩或剪切波速度的速度模型的轨迹分析和传播时间计算，对于每个起始时间差δ和对于每个频率比d以及所述对应的相关信号M(Ri,t)，计算仿真的微地震现象的空间坐标，其中所述两个声信号非线性地相互作用。
13.根据权利要求11的方法，进一步包括应用逆时运算，使所述仿真的微地震现象的信号以传播速度模型在时间上传播回到所述微地震现象起源的时间和位置。
14.根据权利要求13的方法，进一步包括确定所述微地震现象在所述起源位置的相对强度，其中在所述起源位置的相对强度与所述微地震现象的所述起源位置处所述介质的非线性性质成正比。
15.根据权利要求11的方法，进一步包括通过使用基尔霍夫成像、射束成像、波动方程成像或者逆时方法，从包含所述第一声信号和第二声信号的混合区的所述仿真的微地震现象的所述相关信号来产生所述微地震现象或所述介质的非线性性质的三维图像，或者兼而有之，其中，所述相关信号从在多个接收器接收的信号获得，其中所述第一源和所述第二源位于不同位置。
16.根据权利要求11的方法，进一步包括通过把所述第一源放置在多个位置，多次重复地由所述第一源产生所述第一信号，以及通过把所述第二源放置在多个位置，多次重复地由所述第二源产生所述第二信号；以及通过把接收器放置在多个位置，多次重复地接收所述第三信号；以及通过使用基尔霍夫成像、射束成像、波动方程成像或者逆时方法，从包含所述第一声信号和第二声信号的混合区的所述仿真的微地震现象的所述相关信号来产生所述介质的非线性性质的多个三维图像，其中，所述相关信号从在多个接收器接收的信号获得。
17.根据权利要求16的方法，进一步包括通过速度层析成像方法更新所述传播速度模型，使按照从所述第一声源和所述第二声源的多个位置接收的信号的所述多个三维图像的差异最小。
18.根据权利要求17的方法，进一步包括迭代地更新所述传播速度模型和成像直到获得的多个图像之间的差异最小，以及结合所述多个图像以产生所述介质的非线性性质的最终三维图像或所述局部速度比Vp/Vs或者兼而有之。
19.一种用于在介质中从非线性相互作用产生微地震现象以便描述所述介质的特征的系统，所述系统包括：
第一声源，被配置为产生包含排列为时间序列的第一多个脉冲的第一编码声信号，所述第一多个脉冲在时间上分开，每个脉冲都包含在中心频率的被调制信号，其中两个相继脉冲的中心频率不同；
第二声源，被配置为产生包含排列为时间序列的第二多个脉冲的第二编码声信号，所述第二多个脉冲在时间上分开，其中两个相继脉冲的中心之间的时间分离与所述第一多个脉冲中两个对应脉冲的中心之间的时间分离相同，其中在传播所述第二多个脉冲的起始时间与传播所述第一多个脉冲的起始时间之间提供了起始时间差，其中每个脉冲都包含被调制信号并且所述第二多个脉冲中的每个脉冲内所述被调制信号的中心频率是所述第一多个脉冲中的对应脉冲的被调制信号的中心频率的选定分数d；
其中所述第一声源和所述第二声源是可控的，使得所述第一声信号和所述第二声信号的轨迹在所述介质内的混合区中相交;
接收器，被配置为接收包括在所述混合区中由非线性过程从所述第一声信号和所述第二声信号产生的第三信号的检测信号；
处理器，被配置为对所接收的信号进行数据处理，或与编码信号模板相关，或者兼而有之，以提取噪音之上或由线性相互作用过程产生的信号之上，或者兼而有之，由所述非线性混合过程产生的第三信号，以获得在所述第一声信号和第二声信号的混合区中发生的仿真的微地震现象信号，以便基于所述仿真的微地震现象信号，描述所述介质性质的特征或产生所述介质性质的3D图像，或者兼而有之。
20.根据权利要求19的系统，其中所述第一声源被配置为产生所述第一编码声信号，使得所述第一多个脉冲中两个相继脉冲的中心之间的时间分离大于每个脉冲的持续时间。
21.根据权利要求19的系统，其中，由所述非线性混合过程产生并在接收器接收的所述第三信号包括以时间序列到达并在时间上分开的第三多个脉冲，其中两个相继脉冲的中心之间的时间分离与所述第一多个脉冲中两个相继脉冲的中心之间的时间分离相同；
其中所述第三多个脉冲中的每个脉冲都包含被调制信号，该被调制信号的每个脉冲的中心频率都等于所述第一多个脉冲的对应脉冲的中心频率与所述第二多个脉冲的对应脉冲的中心频率之间的差；
其中所述第三多个脉冲的每个脉冲到达接收器的时间是相对于产生所述第一多个脉冲的对应脉冲延迟从所述第一声源到所述混合区中心的传播时间和从所述混合区中心到接收器的传播时间的总和的时间。
22.根据权利要求19的系统，其中，所述处理器被配置为对起始时间差的范围和频率分数d的范围，控制所述第一源和所述第二源重复产生所述第一信号、产生所述第二信号，控制接收器重复接收所述第三信号以及重复执行数据处理或与编码信号模板相关，或者兼而有之，以便从多个混合区产生的仿真微地震现象获得信号。
23.根据权利要求19的系统，其中所述第一声源被配置为产生作为多个脉冲之和的第一信号，其中每个脉冲具有的信号振幅都等于包络函数与被调制信号函数的乘积；以及其中所述第二声源被配置为产生作为多个脉冲之和的第二信号，其中每个脉冲具有的信号振幅都等于包络函数与被调制信号函数的乘积。
24.根据权利要求19的系统，其中所述第一信号包括第一信号u1(t)以及所述第二信号包括第二信号u2(t)，形式如下：
以及
其中m是与脉冲相关联的索引号；
∑表示对索引m=1至M求和，M是等于或大于1的整数；
E1m(t-tm)是所述第一声信号u1的脉冲m的振幅包络；
E2m(t-tm-δ)是所述第二声信号u2的脉冲m的振幅包络；
ωm是所述第一声信号u1的脉冲m的调制信号的中心频率；
(d＊ωm)是所述第二声信号u2的脉冲m的调制信号的中心频率；
d是所述第一多个脉冲中所述频率ωm与所述第二多个脉冲中所述频率d＊ωm之间的频率比，其中d是正实数；
δ是产生第一声信号u1与所述第二声信号u2之间的所述起始时间差；
exp(iωm*(t-tm)是所述第一声信号u1的所述脉冲m内的调制信号；
exp(id*ωm(t-tm-δ)是第二声信号u2的所述脉冲m内的调制信号；
tm是在所述第一声信号u1中产生脉冲m的时间；
tm+δ是在所述第二声信号u2中产生脉冲m的时间；以及
exp(iζm)是在所述第一信号u1或第二信号u2内每个脉冲m的相位项。
25.根据权利要求19的系统，其中，由非线性混合过程产生的所述信号u3(t)形式如下：
其中：
E3m（t-tm-T1-T3）是通过非线性混合所述第一信号和所述第二信号而产生的所述第三多个脉冲信号u3的脉冲的包络；
exp(i*(1-d)*ωm*(t-tm–T1-Τ3))对应于所述第三多个脉冲信号中第m个脉冲内的第三调制信号；
(1-d)*ωm是所述第三多个脉冲中第m个脉冲内的所述调制信号的频率，其是所述第一多个脉冲的第m个脉冲内的所述第一调制信号的中心频率ωm与所述第二个多个脉冲的第m个脉冲内的所述第二调制信号的中心频率d＊ωm之间的差；
d是选定的频率比；
t是所述信号时间；
R是接收器的标志或接收器的位置；
（T1+T3）是所述第三多个脉冲中第m个脉冲到达接收器的时间；以及
ζm是每个脉冲m的相位而exp(iζm)是每个脉冲m的相位项。
26.根据权利要求19的系统，其中，所述处理器被配置为提取由所述非线性混合过程产生的信号，以便通过将所述检测信号与模板信号相关，从所述第一声信号和第二声信号的混合区的所述仿真的微地震现象获得所述信号，其中所述模板信号包括多个脉冲，其中所述模板信号中所述多个脉冲的中心之间的时间分离与所述第一声信号的多个脉冲中的中心之间的时间分离相同，并且其中在所述第三信号中所述多个脉冲的中心频率是第一声信号中所述多个脉冲的中心频率和所述频率分数d的函数。
27.根据权利要求19的系统，其中，所述处理器被配置为提取由所述非线性过程产生的第三信号，以便通过将所述检测信号与模板信号相关，从所述第一声信号和第二声信号的混合区的所述仿真的微地震现象获得所述信号，其中所述模板us形式如下：
其中：
Wm(t-tm)是所述模板信号中第m个脉冲的振幅包络；
exp(i*g(ωm)*(t-tm))对应于所述模板信号的多个脉冲中第m个脉冲内的调制信号；
g(ωm)是在所述模板的多个脉冲中第m个脉冲内，取决于频率分数d的调制信号的选定函数；
t是传播时间；
tm是所述多个脉冲中第m个脉冲被模拟到达接收器的时间；以及
ζm是每个脉冲m的相位而exp(iζm)是每个脉冲m的相位项。
28.根据权利要求29的系统，进一步包括：
其中，所述处理器被配置为通过使在接收器Ri的所述检测信号与编码模板信号us(t)相关，从所述第一声信号和第二声信号的混合区的所述仿真的微地震现象提取所述第三信号；以及
对于所述第二多个脉冲与所述第一多个脉冲之间的多个起始时间差δ以及所述第二多个脉冲中的每个脉冲内的所述调制信号中心频率与所述第一多个脉冲中的每个对应脉冲内的所述调制信号中心频率之间的多个频率比d，重复所述模板信号与所述检测信号之间的相关，对于每个起始时间差δ和每个频率比d产生对每个接收器Ri的相关信号，表示为M(Ri,t)，它包含在根据非线性相互作用的选择规则的特定混合位置由非线性相互作用产生的仿真的微地震现象。
29.根据权利要求28的系统，其中，所述处理器被配置为使用基于压缩或剪切波速度的速度模型的轨迹分析和传播时间计算，对于每个起始时间差δ和对于每个频率比d以及所述对应的相关信号M(Ri,t)，计算仿真的微地震现象的空间坐标，其中所述两个声信号非线性地相互作用。
30.根据权利要求28的系统，其中，所述处理器被配置为应用逆时运算，使所述仿真的微地震现象的信号以传播速度模型在时间上传播回到所述微地震现象起源的时间和位置。
31.根据权利要求30的系统，其中，所述处理器被配置为确定所述微地震现象在所述起源位置的相对强度，其中在所述起源位置的相对强度与所述微地震现象的所述起源位置处所述介质的非线性性质成正比。
32.根据权利要求30的系统，其中，所述处理器被配置为通过使用基尔霍夫成像、射束成像、波动方程成像或者逆时方法，从包含所述第一声信号和第二声信号的混合区的所述仿真的微地震现象的所述相关信号来产生所述微地震现象或所述介质的非线性性质的三维图像，或者兼而有之，其中，所述相关信号从在多个接收器接收的信号获得，其中所述第一源和所述第二源位于不同位置。
33.根据权利要求30的系统，其中，所述处理器被配置为控制所述第一源和所述第二源，通过把所述第一源放置在多个位置，多次重复地由所述第一源产生所述第一信号，以及通过把所述第二源放置在多个位置，多次重复地由所述第二源产生所述第二信号；通过把接收器放置在多个位置，多次重复地接收所述第三信号；以及通过使用基尔霍夫成像、射束成像、波动方程成像或者逆时方法，从包含所述第一声信号和第二声信号的混合区的所述仿真的微地震现象的所述相关信号来产生所述介质的非线性性质的多个三维图像，其中，所述相关信号从在多个接收器接收的信号获得。
34.根据权利要求33的系统，其中，所述处理器被配置为通过速度层析成像方法更新所述传播速度模型，使按照从所述第一声源和所述第二声源的多个位置接收的信号的所述多个三维图像的差异最小。
35.根据权利要求34的系统，其中，所述处理器被配置为迭代地更新所述传播速度模型和进行成像直到获得的多个图像之间的差异最小，以及结合所述多个图像以产生所述介质的非线性性质的最终三维图像或所述局部速度比Vp/Vs或者兼而有之。