2D Multiline Seismic Reflection Tomography With Seismic-Tie Constraint

Abstract

A method, including: generating updated velocity models, each corresponding to one of a plurality of initial velocity models of intersecting 2D seismic survey lines, wherein updates to the plurality of initial velocity models are computed by imposing a seismic-tie regularization constraint on an inversion process that inverts for the updates to the plurality of initial velocity models, and the seismic-tie regularization constraint causes the updated velocity models to have consistent values for depth of seismic reflectivity at intersecting spatial locations.

Claims

1. A method comprising: generating updated velocity models, each corresponding to one of a plurality of initial velocity models of intersecting 2D seismic survey lines, wherein updates to the plurality of initial velocity models are computed by imposing a seismic-tie regularization constraint on an inversion process that inverts for the updates to the plurality of initial velocity models, and the seismic-tie regularization constraint causes the updated velocity models to have consistent values for depth of seismic reflectivity at intersecting spatial locations.

2. The method of claim 1, wherein a system of equations for the 2D seismic survey lines are merged diagonally, instead of vertically, to cause the updated velocity models to have consistent values for depth of seismic reflectivity at intersecting spatial locations.

3. The method of claim 2, wherein the system of equations that are merged diagonally include null matrices on an diagonal opposite to a diagonal for the system of equations for the 2D seismic survey lines.

4. The method of claim 3, wherein the system of equations that are merged diagonally include a trade-off parameter to adjust regularization of the seismic-tie regularization constraint.

5. The method of claim 1, wherein the inversion process determines the updates to the plurality of initial velocity models from residual depth error measurements from common image gathers, raytracing based on the initial velocity models, and the seismic-tie regularization constraint.

6. The method of claim 1, further comprising: from each of the plurality of initial velocity models, generating a 2D seismic reflectivity image from a 2D pre-stack depth migration; determining, for each 2D seismic reflectivity image, how much a selected reflector's depth shifts up or down along offset with respect to particular reference depth; and generating the updates to the plurality of initial velocity models with the seismic-tie regularization constraint and the reflector's depth shifts.

7. The method of claim 6, wherein the particular reference depth is a nearest offset trace.

8. The method of claim 1, wherein the inversion process determines a velocity update using a numerical solver, and the method further includes splitting the velocity update determined by the inversion process into a separate update for each of the plurality of initial velocity models.

9. The method of claim 1, wherein the plurality of intersecting 2D seismic survey lines each has a different seismic-tie regularization constraint.

10. The method of claim 1, further comprising updating at least one of the initial velocity models by adjusting a trade-off parameter included in the seismic-tie regularization constraint.

11. The method of claim 10, wherein as the trade-off parameter increases, greater weight is placed on the seismic-tie constraint.

12. The method of claim 1, further comprising using the updated velocity models to generate a subsurface image of the Earth.

13. The method of claim 1, further comprising drilling a well at a location determined at least in part from the updated velocity models and extracting hydrocarbons from the well.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0025] While the present disclosure is susceptible to various modifications and alternative forms, specific example embodiments thereof have been shown in the drawings and are herein described in detail. It should be understood, however, that the description herein of specific example embodiments is not intended to limit the disclosure to the particular forms disclosed herein, but on the contrary, this disclosure is to cover all modifications and equivalents as defined by the appended claims. It should also be understood that the drawings are not necessarily to scale, emphasis instead being placed upon clearly illustrating principles of exemplary embodiments of the present technological advancement. Moreover, certain dimensions may be exaggerated to help visually convey such principles.

[0026] FIG. 1 illustrates an example of seismic acquisition and a corresponding common image gather, generated by seismic depth migration technique.

[0027] FIG. 2 illustrates the generation of an exemplary 2D stack image from an synthetic model.

[0028] FIG. 3 illustrates a cause of the out-of-plane effect.

[0029] FIG. 4 illustrates an exemplary method of carrying out the present technological advancement.

[0030] FIG. 5 illustrates an exemplary 2D synthetic velocity model.

[0031] FIG. 6 illustrates an exemplary 2D synthetic density model.

[0032] FIG. 7 illustrates an exemplary comparison between a true model and two initial velocity models.

[0033] FIG. 8 illustrates an exemplary comparison between a true model and two updated velocity models without using the present technological advancement.

[0034] FIG. 9 illustrates an exemplary comparison between a true model and two updated velocity models using the present technological advancement.

[0035] FIG. 10 illustrates an example of a computer system useable with the present technological advancement.

