As Kevin White has mentioned, the Pierce oscillator uses just one inverting gate, and to me that feels like exactly the correct approach. Not sure about "series resonance" (I recall reading something about different resonant modes in the crystals, and crystal cutting geometries) but my general idea is, that the crystal, being a two-pin device, needs a clear AC signal applied to its pins, i.e. the two pins should be driven "inverted against each other", i.e. 180 degree phase shift at resonance. The metal case is effectively a third pin, like a loosely coupled reference ground (shielding) but that should not matter much.
I'm wondering what happens if you attach a crystal to a non-inverting amplifier. I am surprised that the circuit oscillates at all :-) I can come up with the following explanation: the two gates in series each have a "propagation delay". So the actual phase shift is not 0 ns, not 360 degrees. Depending on the logic family that you have chosen, the cascade can have anywhere between 4 ns and say 50 ns of propagation delay. Therefore, IMO your circuit resonates at some "upper harmonic mode / frequency" of the crystal, where that time-domain propagation delay of your two gates combined corresponds to a phase shift of 180 degrees (the crystal adds another 180*) and the AC gain within the loop is > 1. Actually, where the gain si >1 and highest of all such harmonic resonant modes. I am inclined to call that oscillations style "a little parasitic" :-)
Actually, my explanation above is slightly wrong. In reality, the crystal (being a resonance-heavy device) likely presents a steep gradient of phase shift around any of its resonant modes (harmonic frequencies). So it's not that "the propagation delay of the gates equals 180 degrees". Possibly not nearly that. It settles at some frequency near a particular resonant mode, where the propagation delay "adds just that wee bit needed for an optimal solution", when combined with the phase shift of the crystal being pulled up...
Thinking again, I'm probably starting to understand what Kevin might mean by series resonance, approximately... perhaps the 4th harmonic indicates a mode where the wave in the crystal tends to drive the two pins in-phase, common mode style? Against the chassis as a "reference center pin" ? But you don't have the crystal chassis grounded... This train of thought is giving me a bit of a headache :-)
The crystal is an electro-mechanical device. For optimal working conditions, it needs loading of its pins against a "common ground" (ideal center signal, maybe?) and the optimal load capacitance is mentioned in the crystal's datasheet. If you don't have one, it will likely be similar for other crystals of the same frequency and geometry. That's where the crystal supposedly resonates happily at its nominal frequency, if you drive it by 180deg-phase-shifted signal. If the "crystal load capacitance" is lower than specced, normally it will pull the crystal a little off the center frequency (within say 400 ppm, IME) - but arguably a much lower capacitive load will allow for the higher harmonic modes to prevail...
In the spring of 2018, I've done my own hobby hack with a crystal, and I have seen some harmonic oscillation too. In my case, it was the 3rd harmonic. And I had a cold joint in the circuit, which did not help debugging ;-) but abstracting from that bug, my solution to move the crystal to its base resonance was: properly load the pins (which I did right from the start) and provide an additional bit of LC filtering in the feedback loop, to make the base frequency more appealing to the loop (in terms of loop gain). My own motives were probably different from yours, but I'm providing a link because the people here have provided a lot of useful comments and help on the topic of crystals.