Good point. I have always said that Production PDF has to be the convolution of a Discovery PDF with some other function. Unfortunately, there is no function that when convolved with another one will give either a Logistic or Lorentzian or Gaussian for that matter. This means that the use of these curves is a contrivance that does not prove anything.

So your point on production rate falling forcing a finite convergence makes a lot of sense because it makes a lot of intuitive sense. It's like pigs sucking out of constantly shrinking trough. The pigs are furiously slurping but as the trough gets smaller, the number of pigs decreases as well.

The convolution of a Gaussian with a Gaussian gives another Gaussian. The sigmas add in quadrature. This is because the product of two Gaussians is a Gaussian and the Fourrier transform of a Gaussian is a Gaussian. It is a special function. I've derived analytically what minimal Gaussian kernels are needed to bring two arbitrarily oriented 2-D Gaussians ellipsiods of differing dimensions into coincidence with the least required distortion to each in an attempt to speed up a program I wrote which derives such kernels for galaxy images obtained with different seeing by means of minimization. The problem there is very degenerate so I wanted to give it a head start. From working on the problem, I can say that you can always take a delta function into another function but after that it gets tricky. Things that look quite narrow sometimes don't have a kernel that takes them into a broader function. But, there is always a pair of kernels that take both into a common function (each other in the trival case) and it is possible to say something about what the best pair looks like.

I think the beginning of the logistics function has a sound basis. If you are making money and you want to make more, you reinvest based on the amount of money you are making. This produces an exponential. The center of the logistics function also has a sound basis: diminishing returns. At some point, putting more effort in isn't getting you anywhere so you put your money somewhere else. I don't know why depletion would look like the logistics function. Perhaps it has to do with wanting to keep some of the people employed who made the whole thing go in the first place. If the coal seam, for example, is running out though, you can only put so many miners at the face. And if an oil well breaks down that is not producing much, you might decide not to fix it. Perhaps depletion is sentimental so it does not get a clear cut off that you might expect.

If you like, I can see if there is a kernel that takes discovery into production via convolution, but convolution is not necessarily the proper method because the time-scale of production probably depends on what has been learned from prior experience. A time varying kernel might be a better representation. It may be possible to take production curves as a function of time and derive a projection that is helpful for estimating the behavior of future or recent discoveries. One catch there though is that the recent curves which are most constraining of future curves may not be complete enough to learn from.

Chris

You are technically correct but this brings up two points:

1. The Gaussian convoluted with a Gaussian is a Gaussian only if it stretches to + and - infinity on the tails. Same thing with a Lorentzian. I see with your background in Fourier optics that this is pretty obvious, but then again spatial diffraction is different than temporal time series in that the former stretches + or - in the limit, whereas temporal has a clear one-sided forcing function due to time causality. (i.e. effects at time at -infinity cannot happen)

So this only means that it is a non-analytic Gaussian and probably approximates it only depending on how bad the temporal tail truncation is. This is even more obvious for Lorentzians convolved with each other, which have even bigger tails.

2. Yet if we may say that it approximates a Gaussian after a convolution with itself, this means that the Discovery profile is a Gaussian! Then we have to explain the rationale for this behavior. And we have another problem to deal with, as we are only switching it from the Production profile to the Discovery profile.

The temporally causal nature of the problem is why I came up with the Dispersive discovery model and the components of the Shock model which use temporally one-sided exponential functions to model the production delays and extraction rates. I am definitely not seeking kernel or windowing functions of the sort you are describing.