DETAILED DESCRIPTION

[0036] Exemplary embodiments are described herein. However, to the extent that the following description is specific to a particular embodiment, this is intended to be for exemplary purposes only and simply provides a description of the exemplary embodiments. Accordingly, the present technological advancement is not limited to the specific embodiments described below, but rather, it includes all alternatives, modifications, and equivalents falling within the true spirit and scope of the appended claims.

[0037] This present technological advancement can adopt the approach of making a proper 2D velocity model for each line, described at the end of Background section, but resolves the seismic mis-tie problem automatically by incorporating a seismic-tie regularization constraint into tomographic inversion.

[0038] Seismic Pre-stack Depth Migration (PSDM) is a technique that depicts subsurface structures by accurately positioning seismic reflectivity between geological strata. A correct velocity of a target area is a fundamental requirement to avoid mis-positioning of reflectivity and deterioration of seismic image quality. In order to obtain a correct velocity, geophysicists often develop a smooth initial velocity model first, and then sequentially improve it with various methods and data. As an automatic method to improve an initial velocity, a seismic reflection tomography method is widely adopted, in which existing velocity is iteratively updated by inverting from seismic data misfit (Woodward et al., 2008). Data misfit is often constructed by measuring depth change of residual moveout curve (101) along offset with respect to a certain depth, typically the nearest trace (102), in the PSDM CIG gather (FIG. 1). For example, FIG. 1 is a schematic illustration that shows seismic acquisition with sources S1 and S2 and receivers R1 and R2, where waves reflect at a common midpoint between the sources and receivers. Seismic signal recorded at receivers are mapped to seismic reflectivity image for subsurface structure, in which Common Image Gather (CIG) image can be generated along different offset (101). The residual moveout curve (101) is curving up in this example; Residual Depth Error (RDE) is negatively increasing with offset, which indicates migration velocity is slower than true velocity. If velocity is correct, CIG is usually flat.

[0039] The amount of gather non-flatness is proportional to deviation of current velocity from a correct value, which can be described as,

∫.sub.raypath1/v ds=d.sub.rde   (1)

where d.sub.rde is RDE for various offsets (103, 104, 105 . . . ) between seismic sources and receivers, v is seismic velocity and s is a segment length of a ray (106 or 107) within a tomographic inversion cell (108 or 109). A ray path can be calculated by various raytracing algorithms (known to those of ordinary skill in the art) based on a given velocity model because ray path is governed by slowness (1/v) of a material, unknowns for the system of Equation (1). Since ray path is a function of velocity, the above system of Equation (1) becomes highly nonlinear.

[0040] With a use of the first order approximation for a nonlinear equation, Equation (1) can be linearized, and described in a matrix form as,

Gm=d.sub.rde   (2)

where G is raytracing kernel matrix, d.sub.rde is data vector of RDE and m is unknown vector of slowness change (1/dv). The unknown m can be iteratively updated by inverting data d.sub.rde with raytracing information. However, due to data uncertainty and nonlinearity of Equation (1), many different models m can equally fit data well within data error tolerance. Since most geological structure is smoothly varying spatially, it is more advisable to choose a velocity model which fits the data well and is spatially smooth. This can be accomplished with a regularization constraint, typically described as the first or second derivative matrix of model parameterization as below,

Rm=d.sub.rc   (3)

where R is the regularization matrix of the first and/or second order derivative combined (Zelt and Barton, 1998), d.sub.rc is data vector of regularization constraint, and determines whether a smoothness constraint is applied to a model change or an updated model itself. Equation (3) can be incorporated into a ray tracing matrix in Equation (2) as below,

[00001]$\begin{array}{cc}\left(\frac{G}{R}\right).Math.\left(\begin{array}{c}\\ m\\ \end{array}\right)=\left(\frac{d_{\mathrm{rde}}}{d_{\mathrm{rc}.Math.}}\right)& \left(4\right)\end{array}$

[0041] Then, Equation (4) can be further simplified into,

Am=d   (5)

where A is a tomography kernel, including ray path information G and smoothing regularization R, m is the unknowns of velocity change with respect to an initial velocity, and d is data vector, typically residual moveout measurement in depth, attached with the right-hand-side of a regularization constraint from Equation (3).

[0042] Assuming two 2D lines are intersecting at one common-mid-point (CMP) location, we can setup an equation for each 2D reflection tomography,

Am.sub.A=d.sub.A   (6)

Bm.sub.B=d.sub.B   (7)

where A and B are tomography kernels, m.sub.A and m.sub.B are unknowns of velocity change with respect to each initial velocity, and vector of d.sub.A and d.sub.B includes RDE in depth and the right-hand-side of smoothing regularization constraint (Equation (3)) for each 2D line.

[0043] In order to make seismic images at an intersecting point consistent, an velocity change needs to make updated velocity to be equivalent at the same spatial location (x,y,z) if geological structure is horizontally flat (e.g. there is no out-of-plane effect). This requirement can be described as below, even in the case where each line has different initial velocity, V.sub.i and V.sub.j respectively

m.sub.i+V.sub.i=m.sub.i+V.sub.j   (8)

where m.sub.i and m.sub.j are velocity change with respect to each initial model at an intersecting point. To develop a new system of equations for 2D multi-line reflection tomography with velocity-tie constraints, Equations (6) and (7) are combined diagonally as described below,

[00002]$\begin{array}{cc}\left(\begin{array}{cc}A& 0\\ 0& B\\ 0.Math.\mathrm{.Math.}.Math..Math.\lambda .Math..Math.\mathrm{.Math.}.Math..Math.0& 0.Math.\mathrm{.Math.}-\lambda .Math..Math.\mathrm{.Math.}.Math..Math.0\end{array}\right).Math.\left(\begin{array}{c}{m}_{A_1}\\ \mathrm{.Math.}\\ m_i\\ \mathrm{.Math.}\\ \underset{\_}{m_{A_N}}\\ m_{B_1}\\ \mathrm{.Math.}\\ m_j\\ \mathrm{.Math.}\\ m_{B_N}\end{array}\right)=\left(\begin{array}{c}{d}_A\\ d_B\\ \lambda \left(V_j-V_i\right)\end{array}\right)& \left(9\right)\end{array}$

where 0 is null matrix and λ is a trade-off parameter to adjust regularization of the velocity-tie constraint with respect to the original system of Equations (6) and (7). A column index of λ corresponds to a row index of velocity unknown at a certain depth located at intersecting points. In particular, all column indices of the second system of Equation (7) need to be shifted by the total number of unknown for the first system of Equation (6): in Equation (9), j+A.sub.n is column index of λ for the second system of Equation (7). If more than two 2D lines are intersecting each other, an additional system of equations can be merged into Equation (9) in the same way as constructed for two intersecting lines. A solution for Equation (9) would satisfy data misfit as well as make velocity be consistent at intersecting points.

[0044] However, when geological structure is steeply dipping with an angle of θ.sub.j along a dip direction, out-of-plane effect cannot be ignored so that equivalent velocity at intersecting points could not accomplish the same depth of a seismic reflectivity image for different lines. FIG. 2 shows three slices of a synthetic 3D velocity model (201, 202, 203), in which two seismic acquisition directions (along dip and strike) are shown by arrow (204) and nail head (205), and a corresponding 2D PSDM stack image for each direction using a correct velocity. A seismic reflector (206) along a strike direction is positioned shallower than a correct depth (207) due to out-of-plane effect.

[0045] Since seismic waves tend to follow the shortest path as possible (Snell's law), 2D velocity from 2D tomography along the strike direction should be corrected to make a same seismic reflector be positioned at a correct same depth. FIG. 3 provides a schematic illustration to show why a seismic reflector in a 2D seismic image along the strike direction (301) is positioned shallower than a correct depth (302) even with a true velocity used for migration. In order to position seismic reflectivity to the correct depth (302), overburden velocity need to be adjusted by a geometric correction parameter, 1/cos θ.

[0046] Based on an assumption that geological structures are locally continuous, 2D velocity along a dipping direction (V.sub.j) of a reflector can be related with velocity along strike direction of the reflector (V.sub.i) as below.

V.sub.i=V.sub.j/cos θ.sub.j   (10)

[0047] Thus, Equation (9) can be modified to compensate for out-of-plane effect on 2D velocity from 2D reflection tomography,

[00003]$\begin{array}{cc}\left(\begin{array}{cc}A& 0\\ 0& B\\ 0.Math.\mathrm{.Math.}.Math..Math.\frac{\lambda }{\cos .Math..Math.{\theta }_i}.Math..Math.\mathrm{.Math.}.Math..Math.0& 0.Math.\mathrm{.Math.}-\frac{\lambda }{\cos .Math..Math.{\theta }_j}.Math.\mathrm{.Math.}.Math..Math.0\end{array}\right).Math.\left(\begin{array}{c}{m}_{A_1}\\ \mathrm{.Math.}\\ m_i\\ \mathrm{.Math.}\\ \underset{\_}{m_{A_N}}\\ m_{B_1}\\ \mathrm{.Math.}\\ m_j\\ \mathrm{.Math.}\\ m_{B_N}\end{array}\right)=\left(\begin{array}{c}{d}_A\\ d_B\\ \lambda \left(\frac{V_j}{\cos .Math..Math.{\theta }_j}-\frac{V_i}{\cos .Math..Math.{\theta }_i}\right)\end{array}\right)& \left(11\right)\end{array}$

where θ.sub.j is a dipping angle of a reflector in a perpendicular direction to seismic acquisition. Compared with Equation (9), Equation (11) would not necessarily result in a same updated velocity at intersecting points, but leads to a better seismic-tie in a PSDM image.

[0048] The unknowns m are typically computed numerically by inverting a kernel and data matrix with an inversion solver. Common numerical solvers, i.e., LSQR and CGLS, efficiently handle a sparse linearized system of equation by using only non-zero elements in the matrix. In Equations (9) and (11), a seismic-tie constraint matrix is much sparser and smaller in size than kernel matrix, such that additional memory requirement and computational cost is negligible with respect to a line-by-line 2D tomography approach.

[0049] Compared with an existing approach to handle velocity mis-tie from multiple 2D line tomography, the present technological advancement can efficiently invert for velocity models for multiple 2D lines simultaneously and provide each line with a proper velocity update to tie migrated seismic traces at intersecting locations. It eliminates a need for post-tomography manual velocity editing, and ensures that seismic traces at intersecting points will be consistent even if initial velocities from intersecting lines differ.

[0050] The present technological advancement can accommodate multiple 2D reflection tomography problems simultaneously with multiple seismic-tie constraints for different lines, or can update velocity for an individual line without any tie constraint by adjusting a trade-off parameter, λ. For example, by setting λ to 0, Equation (9) and (11) provide the same solution as individual 2D reflection tomography is processed separately to fit data and smoothing regularization only. If λ increases, a solution for Equation (9) and (11) will have more emphasis on seismic-tie constraint rather than data and smoothing regularization. This flexibility can be useful where geological structures are rapidly varying laterally so that overburden velocity cannot be assumed to be continuous locally; a trade-off parameter λ needs to be reduced accordingly.

[0051] The present technological advancement, which incorporates multiple 2D seismic datasets into one system of equation, can be applicable to any 2D inversion method (e.g. 2D full waveform inversion (FWI)) or survey of different geophysical measurement (e.g. gravity, magnetic and electromagnetic data) when a material property-tie constraint at intersecting points is needed.

[0052] FIG. 4 provides an example of a method for implementing the present technological advancement. This is an example of two intersecting 2D survey lines, called by Line 1 and Line 2, respectively. More than two lines can also be implemented with the present technological advancement in the same way.

[0053] Step 10 includes developing an initial 2D velocity model for each seismic acquisition line. An initial model usually includes long wavelengths of velocity features for geological structure.

[0054] Step 20 includes preparing seismic data for a PSDM algorithm to generate offset (or angle) CIG gather by doing offset-binning and/or data regularization if necessary.

[0055] Step 30 includes running PSDM, i.e., Kirchhoff, RTM or Beam migration, using seismic data (step 20) and an initial velocity (step 10) or updated velocity from the step 130 to generate offset (or angle) gather for each 2D seismic line.

[0056] Step 40 includes, if the gather from step 30 is not sufficiently flat, proceeding to step 50. If gather is sufficiently flat, the method proceeds to step 140.

[0057] Step 50 includes estimating residual depth error (RDE) in gather for each line. Output from this step (RDE picks and their location information) serves as an input for a subsequent reflection tomography.

[0058] Step 60 includes executing 2D reflection tomography raytracing using an initial or updated velocity for each line (step 10 and 130) based on a certain parameterization of a tomographic inversion grid (see FIG. 1). A tomographic raytracing can include building files of travel time and reflection angle tables over a target update region based on the initial velocity model and assumed location of sources using any one of various known methods of seismic raytracing; tracing two rays back from each location of RDE measurement to sources and receivers on the surface, respectively.

[0059] Step 70 includes estimating dip information of geological structure using a seismic stack image, derived from stacking CIG gather (step 30) along an offset axis, or other geological information.

[0060] Step 80 includes building a raytracing kernel with measurements of ray segment length within parameterized tomographic cells (step 60) and dip information (step 70). If tomographic cells are hit by rays, write out the length of a ray segment as non-zero element of raytracing kernel, and write out column index and row index count files corresponding to non-zero elements. Write out a regularization matrix that imposes a smallest and smoothest model constraint. A data matrix from RDE for each line (step 50) is built and written out as files.

[0061] Step 90 includes estimating CIG location for each line at intersecting points, and corresponding structural dip angle perpendicular to strike direction of geologic structure.

[0062] Step 100 includes building a seismic-tie constraint matrix for Equation (8) based on information about inversion grid, CMP location of intersecting points and initial velocity models for each line; calculating row indices of inversion cells at intersecting points for each line; utilize row indices as column indices for non-zero elements of seismic-tie regularization matrix; generating a non-zero element file with a value of λ/cos θ.sub.j, where λ is a trade-off parameter to control seismic-tie constraint and θ.sub.j is a dipping angle of a reflector in a perpendicular direction to seismic acquisition; generating a row index count file by counting the number of non-zero values for a particular inversion cell (typically an index count is 2 if only a single cell for each line needs to be tied); and generating the right-hand-side matrix for seismic-tie constraint with a value of λ(V.sub.j/cos θ.sub.j−V.sub.i/cos θ.sub.i). If subsurface geologic structure is not steeply dipping around intersecting locations, write out a non-zero element file for seismic-tie regularization matrix with a value of λ, and a right-hand-side file with a value of λ(V.sub.j−V.sub.j).

[0063] Step 110 includes combining A and B diagonally by shifting all column indices of the second system of Equation (8) by the number of velocity unknowns for the first system of Equation (7). Then, merge seismic-tie constraint into a new combined system of equation to build the left-hand side matrix as described in Equation (9) and (11). Also, d.sub.A and d.sub.B (step 48) are combined, and is merged by a right-hand-side of seismic-tie constraint sequentially as described in Equations (9) and (11).

[0064] Step 120 includes inverting the combined matrix in step 110 for velocity change, (m.sub.A, m.sub.B) for both 2D lines.

[0065] Step 130 includes splitting the updated velocity change (m.sub.A, m.sub.B) between each line, and updating a current velocity model for each line by adding m.sub.A and m.sub.B, respectively, to their initial or immediately previous velocity model in order to create an updated or final velocity model (step 140).

[0066] After step 130, steps 30 and 40 are repeated by rerunning PSDM with an updated velocity model to generate new gather for each line (step 30), and measuring gather flatness (step 40). If gather is sufficiently flat, updated velocity models are saved as final velocity models. If gather is not sufficiently flat, the method reiterates through the steps until gather is sufficiently flat. A final velocity model could contain long and short wavelength of velocity features.

[0067] Although FIG. 4 shows an example of two intersecting 2D lines, this workflow can expand to more than two 2D lines using the same implementation above.

[0068] In FIG. 4, it is possible that gather for one line is sufficiently flat while the gather for the other line is not sufficiently flat. If such a situation would arise, the method would proceed as discussed above, but a velocity model that was used to generate sufficiently flat gather would not necessarily be further updated.

[0069] The updated velocity models generated by the process of FIG. 4 can be used to manage hydrocarbon production. As used herein, hydrocarbon management includes hydrocarbon extraction, hydrocarbon production, hydrocarbon exploration, identifying potential hydrocarbon resources, identifying well locations, determining well injection and/or extraction rates, identifying reservoir connectivity, acquiring, disposing of and/or abandoning hydrocarbon resources, reviewing prior hydrocarbon management decisions, and any other hydrocarbon-related acts or activities.

EXAMPLES

[0070] FIG. 5 shows a 2D synthetic velocity model (701) including a low velocity anomaly of ˜320 m/s with respect to a background velocity, linearly increasing with depth from 1600 to 2900 m/s (FIG. 5). FIG. 6 shows a 2D synthetic density model that increases stepwise from 1 to 2.7 g/cm.sup.3 at every 200 m depth spacing. Elastic modeling was carried out to generate a synthetic dataset using the true velocity model (FIG. 5) and density model (FIG. 6) assuming that two seismic acquisitions are perpendicular to each other.

[0071] FIG. 7 shows a comparison of true model (701), initial model I (702) and initial model II (703), and vertical velocity functions (704) of the initial models at CMP 479 (705). The initial model I (702) has a background velocity only without the low velocity anomaly. Initial model II (703) was created from true velocity model (701) with heavy lateral smoothing. Using two different initial velocity models for each survey line (FIG. 7), 2D Kirchhoff PSDM was run to generate CIG gathers with 75-5975 m offsets (200 m spacing) for each model. Non-flatness in the gather was measured to create RDE for each line, and this was fed into tomography raytracing to build kernels and a data matrix (Equation 9). Velocity change is determined after tomographic inversion, and then added to initial velocity model for better gather flatness. Without a seismic-tie constraint, the updated velocity models from two different initial models are different at the intersecting point (FIG. 8), although these two acquisitions are intersecting at the same spatial location. The updated models (802) and (803) are different from each other at the same spatial location.

[0072] The difference amongst (701), (802), and (803) is due to nonlinearity of reflection tomographic inversion, in which velocity change depends upon initial model and data selection. The present technological advancement, simultaneously updates the velocity models for multiple 2D lines with a seismic-tie constraint, which makes updated velocities to be more consistent at the intersecting point. Although there is still small differences between two updated models (902) and (903) at the boundary of the strong velocity anomaly (904) (FIG. 9), updated velocities becomes more close to the true velocity (701) compared with the ones without a seismic-tie constraint.

[0073] FIG. 10 is a block diagram of a computer system 2400 that can be used to execute the present techniques. A central processing unit (CPU) 2402 is coupled to system bus 2404. The CPU 2402 may be any general-purpose CPU, although other types of architectures of CPU 2402 (or other components of exemplary system 2400) may be used as long as CPU 2402 (and other components of system 2400) supports the operations as described herein. Those of ordinary skill in the art will appreciate that, while only a single CPU 2402 is shown in FIG. 10, additional CPUs may be present. Moreover, the computer system 2400 may comprise a networked, multi-processor computer system that may include a hybrid parallel CPU/GPU 2414 system. The CPU 2402 may execute the various logical instructions according to various teachings disclosed herein. For example, the CPU 2402 may execute machine-level instructions for performing processing according to the operational flow described.

[0074] The computer system 2400 may also include computer components such as nontransitory, computer-readable media. Examples of computer-readable media include a random access memory (RAM) 2406, which may be SRAM, DRAM, SDRAM, or the like. The computer system 2400 may also include additional non-transitory, computer-readable media such as a read-only memory (ROM) 2408, which may be PROM, EPROM, EEPROM, or the like. RAM 2406 and ROM 2408 hold user and system data and programs, as is known in the art. The computer system 2400 may also include an input/output (I/O) adapter 2410, 2414 GPUs, a communications adapter 2422, a user interface adapter 2424, and a display adapter 2418.

[0075] The I/O adapter 2410 may connect additional non-transitory, computer-readable media such as a storage device(s) 2412, including, for example, a hard drive, a compact disc (CD) drive, a floppy disk drive, a tape drive, and the like to computer system 2400. The storage device(s) may be used when RAM 2406 is insufficient for the memory requirements associated with storing data for operations of the present techniques. The data storage of the computer system 2400 may be used for storing information and/or other data used or generated as disclosed herein. For example, storage device(s) 2412 may be used to store configuration information or additional plug-ins in accordance with the present techniques. Further, user interface adapter 2424 couples user input devices, such as a keyboard 2428, a pointing device 2426 and/or output devices to the computer system 2400. The display adapter 2418 is driven by the CPU 2402 to control the display driver 2416 and the display on a display device 2420 to, for example, present information to the user regarding available plug-ins.

[0076] The architecture of system 2400 may be varied as desired. For example, any suitable processor-based device may be used, including without limitation personal computers, laptop computers, computer workstations, and multi-processor servers. Moreover, the present technological advancement may be implemented on application specific integrated circuits (ASICs) or very large scale integrated (VLSI) circuits. In fact, persons of ordinary skill in the art may use any number of suitable hardware structures capable of executing logical operations according to the present technological advancement. The term “processing circuit” encompasses a hardware processor (such as those found in the hardware devices noted above), ASICs, and VLSI circuits. Input data to the computer system 2400 may include various plug-ins and library files. Input data may additionally include configuration information.

[0077] The present techniques may be susceptible to various modifications and alternative forms, and the examples discussed above have been shown only by way of example. However, the present techniques are not intended to be limited to the particular examples disclosed herein. Indeed, the present techniques include all alternatives, modifications, and equivalents falling within the spirit and scope of the appended claims.